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SUMMARY 


The  purposes  of  this  research  program  were  to  extend  or  develop  analytical 
methods  for  determining  rotor  blade  aeroelastic  stability  limits  and  to 
perform  stability  calculations  over  a  range  of  design  and  operating  vari¬ 
ables  for  articulated  and  nonarticulated  configurations.  The  usefulness 
of  simpler  analytical  methods  is  investigated  by  comparing  results  with 
operating  boundaries  from  the  more  elaborate  analysis. 

In  the  part  of  the  investigation  presented  in  this  volume,  analytical  study 
was  made  of  the  effects  of  high  forward  speed  on  the  flutter  character¬ 
istics  of  helicopter  rotor  blades. 

Linearized  equations  of  motion  were  used  to  represent  the  dynamics  of  the 
rotor  blades.  The  aerodynamic  forces  were  obtained  by  using  a  fixed 
azimuth  approach  with  fixed-wing  two-dimensional  compressible  flow  aero¬ 
dynamic  coefficients. 

Flutter  speeds  and  flutter  frequencies  were  calculated  for  two  model  heli¬ 
copter-  blades  for  which  experimental  data  were  available.  One  blade  had 
been  tested  in  hover  at  high  tip  Mach  numbers;  the  other,  at  high  advance 
ratios  with  low  tip  Mach  numbers.  Agreement  between  theoretical  and  ex¬ 
perimental  data  was  qualitative. 

A  parametric  study  was  made  of  the  effect  on  flutter  speed  of  varying 
blade  chordwise  elastic  axis  and  center-of-gravity  position,  torsional 
and  flatwise  bending  stiffness,  and  blade  mass  ratio  and  feathering  mess 
moment  of  inertia. 

The  extended  Normal.  Mode  Transient  Analysis  was  used  to  perform  blade  re¬ 
sponse  calculations  corresponding  to  blade  configurations  and  flight  con¬ 
ditions  for  which  classical  flutter  was  calculated. 
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FOREWORD 

The  investigation  presented  .1,1  this  volume  is  part  of  a  more  general  study 
of  rotor  blade  aeroelastic  instabilities,  which  is  contained  in  five 
volumes.  The  work  was  performed  under  Contract  DA  i»l»-177-AMC-332(  )  vi 
the  U.  S.  Army  Aviation  Materiel  Laboratories,  Fort  Eustis,  Virginia. 

Mr.  Joseph  McGarvey  monitored  the  program  for  USAAVLABS. 

The  rotor  blade  classical  flutter  analysis  presented  in  this  volume  was 
developed  at  Sikorsky  Aircraft  by  Mr.  Clifford  J.  Astill,  who  also  con¬ 
ducted  the  study  of  parameter  variations  and  the  comparison  with  available 
test  data.  Mr.  Charles  F.  Niebanck,  also  of  Sikorsky  Aircraft,  conducte 
calculations  with  the  method  presented  in  Volume  I,  and  compared  the  re¬ 
sults  with  corresponding  flutter  calculations,  which  were  carried  out  with 
the  method  presented  in  this  volume. 

Volume  I  of  this  report  contains  the  development  of  the  differential  equa¬ 
tions  of  motion  of  an  elastic  rotor  blade  with  chordwise  mass  unbalance. 

Volume  III  describes  a  stall  flutter  analysis  based  on  the  calculation  of 
aerodynamic  work  for  a  cycle  of  blade  torsional  vibration.  Two-dimensional 
unsteady  airfoil  test  data  were  used  in  the  evaluation  of  the  aerodynamic 
work.  The  analysis  was  used  to  generate  stall  flutter  boundaries. 

Volume  IV  contains  the  results  of  a  study  of  static  torsional  divergence. 

A  set  of  design  charts  and  the  effects  of  a  range  of  parameter  variations 
are  presented.  The  results  of  the  static  divergence  calculation  are  com¬ 
pared  with  results  calculated  by  using  the  method  developed  in  Volume  I. 

Volume  V  presents  the  results  of  a  study  of  flapping  and  coupled  flap-lag 
instability.  The  results  of  a  parametric  study  based  on  a  single-degree- 
of-freedcm  flapping  or  flatwise  bending  analysis  are  presented.  Compari¬ 
sons  are  made  with  results  from  the  more  elaborate  analysis  developed  m 
Volume  I.  The  analysis  developed  in  Volume  I  was  used  to  deterrine  the 
coupled  flap-lag  response  of  a  rotor  to  a  number  of  sudden  control  changes. 
The  results  of  these  calculations  are  also  presented  in  Volume  V. 
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INTRODUCTION 


In  the  past,  classical  flutter  of  helicopter  rotor  blades  not  been  a 
serious  problem  because  it  has  been  customary  to  place  the  blade  elastic 
axis  and  the  center  of  gravity  on  or  near  the  25$  chord  point.  This  means 
that  the  elastic  axis,  the  axis  of  aerodynamic  centers,  and  the  center-of- 
gravity  axis  are  effectively  coincident,  thus  eliminating  coupled  bending- 
torsion  flutter  and  also  single-degree-of-freedom  flutter.  Of  course,  the 
axes  are  never  exactly  coincident;  but  for  small  deviations  from  coinci¬ 
dence,  the  application  of  suitable  flutter  criteria  is  sufficient  to  guard 
against  the  possibility  of  the  occurrence  of  classical  flutter. 

Placing  the  three  axes  in  coincidence  involves  some  penalty  in  weight  and 
increased  complexity.  To  provide  improved  structural  efficiency,  con¬ 
sideration  is  being  given  to  moving  the  axes  apart,  such  as  occurs  with 
fixed  wings.  Once  the  axes  are  no  longer  in  approximate  coincidence,  the 
blade  becomes  very  susceptible  to  classical  flutter  and  at  the  same  time 
not  very  amenable  to  a  flutter  clearance  based  on  the  application  of 
flutter  criteria.  Therefore,  flutter  calculations  have  to,be  undertaken; 
this  requires  the  use  of  an  aerodynamic  theory.  The  ulade  tip  Mach  number 
of  some  current  helicopters  approaches  unity  for  certain  flight  conditions 
so  that  compressibility  effects  are  very  important.  As  yet,  no  general 
unsteady  aerodynamic  theory  exists  for  helicopter  rotor  blades  operating 
at  high  tip  Mach  numbers  and  high  advance  ratios.  Therefore,  when  flutter 
calculations  are  made  for  these  flight  conditions,  some  approximate  aero¬ 
dynamic  theory  has  to  be  employed. 

In  addition  to  the  tendency  for  modern  helicopter  rotor  blades  to  operate 
at  high  tip  Mach  numbers  and  high  advance  ratios,  some  current  designs 
call  for  the  rotor  to  operate  in  unusual  flight  conditions.  For  example, 
the  stopped  rotor  compound  helicopter  requires  that  the  rotor  be  stopped 
at  the  transition  speed,  at  which  time  the  lift  is  transferred  to  a  pair 
of  fixed  wings.  During  the  last  few  revolutions  before  the  rotor  stops, 
each  blade  as  it  passes  through  \j/  =  270°  becomes  effectively  a  fixed 
wing  traveling  in  reverse  flow  at  the  helicopter  transition  speed,  with 
the  elastic  axis  at  75$  chord.  This  immediately  raises  the  possibility  of 
low-speed  flutter.  This  term  is  applied  to  the  single  degree  of  freedom 
classical  torsional  flutter  that  occurs  when  the  still  air  nodal  line  of 
a  fundamental  mode  lies  on  or  near  the  75$  chord  point.  With  the  nodal 
line  in  this  position,  it  is  possible  for  flutter  to  occur  at  very  low- 
airspeeds.  In  the  past,  this  type  of  flutter  was  mainly  of  academic  in¬ 
terest  and  consequently  theoretical  and  experimental  studies  cf  this 
phenomenon  are  noticeably  absent  from  the  literature.  However,  with  the 
advent  of  stopped-rotor  helicopters ,  low-speed  flutter  will  have  to  be 
given  serious  consideration.  The  methods  developed  in  this  study  are 
directly  applicable  for  this  condition. 

Thus,  for  current  and  future  helicopter  rotor  blades,  classical  flutter 
calculations  begin  to  acquire  the  same  importance  as  they  have  in  i.he  past 
for  fixed  wings.  The  present  study  is  part  of  the  current  effort  tho.t  is 
being  made  in  this  direction.  The  study  is  strongly  oriented  toward  pro- 


viding  a  method  for  rapidly  obtaining  estimates  of  rotating  blade  classical 
flutter  speed  and  the  effect  on  the  flutter  speed  of  varying  parameters 
such  as  blade  chordwise  elastic  axis  position.  Such  a  method  is  mandatory 
for  design  use,  where  introduction  of  new  blade  concepts  requires  that  the 
classical  flutter  solutions  be  attempted  with  the  best  available  method  of 
analysis.  Consistent  with  the  approach,  the  fixed  azimuth  method  was 
chosen  in  conjunction  with  the  use  of  a  digital  computer  to  obtain  solu¬ 
tions  to  the  equations  of  motion.  This  is  in  contrast  to  some  recent 
studies  which  have  used  an  analogue  computer  to  integrate  the  equations  of 
motion  with  time  varying  coefficients.  The  fixed  azimuth  approach  is  an 
expedient  way  of  converting  the  time  varying  coefficients  to  fixed  coef¬ 
ficients  which  are  more  convenient  for  digital  computation.  It  is  probably 
the  major  assumption  made  in  deriving  the  equations  of  motion.  With  a 
rotating  blade,  the  first  two  aerodynamic  effects  listed  below  result  in 
time  varying  coefficients  in  the  equations  of  motion. 

1.  The  dynamic  pressure  varying  azimuthally  for  fi>0 

2.  The  indicial  aerodynamic  response  due  to  the  varying 
of  the  velocity  field  azimuthally  for  ul  >  0 

3.  There  is  a  phase  lag  between  the  occurrence  of  maximum 
angi'.  of  incidence  and  the  resulting  development  of  maximum 

‘  lif  e  when  the  blade  executes  simple  harmonic  oscillations . 

This  phase  lag  is  the  result  of  indicial  aerodynamic  response 
to  the  varying  geometric  angle  of  attack;  for  simple  harmonic 
oscillations,  it  is  integrated  to  give  an  in-phase  and 
quadrature  response. 

4.  The  aerodynamic  forces  are  determined  by  multiple 
nonplanar  shed  wakes. 

The  fixed  azimuth  approach  considers  a  particular  azimuth  angle  only.  The 
instantaneous  velocity  of  the  blade  at  this  azimuth  angle  is  considered  to 
apply  for  all  time,  that  is,  -CO  <  t  <  00  ,  and  the  aerodynamic 

forces  are  calculated  in  accordance  with  the  assumption.  Hence,  only  the 
third  effect  is  taken  into  account,  the  lag  effects  being  derived  from 
fixed-wing  theory.  The  neglect  of  the  first  two  effects  causes  the  coef¬ 
ficients  in  the  equations,  of  motion  to  be  time  invariant.  The  analogue 
computer  studies  referred  to  above  considered  only  the  first  effect,  so 
that  steady-state  aerodynamics  were  used.  So  while  the  coefficients  in  the 
equations  of  motion  were  time  varying,  all  but  the  first  effect  was  neg¬ 
lected.  In  this  respect  it  should  be  noted  that  when  steady-state  aero¬ 
dynamic  coefficients  are  used,  it  is  not  possible  to  predict  single-degree- 
of-freedom  classical  flutter  which  depends  on  unsteady  aerodynamic  lag 
effects . 

The  fixed  azimuth  classical  flutter  study  begins  with  the  derivation  of  the 
equations  of  motion,  and  the  presentation  of  the  set  of  assumptions  and 
approximations  upon  which  the  equations  are  based.  These  equations  of 
motion  are  used  to  calculate  the  flutter  speeds  of  sets  of  model  rotor 
blades  for  which  experimencal  *  ..suits  are  available  for  correlation.  Then 


2 


the  method  described  in  this  volume  is  applied  to  the  investigation  of  the 
change  in  flutter  speed  due  to  a  change  in  various  parameters.  Finally, 
the  results  of  a  number  of  calculations  made  with  the  more  elaborate 
Extended  Normal  Mode  Transient  Analysis  described  in  Volume  I  are  presented 
and  compared  with  the  corresponding  fixed  azimuth  classical  flutter  calcu¬ 
lations  . 
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ASSUMPTIONS 


The  equations  of  motion  are  linearized  by  assuming  that  the 
linear  and  singular  displacements  of  a  blade  element  are  small. 

Strip  theory  is  used  to  compute  the  aerodynamic  forces.  This 
implies  that  the  pressure  distribution  on  any  chordwise  strip 
of  a  blade  is  the  same  as  if  that  strip  were  part  of  a  two- 
dimensional  blade  which  performed  the  same  motion  everywhere 
as  the  strip  under  consideration.  Thus,  the  interaction  be¬ 
tween  strips  is  not  accounted  for. 

The  aerodynamic  forces  are  calculated  using  a  fixed  azimuth 
approach.  This  method  considers  a  particular  azimuth  angle 
only.  The  instantaneous  velocity  of  the  blade  at  this  azimuth 
angle  is  considered  to  apply  for  all  time,  that  is,-C0<t  <00 
The  aerodynamic  forces  are  calculated  in  accordance  with  this 
assumption;  as  a  result,  the  coefficients  in  the  equations  of 
motion  become  time  invariant.  Also,  the  indicial  aerodynamic 
response  associated  with  the  changing  velocity  field  is  neg¬ 
lected. 

The  assumed  two-dimensional  oscillatory  aerodynamic  forces 
acting  on  a  strip  are  approximated  by  two-dimensional  fixed-wing 
theory.  This  means  that  the  following  assumptions  are  made: 

a.  The  steady  angle  of  incidence  is  small  and, 
hence,  well  below  the  stall  region. 

b.  Effects  due  to  a  nonplanar  wake  and  inter¬ 
action  between  shed  wakes  are  neglected. 

The  contribution  to  the  generalized  stiffness  of  bending  modes 
due  to  radial  flow  is  neglected.  This  means  that  the  solutions 
in  the  region  of  \Jr  =  0°  and  180°  become  increasingly  less 
accurate  as  advance  ratio  increases. 

Zero  steady  state: 

a.  Blade  twist. 

b.  Blade  pitch  angle. 

c.  Blade  lag  angle. 


EQUATIONS  OF  MOTION 


GENERAL  DISCUSSION 

The  form  in  which  the  equations  are  cast  for  eigenvalue  solution  is  in¬ 
timately  tied  up  with  the  notation  used  to  represent  the  aerodynamic 
forces.  Measuring  the  blade  displacements  and  the  aerodynamic  forces 
at  the  25$  chord,  then  for  simple  harmonic  motion,  the  lift  L  and  moment 
M  acting  on  a  blade  segment  of  width  Ar  of  a  nonrotating  blade  in  in¬ 
compressible  two  dimensional  flow,  are  given,  as  shown  in  Figure  1,  by  the 
following  equation: 


1 

Tfp^Er 


C/4 


I  +  2C(k) 


(1) 


where 


C(k)=  F(k)  +  iG(k)  is  the  Theodorsen  circulation  function 

k  =  u> b/V  is  the  reduced  frequency  parameter 

A  complex  polar  plot  of  C(k)  is  shown  in  Figure  2.  The  following 
points  should  be  noted: 

1.  The  acceleration  (inertia)  coefficients  are  constant;  thus, 
the  aerodynamic  inertias  are  indistinguishable  from  structural 
inertias . 

2.  The  velocity  (damping)  and  displacement  (stiffness)  coefficients 
are  functions  of  the  complex  circulation  function  C(k).  They 
vary  with  reduced  frequency  k,  but  remain  finite  as  k  ranges 
from  zero  (steady  flow)  to  infinity  (zero  forward  velocity). 
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When  these  coefficients  are  incorporated  into  the  flutter  equations,  it  is 
customary  to  modify  them  with  the  relations  for  simple  harmonic  motion; 
namely , 

q  =  icaq 
q  =  -w2q 

i=^T  (2) 


where  q  is  any  displacement,  such  as  a  or  h 

Two  methods  are  widely  used  and  result  in  different  methods  of  presenting 
the  equations  of  motion.  They  are  the  "V-Q  method"  and  the  "decay-rate 
method" . 

The  V-Q  method  is  used  i?  the  United  States.  A  discussion  of  this  method 
appears  for  example,  on  page  2^1  of  Reference  1.  The  aerodynamic  inertia, 
damping  and  displacement  coefficient  matrices  are  combined  to  form  r.  single 
aerodynamic  complex  stiffness  matrix.  The  above  equation  (l)  becomes 


so  that  the  aerodynamic  coefficients  become  infinite  for  the  case  of  steady 
flow.  This  is  due  to  the  algebraic  rearranging  of  the  terms  and  does  not 
imply  that  the  total  aerodynamic  forces  become  infinite.  The  coefficients 
remain  finite  for  the  case  of  zero  forward  speea  (K=CO) 
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With  the  aerodynamic  forces  expressed  in  the  above  manner,  the  equations 
of  moticn  can  be  put  into  the  following  form 


{[a]  -  fe)2<m9>[Bu]}{q}  ■  {0} 


For  a  particular  value  of  the  reduced  frequency  k,  the  coefficient  matrix 
[Ay]  can  be  evaluated,  and  the  eigenvalue  yields  the  flutter  frequency 
oj  aad  the  value  of  the  structural  damping  coefficient  g  (assumed  here 
to  be  the  same  for  each  mode)  required  to  keep  the  system  oscillating  with 
simple  harmonic  motion.  The  forward  speed  V  is  found  using  and  the 
value  of  k  assumed  when  calculating  the  coefficient  matrix  [Ay]  . 

That  is , 


(6) 


The  forward  speed  will  be  different  for  each  mode.  When  the  process  is 
repeated,  assuming  different  values  for  k  gives  the  familiar  V-g  plot,  as 
in  Figure  3. 


The  decay-rate  method  is  used  in  Great  Britain.  The  complex  aerodynamic 
forces  are  rearranged  into  real  coefficient  matrices ,  with  the  relations 
for  simple  harmonic  motion  being  used,  so  that  the  equations  become,  with 
forces  and  displacements  referred  to  the  leading  edge. 
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again  it  should  he  noted  that 


lim 

k-0 


(8) 


Some  of  the  coefficients  become  infinite  for  both  steady  flow  (k-0) 
and  zero  forward  speed  ( k  =  CO )  , 

With  the  aerodynamic  coefficients  expressed  in  the  above  manner,  the 
equations  of  motion  can  be  put  into  the  following  form: 


[A]T{ii}+  (V[S]AEH0+  [BR]  ){q}  +  (V2[C]#ER0  +  [E]){,}  *  {0}  (9) 


V  and  give  the  effect  of  dynamic  pressure  so  that  the  variation  of  the 
aerodynamic  matrices  [b]aero  and  [C]AER0  with  k  is  due  only  to  the 
unsteadiness  of  the  flow.  This  is  in  contrast  to  the  V-g  method  where  the 
variation  of  the  complex  matrix  [A0]  with  k  includes  the  effect  of  dynamic 
pressure.  This  occurs  because  the  dynamic  pressure  is  expressed  in  terms 
of  the  reduced  frequency  k  according  to 


In  the  V-g  method,  the  oj  is  divided  out  and  appears  only  in  the  com¬ 
plex  eigenvalue 

(1 ir)2  ( I  +  >9 )  (11) 

In  the  decay-rate  method,  the  equations  are  solved  by  assuming  a  value  for 
k  which  allows  all  the  matrices  to  be  computed.  That  is,  the  degree  of 
unsteadiness  of  the  flow  at  the  flutter  condition  is  assumed  in  advance. 
Then  for  a  particular  forward  speed  V,  the  total  damping  matrix  V  [b]aer0 
+  [Br]  and  the  total  stiffness  matrix  V2  [c]AER0  +  [E]  are  formed,  and 
the  eigenvalue  solution  is  obtained.  This  process  is  repeated  for  a  number 
of  values  for  V  and  but  using  always  the  same  matrices  [BjAER0  and 
[c]aero  and,  hence ,the  same  assumed  value  for  k.  This  gives  a  plot  of  the 
exponential  decay  rate  C/Cc  for  each  of  the  modes,  as  in  Figure  U.  The 
decay  rate  is  nondimensionalized  with  respect  to  the  critical  damping 
rate  in  the  usual  manner. 


When  using  the  decay-rate  method,  the  convention  of  plotting  C/Cc  with  the 
stable  region  above  the  abscissa  is  usually  followed.  It  is  of  interest 
to  examine  the  relationship  between  the  two  methods.  The  two  methods  will 
yield  the  same  flutter  speed  and  frequency.  However,  the  variation  of  g 
versus  V  will  be  different  from  the  variation  of  C/Cq  versus  V  even 
though  they  cross  the  abscissa  at  the  same  point.  Sometimes  an  attempt  is 
made  to  convert  C/C c  to  C,  using  the  relation 

Q  =  2("fe)  0.2) 

This  is  actually  a  contraction  of  the  equation 

9  =  2  (feX^)  (13) 

which  is  exact  for  a  single-degree-of-freedom  system.  However,  the  above 
conversion  has  no  justification  for  mult i-degree-of- freedom  systems. 

The  question  arose  as  to  which  of  the  above  methods  to  use  in  the  present 
study.  The  velocity  distribution  across  the  blade  varies  with  forward 
speed  for  a  constant  rotor  speed.  If  the  helicopter  forward  velocity 
is  used  in  place  of  V,  the  above  methods  may  be  applied  to  the  rotating 
blade  in  the  same  way  that  they  are  applied  to  calculating  fixed-wing 
flutter  speeds.  However,  the  aerodynamic  matrices  are  now  a  function  of 
the  advance  ratio  p ,  which  has  to  be  known  in  advance  in  order  to  be  able 
to  calculate  the  matrices. 

If  the  V-g  method  is  used,  it  is  necessary  to  assume  a  value  for 
p-  VH/iIR  ;  and  corresponding  to  the  value  of  k  =  cub/VH  chosen,  the 
flutter  solution  gives  a  value  of  for  each  mode.  The  corresponding 

value  of  rotor  speed  is  found  from  the  equation 

=  TZr  (1U) 

Thus,  the  rotor  speed  Cl  will  be  different  for  each  mode.  In  obtaining  a 
Vh“9  plot,  it  is  required  either  that  the  rotor  speed  be  constant  or  that 
it  vary  in  some  prescribed  manner  with  .  In  either  case,  it  is  most 

unlikely  that  Cl  for  the  mode  of  interest,  calculated  using  the  above 
relationship,  will  follow  the  required  variation.  As  it  is  necessary  to 
have  Cl  very  close  to  its  correct  value,  an  iteration  procedure  is  re¬ 
quired  . 

When  applying  the  decay-rate  method,  the  forward  velocity  is  stated;  with 
Cl  known,  then  p  i'  known  in  advance.  Only  k  has  to  be  assumed,  so  that 
the  method  of  solution  is  basically  the  same  for  both  fixed  wings  and  ro¬ 
tating  blades,  because  the  solutions  are  obtained  for  specified  values  cf 
V  and  k,  then  from  the  relation 


9 


(15) 


this,  in  effect,  means  that  the  solution  frequency  ai  has  been  assumed  in 
order  to  obtain  a  solution.  It  must  be  clearly  understood  that  this  im¬ 
plies  only  that  the  degree  of  unsteadiness  of  the  flow  is  assumed  in  ad¬ 
vance.  For  flutter  not  involving  aerodynamic  control  surfaces,  as  is  the 
case  with  the  present  study,  fixed-wing  flutter  solutions  indicate  that 
the  solution  frequency  and  decay  rate  is  relatively  insensitive  to  the 
assumed  degree  of  unsteadiness  of  the  flow.  Hence,  a  reasonable  guess  at 
the  flutter  frequency  gives  solutions  sufficiently  accurate  for  most 
flutter  investigations,  especially  when  the  approximate  location  only  of 
stability  boundaries  is  required,  and  this  comprises  by  far  the  major  part 
of  most  flutter  investigations.  Should  a  very  accurate  value  of  flutter 
speed  be  required,  or,  alternatively,  if  the  assumed  flutter  frequency 
turns  out  to  be  very  inaccurate,  then  one  iteration  using  the  calculated 
flutter  frequency  is  all  that  is  needed  to  obtain  satisfactory  correspondence 
between  the  calculated  and  ^he  assumed  flutter  frequency. 

It  was  decided  that  the  decay-rate  method  of  solution  was  more  direct  than 
the  V-g  method  for  analyzing  rotor  blade  flutter,  and  so  this  method  was 
adopted.  One  modification  to  the  decay-rate  method  is  necessary  when  it 
is  applied  to  rotating  blades ,  and  that  is  to  recalculate  the  aerodynamic 
matrices  for  each  forward  velocity,  as  they  are  now  a  function  of  advance 
ratio. 

Having  chosen  to  use  the  decay-rate  method,  it  was  then  necessary  to  choose 
a  definition  of  the  aerodynamic  coefficients.  Various  notations  are  used 
in  presenting  the  aerodynamic  forces,  but  none  of  these  is  entirely  suit¬ 
able  for  the  present  study.  Some  take  the  factor  i r  outside  the  equation 
so  that  the  stea  y-flow  rate  of  change  of  lift  coefficient  with  angle  of 
attack  becomes  equal  to  2.  Others  take  out  the  factor  2  so  that  Clq  =  ir 
at  k  =  0.  For  convenience,  it  was  decided  to  take  out  a  factor  of  2tt  , 

giving  Clq  =  1.0  at  k  =  0  so  that  this  term  then  becomes  a  measure  of  the 
aerodynamic  efficiency.  The  aerodynamic  coefficients  are  defined  as 
follows  again  with  reference  to  Figure  1: 


(16) 
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Here,  the  aerodynamic  acceleration  (inertia)  terms  have  been  combined  with 
the  displacement  of  coefficients  with  the  following  relation  being  used: 


(IT) 


Although  for  incompressible  flow  any  acceleration  coefficient  Dq  is  a 
constant,  the  equivalent  displacement  coefficient  k2Djj  varies  with- k 
and  becomes  infinite  for  zeio  forward  speed  (k=00)  .  When  the  decay-rate 
method  is  used,  it  is  desirable  that  the  variation  of  the  coefficient  with 
change  in  k  be  a  minimum  in  order  to  keep  the  calculated  flutter  speed 
relatively  insensitive  to  the  value  of  k  assumed  in  obtaining  the  solution. 
Hence,  it  would  be  desirable  to  keep  the  acceleration  coefficients  separate, 
but  the  following  two  factors  weighed  against  this : 

1.  In  compressible  flow,  the  acceleration  coefficients  become 
complex  functions  of  the  reduced  frequency  k. 

2.  Most  tables  of  compressible  coefficients  have  already  lumped 
the  acceleration  coefficients  in  with  the  displacement 
coefficients. 

For  incompressible  flow,  the  above  coefficients  have  the  following  form: 


4  =  F(k) 

(18) 

4  =  i  +  F(k)  +  -Sto) 

(19) 

4  =  ■  k  (G(k)  +  -jj-) 

(20) 

4=  F(k)  -  k(G(k)+  \) 

(21) 

-mfi  =  0 

(22) 

-md=  ± 

(23) 

kz 

~mh  =  “T 

(21.) 

-ma  =  ~  ||k2 

(25) 

ll 

It  should  be  noted  that 


lim  la  =  -co 

k— -0 


(<*) 


and  that 


lim  /h  = 
k— >CO 

lim  la  -  -CO 
k-CO 


(27) 

(28) 


»m  (-mu)  =  -CO 

(29) 

lim  (-mn)  =  -  CO 
-CO' 

(30) 

KINETIC  ENERGY 


Consider  a  set  of  axes  (  X  ,  Y  ,  Z  )  fixed  in  the  helicopter  with  the  Z 
axis  coinciding  with  the  rotor  shaft  axis  and  the  X  axis  pointing  aft, 
as  shown  in  Figure  5-  Then  consider  another  set  of  axes  (  X  ,  Y  ,  Z  )  ro¬ 
tating  with  the  Y  axis  coincident  with  the  blade  under  consideration.  The 
blade  is  assumed  to  have  zero  twist,  zero  steady  pitch  angle  with  the  tip 
path  plane  lying  in  the  XY  plane. 

The  coordinates  of  any  point  (x,y)  on  the  blade  are  given  by 


X  1  r  cos;ilt  -siniifirx* 

yj  Lsin.dt  cos&tj  LyJ  (31) 


rim 


The  kinetic  energy  T  of  the  blade  is  given  by 


where  dm  is  an  elementary  mass  point  at  (  x  ,  y  ,  z  )  on  the  blade. 
T  =  j  J  [x2+  y2+z2+  2 XI  ( xy  -  yx)  +  &2(x2+/)]  dm 


Generalized  coordinates  .  .  .  qn  may  be  used  to  describe  the  displace¬ 
ment  of  the  blade  relative  to  the  rotating  axes.  That  is, 

X  *  X  (q,,  -  -  -  ,qn>  (31, 

y  =  V  (q..--  -,q„)  (35 

z  =  z  (q,  , - ,qn  )  (36 

If  the  above  terms  are  expanded  in  terms  of  a  set  of  Maclaurin's  series, 
the  following  result  is  obtained. 


♦?(&).<*  +  *??(*&)«,«, 


where  x0  ,  y0  and  z0  are  the  values  of  x  ,  y  and  z  when  each  q.^ 
is  equal  to  zero.  The  derivatives  are  also  evaluated  when  q^  equals  zero. 

Substituting  into  the  expression  for  the  kinetic  energy  and  neglecting 
higher  order  terms  gives 

t  •  */■  {?  ?  [(fti  utx  *  «a«&  *  («. »,).  ] 

+  2fl  ?  ^  Qi’i  (toy).  ] 

t  n'  £ 

+  2rf  ,kc 
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V 


This  is 


T  a  \  Z  ZGMij  *ii  Qj  +  SL  Z  Z  Gjj  qj  qj  +  ^  ZZ^ij  qj  qj 


I  1 


i  ] 


'  J 


+  a  Z  *■</ >(aqXi  1 +  y°( aqYi  1] dm 


BLADE 


where 
G 


Mij  =  /[(a5il(lq))o  +  (3q|)o(lq1)o+(l|l(a5:l]  dm 


8LA0E 


a->  =/[(^l-(fU(^i]  - 


(ill) 


(U2) 


(1»3) 


BLADE 


r"  =./I  W)°+  )• +  x°('aq7^)= +  ) dm 


Note  that 


GMij  =GMji 


(M) 

(•♦5) 


Gjj  =  — Gjj 

(1*6) 

Tij  =  TJi 

POTENTIAL  ENERGY 

(UT) 

The  potential  energy  arises  from  the  elastic  deflection  of  the  blade  and 
control  system  and  from  +he  motion  of  the  blade  in  a  gravitational  field. 

It  is  a  function  of  the  generalized  coordinates  q^  and  shown  as  follows: 

U  =  U  (a, , - ,qn) 

(U8) 

i 

;( 


If  Maclaurin's  expansion  is  used,  Eq.  ( 1+8 )  may  be  expanded  to  give 


II 

i  j 


_aiu_ 

aqi  aqj 


qiqj 


(1*9) 


1U 


where 


U0  =  U(0, - ,0) 


(50) 


(51) 


and  where  higher  order  terms  have  been  neglected.  The  gravitational 
forces  have  a  negligible  influence  on  the  perturbation  deflections  q^,  and 
so  it  is  assumed  that  the  generalized  spring  forces 

K..  =  (  d2U  ) 

IJ  \dqjdQj /o  (52) 


do  not  depend  on  gravitational  forces  . 
Hence, 


u  =  Uo  +  Y,  (aq;)o^i  +  2  £  X'  **ij  qi  qj  (53) 

i  1  J 


LAGRANGE  S  EQUATION  OF  MOTION 

The  equations  of  motion  in  the  form  of  Lagrange's  Equation  will  be  used 
which,  for  the  r  th  generalized  coordinate  qr,  may  be  written 


d. 

dt 


di  ,  au 
dqr  <3qr 


Qr 


r=  I, - ,n 


(5»0 


where  Qr  accounts  for  all  forces  not  included  in  T  and  U.  It  may  be 
separated  into  an  initial  steady  component  (Qr)o  and  a  perturbation 
component  AQr  .  Substituting  the  expressions  for  T  and  U  into  the  above 
equrtion  gives 


I  [GMrj  qj  +  fl(Gjr  -  GrJ)  q(  +  ( Krj  -  B2Trj)  q(  ] 

+  {*X  ~  dm  '  (Qr)°  +  AOr 


BLADE 


(55) 
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The  steady  generalized  force  (Qr)o  must  be  in  equilibrium  with  the 
initial  elastic,  gravitational,  and  centrifugal  forces;  that  is. 


“  flZ/  [Xo('5qr)0t  y°('^)o  ]  dm  °  (Qr)c 


Hence,  the  equation  of  motion  reduces  to 


(57) 

The  generalized  force  AQr  is  due  to 

1.  Aerodynamic  forces  acting  on  the  blade 

2.  Mechanical  damping  forces 

They  can  be  determined  from  the  work  done  by  these  forces  in  a  virtual 
displacement  using  the  relationship 

^  (58) 

The  general  form  of  the  equations  of  motion  having  been  obtained,  a  set  of 
generalized  coordinates  can  be  chosen  and  the  general  expressions  can  be 
evaluated. 

A  new  set  of  axis  (  X ,  ,  Yj  ,Z|  )  are  introduced,  which  are  fixed  in  the 
blade  as  shown  in  Figure  6.  The  origin  0^  coincides  with  the  blade  elastic 
axis  at  the  section  under  consideration.  With  each  q^  equal  to  zero, 

0-jX^  coincides  with  OX 

OjY-^  is  parallel  to  0Y 

OjZ^  is  parallel  OZ 

hz  is  the  flatwise  displacement  of  the  elastic  axis 

hx  is  the  edgewise  displacement  of  the  elastic  axis 

a  ia  the  blade  pitch  displacement 

€  is  the  distance  of  the  blade  raid-chord  aft  of  the 

axis  at  the  section  under  consideration 


1 6 
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Substituting  for  Zp  and  Z 
the  modal  deflections  gives 


and  neglecting  third-order  terms  in 


The  third  step  considers  apparent  blade  shortening,  A3 
Let  s  be  the  arc  length  measured  along  the  elastic  axis 


ds  =  dyA  /  + 

V  Vdy /  Vdy J 


With  expansion  by  the  t  nomial  theorem  and  with  higher  order  terms  being 
neglected,  the  total  blade  shortening  at  station  y  is  given  by 

A,  =  /  (ds-dy)  (7 

Jo 

Hence,  the  y  coordinate  of  P  referred  to  the  XYZ  axes  is  given  by 
yp  =  r  -  (A,  +A2  +  A3)  (7 

C 

That  is,  the  coordinates  of  P  referred  to  XYZ  axes  are  given  by 

2 

xP  =  (€ +aEAb +x, ) +hz +z(a  -  xc^-  (■ 

yp  =  r  -  -  [z.  -  x,a]  (■$“•)  -  i P  [  +  fe) 


The  displacements  of  the  point  P  can  now  he  expressed  as  being  due  to 

two  generalized  coordinates  q,  and  q2  .  This  will  give  the  form 

of  the  coupling  terms  between  the  two  generalized  coordinates  and  will 
allow  the  final  equations  of  motion  to  be  generalized  by  inspection  to  in¬ 
clude  any  number  of  generalized  coordinates. 


Put 


hx  (y,t) 

=  g*(y)q,(t)  +  g2y  g2(t) 

(78) 

hz  (y,t) 

=  gf(y)  q,0)  +  g2(y) 

•  (79) 

a(y,t) 

=  F,  (y)  q,(t)  +  F2  (y)  qp) 

(80) 

where  q*(y)  )  q*(y)  and  F|(y)  are  respectively  the  edge¬ 

wise,  flatwise,  and  torsional  mode  shapes  of  the  ith  generalized  coordin¬ 
ate.  Substituting  into  the  expressions  for  xp  ,  yp  and  zp  and 
neglecting  third  and  higher  order  products  of  the  generalized  coordinates 
gives 


xp  =  («+aEnb+X|)  +  (g;  +  (g*  +  z,F2)q2 

2xiFi^i  ~  xiFi  ^QiPz  ”  x|p2  (8l) 
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zp  =  z,  -  ( gf  -  X,F,)  q,  +  (92  “  x,  Fz)  q2 

_?ZlFiq'  ~ZlFlF2Cl|C?2  ""2  2lF2^2  (83) 

Differentiating  the  above  expressions  with  respect  to  q,  and  q2  and 
evaluating  the  terms  at  q,  =  q2  =  0  gives  the  following  results: 

(xP)0  =  €  +  aEAb  +  x,  (81) 
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mxs  =  mx  +  m(es  +  aEAb) 

( 10U) 

mZ  =  ffp bzi  dx>dzi 

( 105) 

Ixx  =  ffpbf  dx«dz, 

(106) 

Xxz  =  X»ZI  dx,  dz, 

(107) 

izz  =  dx«  dz« 

(108) 

Iyy  =  ffpb^Z  f+xf)dx,dz, 

(109) 

!xx+  -^zz 

where  pb  =  blade  structural  density  and  the  integration  JJ  is  over 
the  blade  cross  section  . 

The  following  relation  will  also  be  used: 


j*Rf(x,  ,  z,  ,  r)dm  =  ///f(x, , z,  ,r)/>b  dx,  dz,  dr  (lio) 

0 

JQ  f(x,,z,  ,  r)dm  =  f  [^b  f(x,,z,  ,r)dx,  dz,]dr  (in) 


where  r  is  now  used  instead  of  y  to  refer  to  the  blade  spanwise 
coordinate. 

GENERALIZED  INERTIAS 


Since  the  generalized  inertxa  is  being  developed  with  respect  to  q,  and 
q2  ,  i  =  I  and  j  =  2  in  this  case. 
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0 

0 

0 

9  \ 

0 

0 

-ml 

-ml 

0 

0 

9*, 

0 

ml 

0 

0 

1,, 

0 

1(95) 

0 

mJ 

0 

0 

I « 

0 

£(«;) 

0 

0 

-l*x 

-Li 

0 

0 

* 

0 

0 

0 

0 

0 

0 

dr 


(117) 


GENERALIZED  CENTRIFUGAL  STIFFNESS  COEFFICIENTS 

The  subscripts  of  Tjj  are  i=l  and  j  =  2  ,  as  for  the  GMt2 

G 12  mass  properties  developed  above  in  Eq.  (112)  through  (117). 


V  T<2  1  I*  [©o('§)o+  (4j)o  (-S~)o  +  (xp)o('5q^  l  +(yp)o(lq^)o]  dm 

rij  =  T«2  **  f0\<<  m  +[{(^J’)fi+(^)F.}r  -<^o£At>)f-F2]  mx 

+  [  (af  ^ +  £ f.)-  { ('S1')  Fz  -1-  (^) F] r  ] ™  +  { (%*)(%■)  « 


and 


(118) 


\  du  / 

(M.) 

\  du  / 

(£)] 

\  du^  dr 

(119) 

It  will  be  seen  that  the  mode  shapes  1 

r£gT^ 

^  dr  / 

and 

(#) 

occur  inside  an  integral  with  variable  upper  limit.  This  complicates 
numerical  calculations ,  and  so  the  mode  shapes  are  brought  outside  the 
integral  by  changing  the  order  of  integration  as  follows: 


/•R  r  /*r  1  /*R  r  /*R  1 

/  0(v)  /  <£(u)du  dv  =  /  4>(u)  /  0(v)dv 

Jo  l  Jo  J  •'o  *-  Jr  J 


du 


2l» 


(120) 


so  that 


(121) 


(122) 


(123) 


0  m 
0  0 
0  0 
0  mz 
0  0 


0 

0 

0 

0 

"gi 

0 

0 

mz 

0 

q  i 

rmx 

0 

r 

-rmz 

0 

rmx 

-rmz  [I„ 

“  lzz~(e+aEAb)mx] 

0 

h 

0 

0 

0  J 

0 

(12U) 
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ARTICULATED  RIGID  BLADE  NATURAL  FREQUENCIES 


Before  proceeding  any  further,  it  is  of  interest  to  obtain  the  frequencies 
of  the  rigid  body  flapping,  lagging,  and  pitching  modes  due  to  the  centri¬ 
fugal  force  field.  The  equation  of  motion  for  a  single  degree  of  freedom 
q |  is  given  by 


Mu  q,  -  T|t  q  |  =0 


„  p 

and  for  simple  harmonic  motion  q.  =  -w  q. 


OJ  _  /-T|| 

n  V  M„ 


where 


so  that 


(125) 


(126) 


r  mx  -  ( «  +  OeaM  F.2  mx 


T. 

+  ■»«  -  2 fir)  +  {(t))‘  -  F.1} I.. 


(127) 


cJMi  =  m  rdr 


(128) 


f0  [{faff  +  (gf)2  }  m  -  F,  m*  +  2g*  F,  mz 


(129) 
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Rigid  Blade  Flapping,  Frequency 


For  a  flapping  mode 


,z  . 


0 

r-e 

R-e 


0  <  r  <  e 
e  <  r  <  R 


(130) 

(131) 


dgf 

dr 


0 


R-e 


0  <  r  <  e 
e  <  r  <  R 


g 


X 

I 


sr 


F,  =  0 


where  the  flapping  hinge  is  at  r  =  e.  In  addition,  Ixx  =  0  . 
Then, 


where 


(132) 

(133) 

( 13U) 

(135) 

(136) 

(137) 

(138) 

(139) 
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For  a  lagging  mode. 


9*  ■ 

0 

O  <  r  <  e 

(U*o) 

■ 

r-e 

R-e 

e  <  r  <  R 

(11*1) 

II 

*T  V. 

^■0 

0 

0  <  r<  e 

(ll*2) 

, 

dr 

1 

R-e 

e  <  r<  R 

(11*3) 

9?  = 

II 

11 

0 

(11*1*) 

where  the  lagging  hinge  is  at  r=  e  .  In  addition,  I2Z  =  0  . 
Then, 


(1U5) 


where 


m  (r-e)2  dr 


fR 

So  =  /  m(r-e)dr 
Je 


(11*6) 

(lW) 

(iha) 

( 1U9) 
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Rigid  Blade  Pitching  Frequency 


For  a  pitch  mode  F( 


0 


(150) 


In  addition,  it  is  assumed  either  that  the  center  of  gravity  is  coincident 
with  the  elastic  axis  (i.e.,  x=Q  )  or  that  the  elastic  axis  is  coinci- 


dent  with  the  OY  axis,  (i.e., 

€  +  oEAb  =  0  ) .  Then , 

"T„  = 

Izz  -  Ixx 

(151) 

®M„  = 

Izz  +  Ixx 

(152) 

to  _ 

a  " 

/'  Vlj 

/ 

y  .♦$*) 

(153) 

where 

^xx  = 

i  I“dr 

(15*») 

Izz  = 

(155) 

for  a  plane  lamina,  Ixx  is  zero,  therefore  oj/Q,  -  !  . 


DISSIPATION  FUNCTION 


In  order  to  avoid  the  use  of  complex  coefficients  in  the  equations  of 
motion,  energy  dissipation  is  assumed  to  occur  due  to  viscous  structural 
damping,  which  is  the  usual  assumption  made  when  using  the  decay -rate 
method.  For  the  small  amount  of  structural  damping  present  in  metal  blades 
(5  =  0.03  approximately)  it  makes  little  difference  if  into  the  equations 
of  motion  one  puts  hysteretic  structural  damping  with  g  =  0.03  or  viscous 
structural  damping  with  C/Cc  =  0.015.  Hence,  the  generalized  structural 
damping  force  acting  in  the  rth  generalized  coordinated  can  be  represented 
in  the  form 


AQr  =  -Dgfrqr 


(156) 
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It  is  convenient  to  obtain  Dgrr  from  the  single-degree-of -freedom 

relation 


=  Qr< 


/GMrr  cu? 


‘Mrr=  gr  <*>r 


where  6Mrr  =  generalized  inertia  of  the  mode  used  as  the  r^h 
generalized  coordinate 

mr  =  natural  frequency  of  the  rth  mode 

gr  =  structural  damping  coefficient  of  the  r^h  mode 

ELASTIC  STIFFNESS 

The  elastic  energy  due  to  small  bending  and  torsional  deflections  of  a 
blade  can  be  expressed  as 


( 156) 


Note  that  the  units  of  EIX  ,  EIZ  ,  and  GJ  must  be  consistent  with 
those  of  the  deflections.  The  last  terra  on  the  right,  arises  from  the 
initial  tension  stresses  in  the  blade,  ka  is  the  polar  radius  of 
gyration  of  the  tensile  cross  section.  The  generalized  elastic  stiffness 
can  be  expressed  as  a  part  independent  of  the  rotational  speed  to  and 
a  part  dependent  on  to  .  These  are,  respectively, 


(160) 


Oh 


A  generalized  elastic  tension  stiffness  coefficient  L)2  is  determined 
from  the  potential  energy  Eq.  (162)  by  expressing  the  torsional  deflections 
in  the  modal  form  of  Eq.  ( Qq )  and  then  differentiating  as  indicated  in 
Eq.  Q60I. 


U 


J  I  (w)  {jfa2mrdr}dr 

{/mrdr}dr 


(161) 

(162) 


Hence 


ft'u„  ■  ft'u,  ■ 


(163) 
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o 

o 

o 

o 
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o 

0 

0 

0 

c 
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0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 
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0 
0 
0 
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*IM\ 


v't 

*M) 

Ft 

!£«&» 


dr 


(16U) 


TOTAL  STRUCTURAL  STIFFNESS 

It  will  be  assumed  that  the  generalized  modal  inertia  and  natural  fre¬ 
quency  of  each  of  the  modes  are  known  for  the  blade  rotating  in  a  vacuum 
at  a  particular  rotor  speed  .  Then  the  integration  implied  in 

Eq.  (159)  for  the  formulation  of  the  Kjj  elements  can  be  avoided. 

In  general,  the  mode  shapes  used  as  generalized  coordinates  will  be  con¬ 
sidered  to  be  semirigid,  a  term  used  in  fixed-wing  flutter  to  mean  that 
the  same  mode  shapes  are  used  for  all  forward  speeds.  Irrespective  of 
whether  the  generalized  coordinates  are  arbitrary  modes  or  orthogonal 
modes  at  zero  forward  speed,  they  will  be  arbitrary  modes  to  some  extent 
for  non-zero  forward  speeds.  However,  if  a  sufficient  number  of  generalized 
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coordinates  are  used  in  the  solution,  then  for  any  forward  speed  these 
arbitrary  modes  will  give  the  correct  coupled  mode  shapes.  The  question 
is:  Just  what  comprises  a  sufficient  number?  Fortunately,  the  effect  of 
forward  speed  on  the  mode  shapes  is  small,  so  that  quite  often  two  gener¬ 
alized  coordinates  are  capable  of  giving  a  good  approximation  to  the 
flutter  speed.  In  addition  to  the  effect  of  forward  velocity,  it  is 
necessary  also  to  consider  variations  in  the  mode  shapes  due  to  changes  in 
any  other  parameter  which  affects  the  mode  shapes.  Here  the  assumption 
of  semirigidity  is  not  always  so  satisfactory,  as  it  may  require  the  use 
of  an  excessive  number  of  generalized  coordinates.  In  such  cases,  it  is 
better  to  use  a  reduced  number  of  modes  and  to  assume  the  modes  to  be  semi¬ 
rigid  over  only  a  limited  range  of  variation  of  the  parameter  and  to  use 
new  mode  shapes  for  the  generalized  coordinates  as  necessary. 

The  concept  of  semirigidity  can  be  applied  also  to  rotating  blade  flutter 
where  rotor  speed  enters  as  one  of  che  parameters  affecting  mode  shapes. 

In  this  case,  it  is  desirable  to  be  able  to  assume  semirigidity  with  re¬ 
spect  to  both  forward  speed  and  rotor  speed.  This  is  facilitated  if  each 
mode  involves  displacements  in  1  degree  of  freedom  only;  that  is,  if  each 
mode  is  a  pure  mode.  It  has  been  found  in  practice  that  modes  involving 
only  pure  flatwise  bending,  or  pure  edgewise  bending, or  pure  torsion  have 
shapes  which  vary  little  over  a  wide  range  of  forward  velocity  and  rotor 
speed. 

Considering  the  mode  shapes  tc  be  semirigid,  the  total  structural  stiffness 
acting  in  a  generalized  coordinate  is  composed  of  an  elastic  stiffness 
which  is  constant  and  an  elastic  tension  stiffness  and  centrifugal  stiffness 
which  are  both  proportional  to  Az  .  Suppose  that  the  modes  are  calcu¬ 
lated  for  rotational  velocity  .  Let  a)R,  be  the  natural  frequency 

of  the  ith  mode  and  Mjj  the  corresponding  generalized  inertia.  Then  the 
generalized  total  structural  stiffness  at  the  angular  velocity  is 

given  by 

[Cs]„=  [  SMiK]  =  [k]+4[l]  '  ArM  (165) 

where 

is  a  diagonal  matrix  of  generalized  masses  multiplied 
by  the  corresponding  natural  frequencies  squared, 
appropriate  to  the  rotational  speed  ilR 

is  the  generalized  elastic  bending  and  torsion 
stiffness  matrix 

is  the  generalized  elastic  tension  stiffness  matrix 
is  the  generalized  centrifugal  stiffness  matrix 


32 


coordinates  are  used  in  the  solution,  then  for  any  forward  speed  these 
arbitrary  modes  will  give  the  correct  coupled  mode  shapes .  The  question 
is:  Just  what  comprises  a  sufficient  number?  Fortunately,  the  effect  of 
forward  speed  on  the  mode  shapes  is  small,  so  that  quite  often  two  gener¬ 
alized  coordinates  are  capable  of  giving  a  good  approximation  to  the 
flutter  speed.  In  addition  to  the  effect  of  forward  velocity,  it  is 
necessary  also  to  consider  variations  in  the  mode  shapes  due  to  changes  in 
any  other  parameter  which  affects  the  mode  shapes .  Here  the  assumption 
of  semirigidity  is  not  always  so  satisfactory,  as  it  may  require  the  use 
of  an  excessive  number  of  generalized  coordinates.  In  such  cases,  it  is 
better  to  use  a  reduced  number  of  modes  and  to  assume  the  modes  to  be  semi¬ 
rigid  over  only  a  limited  range  of  variation  of  the  parameter  and  to  use 
new  mode  shapes  for  the  generalized  coordinates  as  necessary. 


The  concept  of  semirigidity  can  be  applied  also  to  rotating  blade  flutter 
where  rotor  speed  enters  as  one  of  che  parameters  affecting  mode  shapes. 
In  this  case,  it  is  desirable  to  be  able  to  assume  semirigidity  with  re¬ 
spect  to  both  forward  speed  and  rotor  speed.  This  is  facilitated  if  each 
mode  involves  displacements  in  1  degree  of  freedom  only;  that  is,  if  each 
mode  is  a  pure  mode.  It  has  been  found  in  practice  that  modes  involving 
only  pure  flatwise  bending,  or  pure  edgewise  bending,or  pure  torsion  have 
shapes  which  vary  little  over  a  wide  range  of  forward  velocity  and  rotor 
speed. 


Considering  the  mode  shapes  to  be  semirigid,  the  total  structural  stiffness 
acting  in  a  generalized  coordinate  is  composed  of  an  elastic  stiffness 
which  is  constant  and  an  elastic  tension  stiffness  and  centrifugal  stiffness 
which  are  both  proportional  to  Slz  .  Suppose  that  the  modes  are  calcu¬ 
lated  for  rotational  velocity  ftp  .  Let  cuR,  be  the  natural  frequency 
of  the  ith  mode  and  Mji  the  corresponding  generalized  inertia.  Then  the 
generalized  total  structural  stiffness  at  the  angular  velocity  SIr  is 
given  by 


[Cs]R*  [GMii"Ri]  =  [K]+4[L]-fl2R[T]  (165) 


where 

[  gm„ <4,  ] 


is  a  diagonal  matrix  of  generalized  masses  multiplied 
by  the  corresponding  natural  frequencies  squared, 
appropriate  to  the  rotational  speed  &r 


fx) 

4M 


is  the  generalized  elastic  bending  and  torsion 
stiffness  matrix 

is  the  generalized  elastic  tension  stiffness  matrix 
is  the  generalized  centrifugal  stiffness  matrix 
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AERODYNAMIC  GENERALIZED  FORCES 


The  velocity  of  the  airstream  relative  to  the  axis  X,  Y,  Z  fixed  in  the 
helicopter  is  specified  by  the  two  components,  as  shown  in  Figure  9. 


vH/  is  the  velocity  in  the  direction  of  the  X  axis 

>/0  is  the  velocity  in  the  direction  of  the  (-Z)  axis 

The  velocity  of  the  airflow  relative  to  the  blade  is  given  in  terms  of  the 
triad  of  unit  vectors  "ei,  e^,  fixed  in  the  blade,  where 

e^  lies  along  the  blade  radial  axis 

is  in  the  plane  of  the  shaft  axis  OZ  and  the 
rotating  axis  OX 

e^  completes  the  right-handed  triad 

Then,  with  reference  to  Figure  9,  the  velocity  of  the  airstream  relative  to 
the  blade  is  •tiven  by 


v  =  (vH  co$/3j  cos yf,  -  v0  sir,#s)ef 


-  [vHsin^  +  &(e+  rcos/3s)]e^ 

-  [  vH  sm/3s  cosi/r  +  vD  cos/3s] 


(170) 


The  air  velocity  component  in  the  e^  direction  represents  radial  flow  along 
the  blades;  for  the  aerodynamic  model  being  considered,  it  does  not  produce 
any  aerodynamic  force.  Hence,  it  is  neglected.  The  magnitude  of  the  air 
velocity  normal  to  the  blade  is  therefore  given  by 


V  =*/[vH  sin^  +  ft(e  +  r  cos/3s)]2  +  [vH  sm^cos  f  +  v0cos/3s]2  (171) 


The  coning  angle  0S  is  included  in  the  derivation  of  the  aerodynamic 
coefficients.  After  these  coefficients  had  been  programmed,  it  was  realized 
that  the  inclusion  of  the  coning  angle  /3S  would  complicate  unnecessarily 
the  derivation.:  of  the  structural  coefficients,  and  so  it  was  deleted. 

Hence,  there  is  the  inconsistency  in  that  the  coning  angle  is  included  in 
the  aerodynamic  coefficients  but  not  in  the  structural.  This  inconsistency 
is  removed  \y  always  specifying  /3S  =  0  . 

With  Reference  to  Figure  10,  the  virtual  work  done  by  the  lift  L  and  the 
moment  Me/*  acting  on  an  airfoil  section  of  width  dr  due  to  a 
virtual  displacement  is  given  by 


(173) 


SW  =  b[-L8(4)  +  8a] 

SW  =  Qh  3  (-£-)  +  Qa8o 


where  and  Qq  are  the  generalized  aerodynamic  forces  due  to 

vertical  displacement  and  pitch  respectively. 


Hence, 


'Qh' 

=  b 

"-L 

Mc/4 

Qa 

b 

mS 

(175) 


The  aerodynamic  coefficients  we^e  previously  defined  by  the  equation 


-L 
Mc/4 
L  b  J 


*  -2^drfe)[D]{8}t(^[D]{8}} 


where 


Also,  by  definition, 


W  -F4 

ta 

"ma. 

a 

i — 
1 

3 

vr 

“ma  J 

W  ■  Rr 


h_ 
b 

a  J  C/4 


where 


hEA  bR  I  fj  (r)  qj  (t) 

=  ZF.l,)q(t) 

'  J  ‘  ) 


a  = 


(176) 


(177) 

(178) 

(179) 

v 

(160) 

(181) 


bR  fj  (r)  =  -g?(r) 


(182) 
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Similarly , 


[gm]  =  r fi  ■  ■ '  *  f*° 

/,••••  ^H0_ 

(191/ 

{q}  =  r  q.  ' 
qND 

w  * 

(192) 

M  -  [(^MNWl'tos) 

(3.93) 

■  (¥)N'N[h]'[qJ 

.Qq. 

(19*+) 

Combining  equations  (I75),  (176),  (l88),  and*(l9lj)  gives 
W  .  -a.,(^^r[G„]'[Bc][H]'fe)[D][H][Bc][G„]{q} 

♦ODWMbJWW 


(195) 

When  the  following  substitutions  are  introduced 
cr  =  A 

P sl  U96) 

dr  =  Idrj  (197) 

W--feW'[Bc][H]'[D][H][B0][GM] 

W*  ^M[8c][h]'[d][h][bc][gm]  (i99) 

[bU  ■  £$)[*]  * 

\v«/  L  J  (200) 

[C]aero  =  jT  °(^-)  [*a]  d?7  (201) 

and  when  integrated  to  find  the  generalized  force  due  to  the  whole  blade, 
the  equations  become 


w-  -2^v^{(r:)[B]AJic]jq}} 


(202) 
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•PROVISION  .FOP,  MASS  BALANCE 

Mass  balance ,  which  is  distributed  radially  along  the  blade,  is  considered 
separately  iron  the  mass  of  the  basic  blade  structure.  When  this  is  done, 
the-  following  expressions  for  the  structural  coefficients  are  obtained  in 
matrix  notation: 


(203) 

[0]  - 

( 20U) 

M  ■  fM'fr'-WH}  M* 

(205) 

[t] 

(206) 

where:  subscript  refers  to  the  basic  blade  structure,  subscript  b  refers 
td  the  mass  balance  weight,  and  tn  is  a  nondimens ional  parameter  to  vary 
the  amount  of  blade  mass  balance.  For  tn  =  0  mass  balance  effects  become 

zero.  The  form  of  - - -  [T[j]  can.  be  obtained  by  inspection  of 

the  expressions  for  [6M  J  ,  [g]  ,  [l]  ,  and  [t]  obtained  previously. 

The  matrix  ^[GMdb]  and  the  other  matrices  on  the  right  hand  sides  of 
Eq.  (203)  through  (206)  are  determined  by  using  expressions  which  ere  com¬ 
pletely  analogous  to  those  found  previously  without  consideration  of  sep¬ 
arate  mass  balance. 

PROVISION  FOR  VARIATIONS  IN  BLADE  PARAMETERS 


The  provision  to  vary  the  following  parameters  is  now  incorporated  into 
^fche  expressions  for  the  coefficients: 

1.  Frequencies  of  the  modes  used  as  generalized  coordinates. 

2.  Blade  chordwise  center-of-gravity  positron. 

3.  Blade  chordwise  elastic  axis  position. 

1».  Blade  chordwise  aerodynamic  center  position. 

Variation  of  Blade  Modal  Frequencies 

The  generalized  structural  stiffness  at  the  rotor  speed  JIr  is  given  by 
[Cs3r  =  [GMjj  cur|J  .  The  modal  frequency  is  considered  in¬ 
creased  by  the  amount  yr  ofvi  where  o;Vj  is  specified  for  each 
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mode  and  yr  is  an  overall  variation  factor.  The  new  structural  stiff¬ 
ness  is  then  given  by 

[cs]R  =  [  6Mii(a>R,  +  yvWvi)Z  ]  (207) 


Variation  of  Blade  Chordwise  Center  of  Gravity  and  Elastic  Axis  Position 

The  structural  model  to  be  used  when  changing  the  center  of  gravity  and 
elastic  axis  positions  is  as  follows: 

The  chordwise  center-of-gravity  position  changes  due  to  a  mass  balance 
weight  which  is  moved  chordwise.  The  total  mass  of  the  blade  remains  un¬ 
changed,  as  do  the  mass  and  inertial  properties  of  the  fixed  portion  of 
the  blade.  The  fixed  or  basic  part  of  the  blade  and  the  balance  weight 
are  entered  into  the  structural  coefficients  as  separate  components,  with 
the  magnitude  of  the  balance  mass  made  variable  through  use  of  the  coef¬ 
ficient  tn 

The  input  mode  shapes ,  frequencies ,  and  generalized  inertias  remain  un¬ 
changed.  The  effect  of  an  elastic  axis  shift  is  to: 

V 

1.  Change  the  mass  moment  and  inertia  about  the  elastic  axis,  but 
not  the  mass  moment  and  inertia  about  the  center  of  gravity. 

This  effect  will  be  taken  into  account  when  blade  inertias  Eire 
calculated. 

2.  Change  the  aerodynamic-center-to-elastic-axis  distance.  This  , 
is  taken  into  account  by  changing  the  value  of  aEA  when 
the  aerodynamic  coefficients  are  calculated. 

With  reference  to  Figure  11, 

tb*x  is  the  distance  the  center  of  gravity  of  the  mass 
balance  weight  is  shifted  aft. 

0 

tebAa£A  is  the  distance  the  elastic  axis  is  shifted  aft. 

Then,  for  the  balance  weight, 

*b  =  *b  +  fb^  "  ■«  b^°EA  (208) 


and  for  the  basic  blade, 

=  *<f>  ~  U  bAaEA 
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(209) 


! 


The  elastic  axis  distance  aft  of  the  mid-chord  is  given  by 

flEA=  aEA+  AaEA 


(210) 


First,  the  changes  in  the  blade  inertias  referring  to  tha 
the  elastic  axis  are  calculated. 


new  position  of 


The  mass  moment  per  unit  span  of  the  balance  weight  prior  to  the 
tion  of  the  coefficient  tn  is 


applica¬ 


nt)  xb 


mb  xb 


mb(  t^x  -  te  bAaEA) 


(211) 


The  seme  quantity  for  the  basic  blade  is 

bAa^ 


(212) 


The  balance  weight  edgewise  mass  moment  of  inertia  per  unit  span  is 

Izz>>  *  x«bCG  +  '"bl*!,')2  +  2<tbas  -t.baoE4)mbxb' 

*  ( tjj Ax  -  t*bAaEA)  (213) 

*zzb  =  Izzb  +  2(tbAx  -  tebAa^Jm^h  +  mb(tbAx  -  tebAaEA)2  (21 1() 

The- corresponding  quantity  for  the  basic  blade  is 

Izz^(  =  Izz<£  -  StebAa^m^  +  m^(tebAaEA)2  (215) 

The  balance  weight  product  of  inertia  per  unit  span  is 

rxzb  =  Ixzbc6  +  rnb[xb  +(tfaAx  -  tebAaEA)]zb  (2l6) 

Ixzb  =  Ixzb  +  mb(tbAx  -  tebAaEA)zb 


(r.-7) 


The  corresponding  quantity  for  the  basic  blade  is 


Ixz  <j> 


Ixz^  -  tebAaEAm^,z^ 


(218) 


The  inertial  properties  mb  Zb  and  IXxb  for  the  balance  weight  and 
mz  and  Ixx  for  the  basic  blade  are  not  changed.  The  changes  in 
the  aerodynamic  coefficients  due  to  shifti.;g  the  elastic  axis  are  now 
calculaued.  This  is  done  by  changing  the  value  of  oEA  used  in  the 
calculation  of  the  [h]  matrix,  which  is  thereby  redefined  as 


[h]  ■ 

1  -  (£+°eA+  teAoEA )' 

0  1 
*  - 

(219) 


Variation  of  Blade  Chordwise  Aerodynamic  Center  Position 

In  addition  to  arbitrary  variation  of  the  position  of  the  aerodynamic 
center  there  is  also  a  shift  in  the  aerodynamic  center  due  to  reverse  flow. 
These  effects  will  be  considered  separately.  Consider  first  an  arbitrary 
variation  in  the  position  of  the  aerodynamic  center .  Suppose  the  aero¬ 
dynamic  center  shifted  aft  a  distance  tobSAC  ,  as  shown  in  Figure  12. 

In  classical  flutter,  changes  in  lift  and  moment  due  tv,  very  small  motions 
of  the  blade  are  considered.  For  steady  flow,  the  aerodynamic  center  is 
not  coincident  with  the  center  of  pressure  except  for  airfoils  with  a  zero 
steady  moment  at  zero  lift.  However,  when  considering  incremental  changes, 
the  aerodynamic  center  and  the  center  of  pressure  become  coincident.  In 
flutter,  the  two  terms  are  used  interchangeably  and  usually  refer  to  the 
point  through  which  the  incremental  lift  force  acts,  for  nonoscillatory 
motion.  This  definition  of  the  aerodynamic  center  means  that  for  oscilla¬ 
tory  motions  there  will  in  general  be  an  oscillatory  moment  acting  at  the 
aerodynamic  center  (center  of  pressure).  A  shift  in  the  aerodynamic  cen¬ 
ter  is  defined  to  be  such  that  the  forces  and  moments  acting  at  the  old 
position  of  the  aerodynamic  center  are  transferred  unchanged  to  the  new 
position  of  the  aerodynamic  center. 


That  is, 


il 

l  b 

where  the  quantities  on  the  right  hand  side  are  applied  at  a  distance  Ab 
aft  of  the  25#  chord  and  the  quantifies  at  the  left  hand  side  are  applied 
at  a  distance  (Ab  +  tabSAc)  aft  of  the  25#  chord,  as  shown  on  Figure  12. 
The  position  of  the  aerodynamic  center  will  be  denoted  henceforth  by  using 
the  superscript  notation  above. 


1  u 


b 


(220) 


1*1 


‘  v 


Substituting  Eqs.  f 2i£4)  and  (225)  into  (223)  and  setting  the  respective 


coefficients  of 


w 


and 


equal,  gives 


w 

[*]"■  [A]  [6] 


(226) 


where 


[of  -  MM 

[A]  .  f  ■  0 

j'oS  ' 


(227) 


(228) 


Consider  now  the  effect  of  reverse  flow.  Only  those  quantities  that  affect 
the  aerodynamic  coefficients  are  involved,  since  the  structural  coefficients 
are  the  same  for  forward  and  reverse  flow.  Inspection  of  Figure  13  shows 
the  following  relations  between  forward  and  reverse  flow: 


•eAasA f-  ~ 

(°L+  teAaEA)F 

(229) 

MR  • 

(h  ea)F 

(230) 

(«)*  * 

-(a)F 

(231) 

where  superscript  F  refers  to  forward  flow  and  superscript  R  refers  to  re¬ 
verse  flow  . 

The  change  in  sign  of  the  pitch  angle  a  shows  up  as  a  change  in  sign 
of  the  torsion  mode  shape  Fj  through  the  equation 


Q  =  ^  Fj  qj  (232) 

Reverse  flow  is  3aid  to  occur  when  the  component  of  the  air  velocity  in  the 
(  -e£  )  direction,  as  previously  defined,  is  negative.  That  is, 

vx  =  sirnJ/vH  +  &  (e  +  rcos/3s)  <°  (233) 


l<3 


Hence,  when  Vx  <  0  ,  reverse  flow  is  allowed  for  in  the  following  way: 

°  1.  Change  the  sign  of  (a£A  +  ^Aa^)  in  the  matrix 

[H]  so  that  it  becomes 


[h]  =  1  -  [r "  (°L+  feAaEA)] 

0  I 


(23*0 


2.  Change  the  sign  of  Fi  in  the  second  row  of  the  matrix 
so  that  it  becomes 


[gm]  8  f • 

^ND 

(~Fi) - 

”  "  (~Fnc) 

(235) 


Similar  changes  are  also  made  in  [h]7  and  [g^]  1  . 

SPAIWISE  INTEGRATION  OF  THE  COEFFICIENTS 

So  far,  the  generalized  forces  acting  in  a  mode  have  been  expressed  in 
terms  of  integrations  across  the  blade  span.  The  same  procedure  has  been 
followed  in  the  definition  of  the  generalized  mass  and  stiffness  coef¬ 
ficients.  In  practice,  the  integrations  are  replaced  by  finite  summations, 
using  the  trapezoidal  rule.  The  blade  is  divided  into  segments  or  strips, 
each  of  which  extends  from  the  blade  leading  edge  to  the  trailing  edge, 
and  whose  radial  or  spanvise  dimension  is  a  suitable  fraction  of  the  blade 
rddius.  The  segments  are  labeled  k=l  to  k =  Ns  from  root  to  tip,  re¬ 
spectively.  Considering  the  structural  coefficients  and,  for  convenience, 
using  the  form  of  the  coefficients  written  in  matrix  notation  gives 


/  [*]  A[«]|>]*  ■  £  M/HH,  Ar« 


(236) 


The  matrix  [a]  represents  any  of  the  structural  coefficient  matrices. 
In  addition, 


f " 

M  (r )  =  / 


(237) 


is  replaced  by 


=  4-  m  (rK+  4"ArK)ArK+ ^  m  rn^rn 


(238) 


The  values  of  m  in  Eq..  (238)  are  expressed  in  mass  per  unit  span  and 
are  appropriate  to  each  finite  segment.  When  the  integrals  are  replaced 
by  finite  summations ,  it  is  convenient  to  change  the  definitions  of  the 
mass  terms  m  ,  mx  etc. ,  from  "per  unit  span"  to  "per  segment"..  The 
definition  of  Mk  is  not  changed,  since  it  refers  not  only  to  a  local 
part  of  the  blade,  but  also  to  a  summation  involving  all  parts  of  the 
blade  outboard  of  a  particular  segment.  Because  of  this,  the  value  of  Mk 
cannot  be  given  easily  as  input  data  on  the  basis  of  a  contribution  from 
each  segment,  and  its  value  at  each  segment  is  computed  by  the  program. 

Let  m  ,  fnx  etc.,  be  the  inertias  per  unit  span  when  integrating.  Let 
mk  ,  mkx  etc.,  be  the  inertias  per  blade  segment  when  a  finite  summa¬ 
tion  is  used.  Let  [a]  be  per  unit  span.  Let  [5]  be  per  blade  segment. 
Then,  the  above  equations  become 

(239) 


Mk  =  i  mK(rK+5-ArK) 


*2 


n=K+l 


(2l»0) 


and  any  element  of  [a]k  containing  Mk  is  multiplied  by  Ar^ 
When  the  aerodynamic  coefficients  are  considered. 


(2Ul) 


(2142) 


where 


[4  ■  (!)!>]„'  [Bc]k  [h]'  [a]„  [6]k  [H]k  [Bc]k  [g„]( 


(2l»3) 


[4  =  [gm]'  [Bc]k  [h  ]'  [A]k  [d]k  [H]k  [Bc]k  [g„]. 


(2UU) 


NONDIMENSIONALIZING  THE  EQUATIONS  OF  MOTION 

The  equations  of  motion  as  developed  so  far  can  he  put  in  the  fora 

w  m  *  wm 

-  -  2l%.bR  f  1  {(' *)  [8  IJl)-}  U 1  }} 

where  A»  BS  and  C$  are  the  total  structural  coef¬ 
ficients.  L  J 

I  /bR\2 

The  equations  are  now  multiplied  by  27r/jSLbJ  JL  \vn/  and  the 
nondimensional  time  r  is  introduced  by  the  relation 


(2U5) 


so  that 


A  -  (Y*.)  A 

dt  '  bR  •  dr 


■  - 


=  (Yr.)2  _i_2 

dtfc  VbR/  dr2 


(2U6) 


(2U7) 


After  some  rearranging,  the  equations  become 

Mav]  {^}+  {(vj)  *  t»  !«.}{&} 


(2U8) 


FINAL  FORM  OF  THE  EQUATIONS  OF  MOTION 


The  final  form  of  the  equations  of  motion  can  now  be  stated,  as  follows: 


[*],{£}♦  [B]r{£Mc]T{<}=  {0} 


(2U9) 


k6 


where 


|q  j-  is  u  column  matrix  of  generalized  coordinates  q, . qN0 

r  =  (^)  t  (250) 

Wf  [»]-!«[']*['L  'W1) 

[c]t  =°[Cs]r  *  -  125)  {°[l]  -°[t]}  +  [CLR0  (252) 


=  (  gf, g5j  +  gJ,gJj)  m<#>*  -  ( gkI Fkj  +  g^ Fki )  m#x# 

+  (gU)Fki+  gJjFw)^kZ^k  +(^r  ^)(IxxV 

+  (^ds|,M^)(i„v 

+  (¥  +  (Fki  Fk))(lyy)* 


(261) 


(262) 


9k*  • 

- 

-  - 

-  gx*D 

£ 

*Iho 

afxte) 

ffkkN,,) 

d^W*) 

dF^KNo) 

Fki 

Fkk 

FkND 

dF  M 

i  (*j 

• 

i,  (F*J 

rn^,k 

0 

0 

0 

(-mx)^k 

0 

0 

m$k 

0 

0 

(mz)^x 

0 

0 

0 

(Iwl* 

(Jxz)^ 

0 

0 

0 

0 

(*Xz)$k 

(^zr)^k 

0 

0 

-(mx)$k 

(mf)^k 

0 

0 

M^k 

0 

0 

0 

0 

0 

0 

0 

(263) 


(261*) 


1*8 


(%k)ij  8  (o{|95j)m+k  +  (gKV:ki  +gkX]Fki)m^^k 

+  (%*■  Fk)  +  ^  FKi  )  r# 

+  V“HT  “5^)  ^xxJfk  -^<#.kArk] 

+  TT^i  [(hz)<t>*  -  Jb&br k  ] 

A  /'dg,*,  dgzk,  dg£  dgzki\ 

+  \"ar"3r  + 

+  ^ki  ^kj  [(^xx)^k-“  (lzz)^,k_  (c  +aEAb)m<^x^k] 


(271) 


h,  ■  w:  h>  w, 


(272) 


000  0  0  0 

0  0  0  (m*4k  0 

0  0  ^(l»»)^k~i/^t^kAr|lj  (Ixz)^k  (rmk)^  0 

0  0  (I  »r)^v  [(I:z)^-^*ArJ  -(rmz)^  0 

0  (mZ)$*  (rmX)^k  -(rmz)^k  [(I»»)^k"(I«Vk"(«+a£»bXmX)^kJ  0 
0  0  •  0  0  0  0 


(273) 
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The  expressions  for  the  mass  balance  coefficients  are  identical  to  the 
expressions  given  above  except  that  subscript  <p  is  replaced  by  b. 


[4  “  <*ifc  fe)  b];  [bc]r[h]'  [a],  [&]k  [h]„  [bc]„  [gm]„ 
[B4  K  [a],  [d](  [h]„  [Bc]k  [gm]„ 


-[^  i(aEAK+te^aEAK)] 
I 


(27*4) 

(275) 

i 

(276) 

(277) 

'(278) 

(279) 

(280) 

(281) 
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Reverse  flow: 


=  sin  \J/  vH  +  XI  (e  +  rKcos  /3S)  <  o  (282) 

Then  use  the  negative  signs  in  the  expressions  for  [h]^  and  [Gm] ^  . 

Examples  of  the  application  of  the  reverse  flow  relations  are  given  in 
Table  1. 

The  cases  listed,  except  those  at  MN=  0.3,  show  satisfactory  agreement 
betwee  -  the  coefficients  from  diff en nt  references ,  and  the  reverse-flow 
relations  are  satisfied.  At  MN=  0.8,  Reference  7  satisfied  the  reverse 
flow  relation  for  calculating  /h  but  not  the  relation  for  calculating 
.  For  this  case,  the  greatest  discrepancy  was  in  the  coefficient 
(-m<i  )  ,  a  few  values  of  which  are  given  in  Table  II. 

A  large  proportion  of  the  tabulated  coefficients  come  from  Reference  7 , 
and  it  is  unfortunate  that  these  coefficients  in  many  cases  satisfy  only 
partially  the  reverse-flow  relations. 

When  first  published  (Reference  5)>  they  were  criticized  for  not  satisfying 
these  relations.  Subsequently,  an  amended  tabulation  appeared  (Reference 
6),  and  it  was  claimed  that  this  deficiency  had  been  corrected. 

About  the  same  time,  Reference  7  appeared;  it  gave  a  much  more  complete 
table  of  coefficients,  but  without  claiming  that  they  satisfied  the 
reverse-flow  relations.  It  was  assumed  that  they  did;  indeed,  comparison 
between  the  coefficients  tabulated  in  Reference  7  and  the  revised  coef¬ 
ficients  in  Reference  t  show  them  to  agree  exactly  or  with  differences 
confined  to  the  last  decimal  place.  The  matter  was  not  pursued  any 
further ,  though  it  would  certainly  be  of  interest  to  do  so. 

With  the  reverse-flow  relations,  it  was  possible  to  sort  out  the  coef¬ 
ficients  satisfactorily,  and  a  table  of  coefficients  was  prepared.  These 
coefficients  are  given  in  Table  XII.  It  will  be  noticed  that  the  coef¬ 
ficients  for  M  =  1.05  exhibit  small  oscillations  with  respect  to  k. 

N 
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AERODYNAMIC  COEFFICIENTS 


GENERAL  DISCUSSION  OF  AERODYNAMIC  COEFFICIENTS 


The  aerodynamic  forces  acting  on  a  rotor  blade  appear 
motion  in  the  form  of  the  matrices 


Mr 


ices  [6] 


and 


k 


la 


M 


in  the  equations  of 
where 


k 


(283) 


(284) 


That  is,  at  each  blade  station,  eight  aerodynamic  coefficients  are  required 
to  specify  the  aerodynamic  forces  acting  on  the  blade  segment.  Each  of 
these  coefficients  is  in  general  a  function  of  the  reduced  frequency  param¬ 
eter  k=  cnb/V y  and  the  Mach  number.  Depending  on  the  forward  velocity 
and  the  advance  ratio,  the  value  of  k  and  Mn  at  any  station  will 

be,  in  practice,  within  the  range  0  <  k  <  CO  ‘  and  0<  MN<  I. 

For  certain  combinations  of  VH  and  ^  ,  k  and  Mm  will  have  large 

variations  across  the  blade,  so  that  it  is  necessary  to  form  a  table -of 
aerodynamic  coefficients  covering  a  wide  range  of  k  and  MN  .  The  values 
of  k  and  MN  for  each  segment  are  calculated,  and  by  interpolation  the 
corresponding  value  of  each  of  the  eight  aerodynamic  coefficients  is 
obtained. 


A  search  of  the  literature  located  tables  of  incompressible  and  compressible 
main  surface  aerodynamic  coefficients  in  References  2,  3,  4,  5,  6,  7,  8, 
and  9.  The  subsonic  compressible  coefficients  were  obtained  either  by 
direct  iteration  of  Possio's  integral  equation  or  by  expanding  the  solution 
as  an  infinite  series  of  orthogonal  Mathieu  functions.  Either  method  of 
solution  has  the  characteristic  that  the  convergence  becomes  poorer  as  k 
and  MN  increase.  Because  of  this,  the  range  of  k  for  which  tabulated 
derivatives  are  available  decreases  as  Mn  increases.  This  is  normally 
not  a  restriction,  since  the  range  of  k  values  across  the  blade  also  de¬ 
creases  as  the  blade  tip  Mach  number  increases.  However,  the  particular 
method  used  to  interpolate  for  the  coefficients  at  a  blade  station  requires 
that  the  coefficients  be  tabulated  for  the  same  values  of  k  for  all  Mach 
numbers.  Because  of  this,  full  use  was  made  of  the  available  tabulated 
coefficients  for  the  higher  Mach  numbers  but  not  for  the  low  Mach  numbers. 
This  restriction  is  most  severe  for  incompressible  flow,  where  the  aero¬ 
dynamic  coefficients  are  available  over  the  interval  0<  k  <  GO  ,  This 
semi-infinite  interval  can  be  tabulated  by  transforming  it  into  finite  in¬ 
terval  0  <  k  <  I  by  a  transformation  such  as  k  =  k/k+p  and  by  tab¬ 
ulating  the  coefficients  for  all  Mach  numbers  versus  k  instead  of  k. 
However,  for  the  present  study,  each  of  the  eight  coefficients  is  tabulated 
for  Mn  =  0,  .35,  .50,  .60,  .70,  .80,  .90,  .95,  1.0,  1.05  at  each  of  27 
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different  values  of  k  "between  k  =  0  and  k  =  0.7. 

For  reasons  which  will  he  discussed  below,  the  last  row  of  the  tabulation 
was  repeated  twice,  giving  29  k  values.  Thus,  the  tabulation  consists  of 
2320  values  of  the  coefficients.  It  was  stated  previously  that  the  vari¬ 
ation  of  the  aerodynamic  acceleration  coefficients  with  change  in  k  for 
compressible  flow  was  one  of  the  reasons  for  including  the  acceleration 
coefficients  with  the  displacement  coefficients .  If  they  had  been  kept 
separate,  it  would  have  been  necessary  to  store  3^80  instead  of  2320  values 
of  the  coefficients  in  the  high-speed  store  of  the  computer  when,  even  with 
2320  values,  space  was  at  a  premium. 

The  value  of  some  of  the  coefficients  approaches  infinity  as  k  0 
or  00  .  These  limiting  cases  will  be  considered  for  incompressible 

flow,  but  the  same  principles  apply  for  all  Mach  numbers. 

The  case  k  =  0  is  considered  first.  The  lift  due  to  pitch  rate  in  this 
case  is  given  by 

L  =  27T^b3Ar(-£)  4.d  (285) 

where 

4  =  O'  *  (286) 

"Tor  simple  harmonic  motion,  this  becomes 

L  =  '27r/>$L  b3Ar  |(t)  [i+F(k)]d  +  ($2  [i  (sjpG(k)j a  j-  (287) 

where 

oj  is  the  solution  frequency 

a»0  is  the  solution  frequency  assumed  when  calculating 

The  case  of  k  =  gjq  b  /  V.  -  0  ,  that  is,  steady- flow  aerodynamic  is  now 

considered.  In  practice,  steady  and  oscillatory  flows  occur  with  oj/ojq  =  I  ; 
and  considering,  the  fact  that  0.5  <  F(k)  <  1  and  -0.188  <  G(k)  <  0 
the  above,  equations  show  that  the  theory  predicts  finite  airloads  for  all 
values  of  k.  It  is  the  artifice  of  using  gj/cj0  5*1  in  particular 
tta  =  O  ,  for  the  aerodynamic  associated  with  vibration  modes  oj  £  0 
that  leads  to  apparently  infinite  airloads.  If  oj/ojq  -  I  and  since 
F  ( k )  =  I  ,  G  ( k'=  0.,  and  a  =  0  for  w  =  te o  =  0  the  above  equation 
gives  L  =  0  at  .  k  =  0,  so  that  the  contribution  to  the  steady-state 
loading  due  to  la  is  jero.  This  same  result  can  i>v  ob- 
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tained  with  (jj/oJa  &  I  if  L  J k  =  0  =  -  °0  is  replaced  by  [/aJk=0  s  ®  * 
That  is,  it  is  quite  valid  to  tabulate  /< j  =  0  at  k  =  0,  and  this  was 
done  in  the  table  of  aerodynamic  coefficients. 

The  difficulty  that  occurs  in  the  case  k  -  00  has  been  discussed  before 
and  is  due  to  lumping  the  acceleration  coefficients  with  the  displacement 
coefficients.  If  the  components  due  to  inertia  were  neglected  then  the 
displacement  coefficients  such  as  -ma  would  remain  finite  as,  k  —  CO 
because 


-ma  = 


(•288) 


is  due  solely  to  the  inertia  coefficient  md  =3/16  .  The  effect  of 

aerodynamic  inertia  is  generally  considered  to  be  of  secondary  importance, 
and  the  contribution  of  the  aerodynamic  inertias  will  be  progressively  re¬ 
duced  to  zero  by  making  all  coefficients  constant  above  Mfg  =  0.-7.  This 
avoids  the  uncertainty  of  extrapolating  curves  from  k  =  0.7  to  =  3,  say, 
when  the  curve  at  k  =  0.7  is  almost  ■vertical.  To  keep  the  coefficients 
constant,  the  tabulated  values  of  the  coefficients  at  k  =  0.7  are  repeated 
twice  and  labeled  k  =  0.8  and  k  =  0.9  respectively.  This  causes  the  linear 
extrapolation  routine  to  calculate  values  for  the  coefficients  which  are 
the  same  as  those  at  k  =  0.7*  The  effect  of  progressively  neglecting  the 
aerodynamic  inertias  above  k  =  0.7  is  further  reduced  by  noting  -the  follow¬ 
ing  from  the  relation  V:  ub/k  :  as  k  increases ,  V  decreases , 

so  that  the  region  k  >  0.7  is  fortunately  associated  with  very  low 
dynamic  pressures,  which  cause  the  coefficients  for  k  >0.7  to  make 
only  a  very  small  contribution  to  the  equations  of  motion. 

Two  dimensional  theory  gives  infinite  airloads  for  steady  flow  at  ~  1. 

This  singularity  disappears  for  oscillatory  flow  (  k  ^  0  )  or  with  the 
inclusion  of  finite  span  effects.  The  choice  has  to  be  made  between  tab¬ 
ulating  coefficients  at  Mjj  =  I  ,  k  =  0  for  a  particular  aspect  ratio  or 
simply  making  all  coefficients  zero.  For  simplicity,  the  latter  was  decided 
upon  so  that  the  tabulated  coefficients  cannot  be  used  for  steady  flow  at 
Mn  =  I 

NOTATION 


It  is  necessary  to  transform  the  tabulated  coefficients  available  in  the 
previously  stated  references  so  that  they  conform  to  the  notation  for  this 
study.  These  transformations  are  now  reviewed. 

Each  reference,  in  general,  uses  different  symbols  when  defining  the  coef¬ 
ficients,  so  for  convenience  the  defining  equations  will  be  stated  in  terms 
of  the  notations  already  adopted  for  lift,  pitching  moment,  semichord,  etc. 
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The  coefficients  from  Reference  2  will  be  converted  first  and  are  defined 

by  ■ 


( /* +  i  l ‘I )  {li+ila) 

-(m'+im")  -(n£  +  im") 


(28 9) 


with  the  relations 

2/3bV2Ar  =  (27T/3b3Ar)[^(^-)  ^  jw 

2/>bV2Ar  =  (27T/)b3Ar)  [(-£f^] 
and  for  simple  harmonic  motion. 


iojct  =  a 


also 


Then  the  equations  become 


(290) 

(291) 

(292) 
(293) 

(29*0 

(295) 


(296) 
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If  this  is  compared  vith  the  defining  equation. 


the  following  relations  between  the  coefficients  are  obtained: 
The  damping  coefficients  are: 


The  stiffness  coefficients  are 


(297) 

(298) 

(299) 

(300) 

(301) 

(302) 

(303) 

'30)*) 

(305) 
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(306) 


-  rrth  =  rr  [(-"4)- 

_ma  =¥{2(-m2)  +  ^  ]} 


/"  /"  »  ••  /  / 

The  coefficients  l z  ,  la  ,  -mz  ,  -ma  ,  oi  ,  t-a 

and  -ma  are  tabulated  in  Tables  2  through  2d,  pages  3h  through 

Reference  2. 


The  coefficients  in-  Reference  3  are  defined  by 


-1 

’  -  L~ 

(L  i+  i  Lz) 

(  L3+  i  L4) 

4/>bV2k2Ar 

id 

b 

EA 

(M,  +  i  M2) 

(  M3  +  i  M4) 

EA 


and 


_ -1 

~-L 

= 

’( 

lJ-hl,)  | 

1 - 

+ 

Jl 

L- 

4  /jb  V2  k2Ar 

M 

T 

LE 

K  +  iM,') 

(M'  +  iM'4) 

b 

a 

but-  only  the  primed  derivatives  are  tabulated. 

The  following  relations  are  used: 

4/>bV2  k2Ar  =  (27r/>b3Ar)  [(^f*  ]w 

4/>bV2  k2Ar  =  [zirp  b3Ar.)  [(£)*£'k2] 
For  simple  harmonic  motion, 


£ 

b 


icua 


a 


(307) 

l 

-mz  , 
of 

(308) 

(309) 

(310) 

(311) 

(312) 
(313) 


58 


M 

LE 


'  -m 


l  a 


Then  the  following  relations  "between  the  coefficients  are  obtained: 


The  damping  coefficients  are 


k  &  1  =  ^  [  i  0 Ilk2 La  k2L4l f I 


mh  "ma 


4  I J  Lk2  Mi  k2  Mi  0  I 


^  s  ih  (  ^  ^*2 ) 


4  =  If  (kZ^4— g-k2  Lg) 


H  ”  7Tk 


(k2Mj4k 2  4) 


"ma=  wTT  (k2f^  -- 2"[k2M2+k2L^]+  ^k2  Lg) 


The  stiffness  coefficients  are 


L  L~\  =  4  f  I  ol  fk2L'  k2  La  1  fl 


I  k2  m[  k2  Mg  0  I 


(31*») 


(315) 


(316) 


(317) 


(318) 


(319) 


(320) 


(321) 


V 


/h  =-f(k2L',) 


la  -'r  (k2L3  -  i  k2  L,i ) 

™h  *  #(kSM,'-|k2L') 

"ma  (h2M3-£[k2Mi  +  k2|4]  +  Tk2Li) 


(322) 


(323) 


(321.) 


(325) 


The  coefficients  k2  ,  k2  L4  ,  k2  Mg  ,  k2  M4  ,  k2  l(  ,  k2  l_3, 
kzM|  and  k2  M3  are  tabulated  in  Table  I,  page  9  of  Reference  3. 

The  coefficients  in  References  4  and  9  are  defined  by 


__  ~  -L 
kr 


-  (Lhn+'Lhj]  (LaR+  iLai)  ‘F 
C/4  -(MhR+lMhi)  (MctR+  'Mai)  a  C/4- 


The  following  relations  are  used: 


7r/3b3u»?Ar  =  "  {zirp  b3Ar)  GB*"!’]  w 
vp  b3w2Ar  =  -  (277/3b3Ar) 


For  simple  harmonic  motion, 


la,(^)=  i 


icua  =  a 

Then,  the  following  relations  between  the  coefficients  are  obtained; 
The  damping  coefficients  are 


(326) 


(327) 


(328) 


(329) 


(330) 


4  4  ”1  =  -(if  Lhx  <-ai 


-mh  -"6 


(331) 
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The  following  relations  Eire  used: 

r/>  bVzAr  =  -(2irpb3Ar)  [-(g)  « 

7r/>bV2Ar  =  -(27rpb3Ar) 


(342) 

(343) 


-L~ 

- 

& 

C/4 

0 


I 

7 


-L 

M 

L  b  Jc/2 


(344) 


V 

= 

"1 

1  *1 

Y 

b 

a 

C/2 

0 

1 

a 

C/4 


Then  the  following  relations  between  the  coefficients  are  obtained: 
The  damping  coefficients  are 


(345) 


k  4 

-mji  “ma 


=  ~(w) 


0 

I 


n 

ka 


II 

m0 


//" 

<b 

n 


rtn 


1 

2 


k  ~  “(2k)  ka 


4  =  ~ (Wk)  ( kb  +  Ika  ) 

(■^R-)(m?  +  ^ka) 

(5k)(mb  +  i[ma  +  kt>]  +(i)  ka  ) 


-me  *  - 


■ma  = 


(346) 

(347) 

(348) 

(349) 

(350) 
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The  stiffness  ocefficients  are 


- 1 

rr 

Q 

_ I 

u 

1 

SB 

!  0 

",  /  j  ~ 

ka  kfa 

1  i 

1 

/  / 

_-mh  -ma 

L 

ma  mb 

0  . 

4  =  “(2)  k°  (352) 

4  =  "'('2)  (kb  +  2  ka  )  (353) 

~mh  =  -(r)(ma  +  2ka)  (35U) 

~ma  =  “(i)(mi  +  rK+  kb]  +  (if  ka  )  (355) 


The  coefficients  in  Reference  8  are  defined  by 


4/>  b  V2k2Ar 


((-I  +  iLg) 

(M,'+iM2) 


(L-3  +  iL-4} 

(mJ+ImJ) 


LE 


(356) 


Note  the  similarity  to  vhe  definition  used  in  Reference  3.  By  use  of  the 
same  procedure  as  that  for  Reference  3,  the  following  relations  between 
the  coefficients  can  be  obtained: 


The  damping  coefficients  are 


4 

4' 

ojJ^ 

If 

1  0 

^2 

l'," 

p  -*1 

_~mK 

“™a_ 

_"2 

Mg 

_°  'J 

(35T ) 


4  =  »  (L.) 

4  =  ^(l,-*l2) 


(358) 

(359) 
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"mh  =  %  (^2-^2) 


(360) 

(361) 


The  stiffness  coefficients  are 


'  -  (#) 

1  0 

L,  Ls 

1  "i 

_4  L 

0  1 

4  ■  ¥  U,) 

4  =  (l3  -  iLi) 

-mh  =  -^  (M,'-  gL, ) 

_ma  =  (M3-  3  [M|  +L3]  +  j-L|  ) 


(362) 

(363) 

(3610 

(365) 

(366) 


Not  only  arc  the  available  tabulated  coefficients  in  differing  notations, 
but  they  also  cover  varying  ranges  of  k  and  .A  computer  program 

was  written  to  transform  the  coefficients  into  the  required  notation  and  to 
assemble  the  coefficients  in  order  of  k  and  MN  .  There  is  a  degree 

of  overlapping  in  the  resulting  tabulation  because,  for  particular  combin¬ 
ations  of  k  and  M  ,  coefficients  are  available  from  more  than  one 
source.  This  overlapping  permits  a  check  to  be  made  of  the  accuracy  of  the 
coefficients. 

From  an  inspection  of  the  table,  it  would  be  seen  that  although  good  agree¬ 
ment  was  obtained  in  the  majority  of  cases,  a  significant  minority  of  cases 
remained  in  which  discrepancies  of  various  magnitudes  existed  between  the 
coefficients.  In  particular,  for  Mn=  0.8,  the  discrepancies  were  large 
and  it  was  not  obvious  as  to  which  were  the  correct  coefficients.  In  an 
attempt  to  sort  out  the  accurate  from  the  inaccurate  coefficients,  it  was 

decided  to  apply  reciprocal  flow  relations  to  all  of  the  available  coef¬ 

ficients. 


6U 


RECII ROCAL  FLOW  RELATIONS  FOR  UNSTEADY ,  COMPRESSIBLE  FLOW 


For  an  airfoil  executing  small  oscillation  in  a  nonstationary  compressible 
inviscid  flow, the  relationship  between  the  pressure  acting  on  the  airfoil 
in  direct  flow  and  the  pressure  acting  on  the  airfoil  when  it  travels  in 
the  reverse  direction  is  given  by  Reference  10.  From  Figure  l4, 

ff  P  "  ds  =  ff  P  w  dS  (367) 

PLANFORM  PLANFORM 


where 

P  is  the  pressure  difference  across  the 
airfoil  in  direct  flow 

P  is  the  pressure  difference  across  the 
airfoil  in  reverse  flow 

w  is  the  vertical  velocity  in  direct  flow 
of  the  air  in  contact  with  the  airfoil 

w  is  the  vertical  velocity  in  reverse  flow 
of  the  air  in  contact  with  the  airfoil 


and  where  Eqs.  (3£8)>  (370),  (372) »  and  (374)  define  direct-flow  quantities 
and  Eqs.  (3 69),  (371),  (373),  and  (375)  are  reverse  flow  quantities. 


and  for  simple  harmonic  motion 


Hence, 


w  =  (vax  +  at  )  zs 

(368) 

w  =  ("v  Ji  *  If ) 

(369) 

motion 

zs  =  Zsae,Wt 

(370) 

zs  = 

(371) 

w  =  wa  eiwt 

(372) 

w  =  wa  e iajt 

(373) 

wa  =  (v^  +  'w)zsa 

(374) 

wa  =  (~v|j  +  i«)  zSa 

(375) 
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in  which  case  the  reciprocal  or  reverse  flow  relation  becomes 


JJ  Pl"a<fc  ‘II  P«, 


ds 


PLANFORM 


(376) 


PLANFORM 


The  application  of  this  relation  to  a  two-dimensional  airfoil  is  now  con¬ 
sidered. 

■Most  unsteady  aerodynamic  theories,  in  particular  those  on  which  the  tab¬ 
ulated  coefficients  are  based,  neglect  airfoil  thickness  effects  and  so 
replace  the  airfoil  by  a  flat  plate  of  zero  thickness.  This,  of  course, 
gives  rise  to  singularities  at  the  leading  and  trailing  edges,  the  singu¬ 
larity  at  the  trailing  edge  being  removed  in  the  usual  manner  by  the  appli¬ 
cation  of  the  Kutta  condition.  With  this  model,  a  two-dimensional  airfoil 
presents  the  same  shape  to  the  airstream  irrespective  of  whether  it  i~  in 
direct  or  reverse  flow.  Hence,  if  the  airfoil  executes  the  same  motion  in 
both  direct  and  reverse  flow,  the  pressure  distribution  ever  the  airfoil 
is  the  same  for  both  flows;  that  is,  p  =  p  ,  and  the  airfoil  will  ex¬ 
perience  the  same  lift  and  moment. 

With  Reference  to  Figure  15 ,  the  following  Eqs.  (377)  through  (381)  are 
obtained  for  direct  flow. 


zs=  -  (hc/4+xa)e,a,t  (377) 

Hence^ 

zsa=  ”  (hc/4  +  *a)  (378) 

(379) 

(380) 

vanish, 

«  -v[(TW(l+T)  °]  (381) 


since  4^  and 

OX  dx 


So  that 


wa  =  (V^  +  '<*>)  ZSQ 

wa  =  -  £  Va  +  ico(  hc/4  +  xce)] 
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With  reference  to  Figure  15,  the  following  equations  (382)  through  (386) 
are  obtained  for  reverse  flow. 


Hence 


so  that 


reverse  flow: 


S  =  “(hC/4-xa)eia,t 

(382) 

;a=  -(hc/4-xa) 

(383) 

1  =  (*“^5x  +  'w)^sa 

(3814) 

=  -  [  Va  +  iu  (  hc/4  -  xa)] 

(385) 

=  -v[(^)hc/4^(l-Wa] 

(386) 

consider  pure  translational  motion  in  direct 

and 

wa  =  -V  (-ty*-)  hc/4 

(387) 

wa  =  -  V  (-jjf4-)  hc/4 

(388) 

when  the  reverse  flow  theorem  for  two  dimensional  flow  is  applied, 


f  Pw0  dx  =  f  P 

/'iirtnn  r. _ 


CHORD 


P  wadx 

CHORD 


(389) 
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Hence , 


i Lp['v(^)i^h 


defining 


L"  =  J  Pdx 

CHORD 


then  the  equation  gives 


Lh 


/> 


P  dx 

CHORD 


Lh  =  Lh 


(390) 


(391) 


(392) 


This  states  that  the  lift  due  to  pure  translational  motion  is  the  same  in 
direct  and  reverse  flow,  a  result  deduced  previously  from  physical  consid¬ 
erations  . 


Pure  translational  motion  in  direct  flow,  and  pitching  about  the  T?%  chord 
point  in  reverse  flow  are  expressed  as  follows: 


wa  = 

-V  (7)hc/4 

(393) 

Wa  = 

-V(l-  Vx)a 

(39U) 

If  Wa  and 

wa  are  substituted  into  the  reverse  flow  theorem 

/  p{-v(i 

WORD 

1  ~  x)a  ]■  dx  = 

/p{-v(Jffl)hc/4}dX 

CHORD 

(395) 

defining 

h 

r 

LC/4  = 

J  Pdx 

(396) 

CHORD 

-a(*<) 

f  _ 

L3C/4  = 

J  Pdx 

CHORD 

(397) 
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(398) 


J  Px  dx 

CHORD 


Then  the  equation  gives 


where 


ik 


-a  ft.) 

^■3  C/4 


h 

LC/4 

h 

^C/4 


is  the  lift  at  25%  chord  due  to  p  ire 
translational  motion  in  direct  flow 

is  the  moment  at  25%  chord  due  to  pure 
translational  motion  in  direct  flow 


-a(*4)  is  the  lift  at  15%  chord  due  to  pitching 
*-3C/4  motion  about  15%  chord  in  reverse  flow 


(399) 


As  shown  previously  for  a  two-dimensional  wing, 

-a(%,)  a{%») 

L3C/4  =  L3C/4 


Hence,  the  reverse-flow  relation  becomes: 


[  k/4  + 


a(^») 

I-3C/4 


(1*00) 


(1*01) 


h  h  . 

Expressions  can  now  be  obtained  for  Lc/jp  ,  Mc/4  and  L3C/4 

in  terms  of  the  aerodynamic  coefficients  which  were  defined  by  the  equation. 


-I 

27r/)b3Ar 


(1*02) 


If  the  relations  for  simple  harmonic  motion  are  applied,  the  equation 
can  be  put  in  the  form 


-I  _ 

2  Tjpv2  bAr 


M 

b 


C/4 


ik4  +  4 

ik(-mfi)-mh 


i  k  /q  *■  la 


h. 

b 


ik(-m6)  -mQ 


(1*03) 
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The  equation  can  now  be  transformed  so  that  it  is  in  terms  of  forces  and 
displacements  at  the  75#  chord  point. 


2irpVl  bAr 


M 

L  b  J3C/4 


I  0 
-I  I 


i  k  4  +  4 

ik(-mt5)-nnh 


ik4  +  4 

ik  (-ma)-ma 


I  - 1 
0  I 


3C/4 


(UoU) 


From  these  two  equations,  it  can  be  shewn  that 
h 


lc/4  =  27r/)VabAr  [ik4  +4]  0^) 

-27rpV2bAr  [ik(-mh)  -mh] 


d(^()  2  r  T 

I-3C/4  *  27r/>V2bAr  [l  0] 


LS='  27r,V2bAr[ik(4-4)  +  (4-4)  ]a 


ik4  +  4 

ik(-mh)-mh  ik(-ma)-ma 


ik4+  4 


lr-in 


(1*05) 

(l*0O 

(U07) 

(1.08) 


Substituting  into  the  reverse-flow  relation  and  dividing  through  by 
2  777)  V2  bAr  gives 


{0*4  +4)  -  ik[ik  (-mh) -mh]  }(%a)a 


=  ik  [  ik  (4 - 4 )  +  (4 " 4)  ](%* 


(1*09) 


Separating  the  equation  into  real  and  imaginary  parts  gives 

{  4  +  k2  [-4  +  4  -  mh  ]j 

+  ik  {  4  -  [  “4  +  4  -  mh  ]  }  =  0 


(1*10) 


Equating  the  real  and  imaginary  parts  to  zero  gives  the  reverse-flow  re¬ 
lation  in  terms  of  the  aerodynamic  coefficients. 


4  =  -4  +  4  ~ 


mi 


4  =  -k2  [-4  +  4-mfi  ] 


(i*n) 


(1*12) 
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It  will  be  noticed  that  the  coefficients  (“ma)  and  (_f^a) 

are  missing  from  the  above  equations, so  that  the  reverse  flow  relation 
provides  no  check  on  the  tabulated  values  of  these  coefficients.  Although 
the  reverse  flow  relation  obtained  above  corresponds  to  that  given  in 
References  11  and  10,  it  was  hoped  that  by  considering  additional  motions 
in  direct  and  reverse  flow,  additional  reverse  flow  relations  containing 
the  missing  coefficients  would  be  obtained.  To  this  end  the  following 
motions  were  considered,  and  the  corresponding  reverse-flow  relations  were 
found  to  be  as  follows: 


Pitching  about  the  25#  chord  point  in  direct  flow; 
pure  translational  motion  in  reverse  flow. 


a(W 

ik  Lc/4 


(*)■  [ 


h 

L3C/4 


+  ik^' 


(M3) 


Pitching  about  the  25#  chord  point  in  direct  flow; 
pitching  about  the  75#  chord  point  in  reverse  flow. 


(MU) 


However,  these  additional  relations  gave  the  same  relations  between  the 
coefficients  as  obtained  before.  The  search  was  not  continued  beyond 
this  stage  but  it  would  be  very  useful  to  obtain  a  reverse-flow  relation 
involving  (-nrijj)  and  (-  ma  )  • 

It  is  possible  to  state  the  reverse-flow  relations  in  the  notations  of  the 
references  from  which  the  tabulated  coefficients  were  obtained.  For 
reference,  they  are  listed  below: 

Reference  2 


ll-  -2k  K"Wa"-mz"] 


(M5) 


(Ul6) 
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Reference  3 


k2  L.',  =  -k  [-2 k2 L2  +  k2t_'4  +  k2M2  ]  (M7) 

k2  L2  =  k  [-2  k2  L7,  +  k2L.3  +  k*  m(  ]  (i»l8) 

Reference  1+ 


Reference  9 


= 

-k  [  ~Lhl  +  Lai  +  Mhx 

(U19) 

*-ht 

^  [“  *“hR  +  LaR  + 

1  (1420) 

Reference  5>  6 

Reference  7 

/ 

ka  = 

.11  II  . 

-k  ( kb  + mQ ) 

(b2l) 

k"  = 

k  (  k^  +  m'a ) 

(1<22) 

Reference  8 

■L,  = 

-  k  [~2LZ  +  L4  +  M2  ] 

(1423) 

*-2  = 

k  [  -2L,  +  l'3  +  M,'  ] 

(U2U) 
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GENERAL  PROCEDURES  FOR  OBTAINING  FLUTTER  SOLUTIONS 


The  following  is  a  brief  summary  of  the  steps  involved  in  obtaining  flutter 
solutions,  and  it  includes  an  outline  of  the  computer  program. 

The  geometry  of  the  blade  is  stated,  including  the  spanwise  breakdown  of 
the  blade  into  segments,  as  required  for  the  numerical  integration.  Then 
the  mass  properties  of  the  fixed  portion  of  the  blade  and  of  the  mass 
balance  weight,  if  any,  are  input.  The  remaining  structural,  data  to  be 
stated  include  such  things  as  the  position  of  the  mid-chord  of  the  blade 
aft  of  a  radial  reference  axis  and  the  chordwise  position  of  the  blade 
elastic  axis  and  center  of  gravity  in  relation  to  the  mid-chord. 

Up  to  10  generalized  coordinates  may  be  chosen  to  describe  the  vibration 
characteristics  of  the  blade.  The  modes  used  as  generalized  coordinates 
are  assumed  to  have  fully  coupled  flatwise,  edgewise, and  torsional  displace¬ 
ments.  These  modes  can  be  elastic  or  rigid-body  modes.  If  modes  are  chosen 
which  have  one  or  two  of  these  freedoms  missing,  the  flatwise,  edgewise, and 
torsional  displacements  must  still  be  specified  for  each  mode,  with  zero 
displacements  being  input  for  the  missing  freedoms.  The  natural  frequency, 
structural  damping  and  generalized  modal  inertia  must  also  be  stated  for 
each  mode. 

All  coefficients  in  the  equations  of  motion  are  calculated  by  the  i>rogram. 
The  structural  and  aerodynamic  forces  on  each  blade  segment  are  calculated 
and  then  summed  across  the  blade  span  to  find  the  generalized  force  acting 
in  each  of  the  generalized  coordinates.  Th>  structural  coefficients  are 
stored  after  being  calculated  and  are  recalculated  only  if  the  mass  data  or 
mode  shapes  are  altered.  However,  the  aerodynamic  coefficients  are  recal¬ 
culated  for  each  eigenvalue  solution.  To  calculate  the  aerodynamic  coef¬ 
ficients,  the  reduced  frequency  k  =  u>b/V  (which  is  based  on  an  assumed 
flutter  frequency  cu  )  and  the  Mach  number  at  each  station  are  found.  Then, 
from  the  table  of  aerodynamic  coefficients  contained  in  the  program,  the 
values  for  the  coefficients  at  each  blade  station  are  found  by  interpola¬ 
tion. 

The  major  variable  in  the  program  is  helicopter  forward  speed.  Up  to  10  veloc¬ 
ities  may  be  input  at  a  time.  The  secondary  ■»,ariables  are: 

1.  Modal  frequency. 

2.  Amount  of  mass  balance. 

3.  Position  of  the  mass  balance  weight. 

.  Aerodynamic  center . 

5.  Elastic  axis  position. 

6.  Assumed  flutter  frequency. 
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7.  The  set  of  generalized  coordinates  used  for  the  solution. 

Item  (7)  needs  some  explanation,  since  generalized  coordinates  are  not  vari¬ 
ables  in  the  same  sense  as  mass  balance.  Suppose  a  set  0  =  (  q ,  qn  ) 

of  generalized  coordinates  is  input,  where  n  <  10  .  Then,  any  subset 

of  Q  may  be  chosen  for  solution,  the  generalized  coordinates  selected 
being  stated  as  part  of  the  input  data.  This  facility  is  used  when  it  is 
desired  to  find  out  the  contribution  tnat  certain  generalized  coordinates 
are  making  to  the  flutter  mechanism.  In  particular,  if  uncertainty  exists 
as  to  how  many  vibration  modes  to  input  to  the  program,  then  extra  modes 
may  be  included.  After  initial  solutions  are  made  to  determine  the  contri¬ 
bution  to  the  flutter  mechanism  of  these  extia  modes,  they  may  be  dropped 
from  all  subsequent  calculations  if  their  contribution  is  found  to  be  small. 

In  addition  to  the  above  variables ,  all  or  any  part  of  the  input  data  may 
be  changed.  However,  the  process  is  not  as  straightforward  as  changing  the 
secondary  variables. 

After  each  eigenvalue  solution  corresponding  to  a  particular  value  of  for¬ 
ward  speed,  the  results  are  printed  out  for  each  eigenvector  in  the  follow¬ 
ing  form: 

Frequency  in  cycles  per  minute. 

Decay  rate  C/Cc  . 

Strength  vectors. 

If  the  decay  rate  C/Cc  is  positive,  then  the  mode  is  stable;  if  negative, 
then  the  mode  is  unstable  and  is  undergoing  divergent  oscillations.  If,  as 
airspeed  increases ,  the  frequency  of  a  mode  drops  to  zero,  this  means  that 
there  is  a  divergent  aerodynamic  force  acting  against  the  elastic  and  cen¬ 
trifugal  restoring  forces  so  that  the  net  restoring  force  is  reduced  to 
zero.  For  higher  airspeeds,  the  modal  frequency  will  remain  at  zero  but 
the  divergent  aerodynamic  forces  are  now  predominant  and  the  mode  suffers 
static  divergence.  The  strength  vectors,  or  eigenvectors,  give  the  relative 
proportion  of  the  generalized  coordinates  in  each  coupled  mode  corresponding 
to  a  particular  eigenvalue,  and  also  the  phase  angle  between  the  generalized 
coordinates.  The  eigenvectors  quite  often  assist  in  the  interpretation  of 
flutter  solutions. 

After  the  eigenvalue  solutions  have  been  obtained  for  the  specified  set  of 
forward  speeds,  a  digital  plot  is  made  of  the  results.  Three  plots  are 
made  in  the  following  order: 

1.  Frequency  to  versus  decay  rate  C/Cc  , 
as  shown  in  Figure  1 6. 

2.  Decay  rate  C/Cc  versus  forward  velocity  VH 
as  shown  in  Figure  17. 


3.  Frequency  u>  versus  forward  velocity  VH 

as  shown  in  Figure  17. 

The  first  plot  is  used  to  identify  the  eigenvalues.  The  eigenvalues  cor¬ 
responding  to  the  first  forward  veiocity  are  plotted  as  1,  those  correspond¬ 
ing  to  the  second  forward  velocity  as  2,  and  so  on  up  to  the  tenth  forward 
velocity  (if  input),  the  corresponding  eigenvalues  being  printed  as  Q 
The  first  plot  is  thus  a  root  locus  plot  with  forward  velocity  as  the 
running  parameter.  The  eigenvalues  are  identified  by  their  values  at  the 
lowest  forward  velocity,  and  then  the  variation  of  the  eigenvalues  with 
increasing  forward  velocity  can  be  traced.  Usually  the  numbers  correspond¬ 
ing  to  a  particular  mode  will  be  grouped  together  as  shown  and  can  be  iden¬ 
tified  by  inspection.  If  ambiguity  exists,  it  can  often  be  resolved  by 
noting  the  variation  of  the  eigenvectors  with  increasing  forward  velocity; 
if  this  fails,  the  ambiguity  can  always  be  resolved  by  obtaining  further 
eigenvalue  solutions,  with  sufficiently  small  increments  in  the  forward 
velocity.  However,  a  reasonable  estimate  of  the  flutter  speed  can  often  be 
obtained  in  spite  of  any  ambiguity,  in  which  case  the  ambiguity  need  not  be 
resolved . 

In  the  second  and  third  digital  plots,  the  eigenvalues  are  plotted  as  an 
asterisk.  When  the  eigenvalues  have  been  identified  by  using  the  first 
plot,  the  corresponding  points  in  the  remaining  two  plots  may  be  identified. 
This  enables  curves  to  be  drawn,  which  give  the  variation  of  decay  rate  and 
frequency  with  changing  forward  velocity  for  each  mode.  Flutter  occurs  if 
a  curve  falls  below  the  axis  C/Cc  =  0  .  Static  divergence  occurs  if  a 
curve  meets  the  axis  oj  =  0  and  C/Cc  <  0 

Sometimes,  two  or  more  eigenvalues  have  the  same  value  within  the  accuracy 
of  resolution  of  the  digital  plot.  In  this  case,  they  are  plotted  as  0. 
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DOCUMENTATION  OF  COMPUTER  PROGRAM 


Items  completing  the  following  documentation  of  the  classical  flutter  com¬ 
puter  program  are  available  at  the  U.  S.  Army  Aviation  Materiel  Labora¬ 
tories,  Fort  Eustis,  Virginia. 

1.  Short  description  of  problem. 

2.  Description  of  solution  method. 

3.  Equations  used  in  computer  solution  are  given  by 
reference  to  this  report. 

4.  Definitions  of  input  symbols  in  program. 

5.  Definitions  of  generated  symbols  in  program. 

6.  Definitions  of  output  symbols  in  program. 

7.  Listing  of  program  source  deck, 
sample  input,  and  sample  output. 

8.  Running  instructions  for  program  and 
estimate  of  running  time. 

9.  Statement  of  type  and  configuration  of  computer  used  in 
the  program  development,  including  any  special  features. 

10.  Name  and  level  of  programming  language  used  in  the  program. 

Copies  of  the  program  source  deck  and  a  sample  input  deck  are  also  available 
at  USAAVLABS. 

The  following  basic  information  is  supplied  here  to  permit  potential  users 
to  evaluate  the  usefulness  of  the  program  for  their  purposes. 

The -program  determines  rotor  blade  flutter  speeds,  frequencies,  and  mode 
shapes ,  using  classical  flutter  theory  in  conjunction  with  the  fixed  azimuth 
assumption. 

The  generalized  structural  inertia,  gyroscopic  coupling,  and  centrifugal 
stiffness  coefficients  are  obtained  for  use  in  the  Lagrange  equations  of 
motion,  with  arbitrary  modes  being  used  as  degrees  of  freedom.  The  gen¬ 
eralized  aerodynamic  coefficients  are  obtained  by  using  standard  classical 
flutter  theory.  The  resulting  equations  of  motion  form  an  eigenvalue 
problem,  which  is  solved  by  Tarnove's  method.  The  eigenvalues  are  complex, 
in  general,  and  define  flutter  frequency  and  modal  damping.  The  condition 
that  modal  damping  is  zero  defines  the  flutter  speeds,  and  the  eigenvectors 
define  the  flutter  mode  shapes. 
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The  program  is  designated  Ei»80  and  is  intended  for  use  on  the  UNIVAC 
1108.  The  source  program  is  written  in  the  FORTRAN  IV  language.  It  is 
divided  into  two  parts,  El*80A  and  El*80B,  which  can  be  executed  independent?.# 
E180A  solves  the  actual  flutter  equations,  while  E^80B  provides  punched  card 
aerodynamic  coefficient  data  for  El*80A.  Each  separate  case  of  EU80B  re¬ 
quires  20  seconds  of  machine  time,  and  EWOA  requires  about  .1  second  for 
each  eigenvalue  solution  with  a  three-degree-of- freedom  case. 


COMPUTER  SOLUTION  TEST  CASE 


A  test  case  was  prepared  to  check  the  computer  solution.  It  is  taken  from 
Reference  12,  and  involves  the  pitch-plunge  flutter  of  a  two-dimensional 
airfoil.  The  equations  of  motion  can  he  put  into  the  form 


where 

qh  =  generalized  coordinate  of  the  mode  consisting 
of  pure  vertical  translation 

=  generalized  coordinate  of  the  mode  consisting 
of  pure  pitch  about  the  pitch  axis 

rar  =  distance  of  the  elastic  axis  aft  of  the  aero¬ 
dynamic  center 

ra  =  radius  of  gyration  about  the  pitch  axis 

wh  =  natural  frequency  of  the  q^  generalized 
coordinate 

wa  =  natural  frequency  of  the  qQ  generalized 
coordinate 


(1*26) 

( U27 ) 

(1428) 


bR  =  b 


( 1*29) 


It  will  be  noticed  that  there  is  no  structural  or  aerodynamic  damping 
present  in  the  equations  of  motion  and  that  the  asymmetry  of  the  stiffness 
matrix  is  clearly  apparent.  This  asymmetry  has  the  potential  to  cause  the 
frequencies  of  the  undamped  vibration  modes  to  coalesce.  Flutter  occurs 
at  the  airspeed  at  which  the  frequencies  of  the  undamped  modes  become 
equal.  The  relationship  of  the  various  axes  is  shown  in  Figure  18.  The 
airfoil  data  are  listed  in  Table  III. 
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Substituting  these  airfoil  data  into  the  equations  of  motion  yields 


1,2516 

-1.8774 


-1.8774 

14.125 


+ 


.57151  xIO'4 
0 


-.47962  xlO'2Z 
.016123(1-  .I79Z) 


(*♦30) 


where  Z  =  .03968  (V/33.8)2  with  V  in  knots. 

The  airfoil  data  were  input  to  the  computer  program,  and  the  above  coef¬ 
ficients  were  compared  with  those  calculated  by  the  program  for  a  range  of 
values  of  V.  In  all  cases,  the  coefficients  were  identical  within  the 
accuracy  of  the  desk  machine  calculations. 

The  eigenvalue  solution  to  the  test  case  is  given  in  Figure  5  of  Reference 
12.  The  results  were  converted  into  the  notation  of  the  computer  program 
output,  and  the  two  methods  of  solution  were  compared  in  Figures  19  and  20. 
Within  the  accuracy  with  which  solutions  could  be  read  from  Figure  5  of 
Reference  12,  the  two  methods  of  solution  give  identical  results.  In 
particular,  the  flutter  and  divergence  speeds  as  quoted  in  Reference  12  are 
compared  in  Table  IV  with  those  obtained  by  the  computer  program. 

To  run  the  test  case  with  the  computer  program,  it  was  necessary  to  use  a 
special  tabulation  of  the  aerodynamic  coefficients.  The  tabulation  con¬ 
sisted  of  the  set  of  coefficients  given  in  Table  V. 

This  set  is  the  same  as  would  be  used  for  a  solution  using  conventional 
aerodynamic  coefficients  at  =  0,  k  =  0  except  that  l\ \  =  0  instead, 

of  1.  This  means  that  the  lift-curve  slope  of  dCL/da  =  2 it  for  a  uni¬ 
form  vertical  translational  velocity  is  replaced  by  zero,  thus  eliminating 
all  aerodynamic  damping.  It  is  use  of  / ii  =  0  which  gives  the  flutter 
solutions  in  Figures  19  and  20  their  unusual  appearance. 

The  test  case  was  also  run  with  the  tabulated  set  of  2320  aerodynamic  coef¬ 
ficients  previously  described .  This  served  as  a  partial  check  on  the  use 
of  this  table  of  coefficients  in  the  computer  program.  The  solution  was 
run  specifying  Mn  =  0,  k  =  0  and  the  results  are  compared  in  Figures  21, 
and  22  with  the  computer  solution  obtained  through  simplified  aerodynamics. 
In  passing,  it  may  be  noted  that  the  simplified  theory  is  not  a  good  approx¬ 
imation  for  this  test  case.  This  is  probably  due  to  the  fact  that  the 
flutter  is  of  the  milder  type  and  that  the  two  modal  frequencies  do  not 
cross . 

This  was  the  only  test  case  in  which  a  check  was  made  of  the  solution  fre¬ 
quencies  and  decay  rates.  However,  many  more  check  cases  were  run  to 
verify  that  the  coefficients  in  the  equations  of  motion  were  being  calcu¬ 
lated  correctly  and  that  the  various  facilities  available  in  the  program 
were  working  correctly.  The  number  of  generalized  coordinates  used  in  these 

checks  varied  from  2  to  10,  the  latter  being  the  maximum  number  for  which 
the  flutter  solution  was  programmed. 
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CORRELATION  OF  CLASSICAL  FLUTTER  CALCULATIONS 
WITH  AVAILABLE  EXPERIMENTAL  DATA 


DISCUSSION  OF  EXPERIMENTAL  DATA  AVAILABLE  FOR  CORRELATION 


Correlation  with  experimental  "blade  flutter  date,  obtained  at  high  Mach 
numbers  and  high  advance  ratios  had  been  intended.  However,  no  such  data 
could  be  located  in  the  l"‘terature  on  the  subject.  Instead,  the  choice 
lay  between  results  obtained  by  NASA  Langley,  Referent..*  13,  for  a  model 
Totor  in  hover  (  /z  =  0}  for  high  Mach  numbers,  and  results  obtained  by 
Cornell  Aeronautical  Laboratory,  Reference  16,  for  a  model  rotor  in  essen¬ 
tially  incompressible  flow  at  high  advance  ratios .  As  the  emphasis  of  the 
study  was  on  the  effect  of  compressibility,  it  was  decided  to  use  the  re¬ 
sults  obtained  by  NASA  Langley  for  hover.  The  Cornell  Aeronautical  Labora¬ 
tory  results  were  used  for  further  correlation. 

Correlation  could  be  attempted  with  only  part  of  the  NASA  Langley  data,  as 
a  large  part  of  it  was  obtained  with  high  pitch  angle  settings,  so  that 
the  blade  was  either  undergoing  stall  flutter  or  at  least  operating  at 
angles  of  incidence  for  which  classical  flutter  theory  is  not  valid. 

CORRELATION  WITH  NASA  LANGLEY  EXPERIMENTAL  RESULTS 

The  correlation  was  divided  into  three  parts,  each  part  using  a  different 
blade  configuration.  The  blades  used  were  designated  NACA  23012  l(R), 

NACA  23012  2(F)  and  NACA  23012  2(R) .  The  main  rotor  and  blade  character¬ 
istics  are  stated  below. 

The  tests  were  conducted  i*.  the  Langley  vacuum  sphere  with  a  single  artic¬ 
ulated  blade.  Thus,  the  blade  operated  at  the  hover  condition.  The  density 
and  speed  of  sound  of  the  testing  medium  could  be  varied.  There  was  no  lag 
hinge  and  no  data  were  given  for  edgewise  motions;  so  the  blade  was  assumed 
to  be  rigid  in  th  •  edgewise  direction.  All  blades  tested,  for  which  corre¬ 
lation  was  attempted,  had  an  NACA  23012  airfoil  section.  The  blades  were 
considered  to  be  rigid  in  flatwise  bending  and  torsion  for  the  inner  8 
inches  of  radius.  The  remaining  38  inches  of  radius  was  considered  as  a 
uniform  cantilever  having  the  properties  given  in  Table  VI. 

Difficulty  was  encountered  in  getting  the  calculated  nonrotating  elastic 
bending  frequencies  to  correlate  with  the  measured  values  stated  in  Table 
VI.  The  cantilever  torsion  frequencies  presented  no  difficulty.  Even¬ 
tually,  it  was  decided  to  increase  the  blade  El  so  as  to  give  a  first 
bending  frequency  equal  to  the  measured  value  and  to  ignore  any  difference 
between  calculated  and  measured  second  bending  frequencies.  The  required 
values  for  El  are  those  stated  in  the  table.  This  approach  was  just¬ 
ified  in  that  flutter  occurred  between  the  flapping  and  torsion  modes,  with 
very  little  participation  of  the  elastic  bending  modes. 

Eefore  proceeding  with  a  discussion  of  the  correlation  results,  mention  will 
be  made  of  the  convention  used  to  refer  to  the  coupled  blade  modes  obtained 
from  the  eigenvalue  solution.  The  labeling  of  the  coupled  modes  follows  the 
conventional  practice  of  referring  to  a  mode  by  the  name  of  the  generalized 
coordinate  with  which  it  becomes  identical  as  the  airspeed  approaches  zero. 
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This  assumes ,  of  course,  that  the  generalized  coordinates  are  normal  modes. 
Otherwise,  the  coupled  mode  will  always  be  a  combination  of  two  or  more  gen¬ 
eralized  coordinates.  As  airspeed  increases,  the  proportion  of  the  general¬ 
ized  coordinates  in  any  coupled  mode  will  vary.  Hence,  a  "torsion"  mode 
will  start  out  as  being  predominantly  torsional;  but  by  the  time  flutter  or 
static  divergence  occurs,  it  may  be  predominantly  flapping.  Hence,  the 
statement  that  torsional  divergence  occurs  is  not  meant  to  imply  that  the 
mode  sh  .pe  was  torsions!  when  the  instability  occurred  but  rather  that  the 
coupled  mode  that  suffered  instability  started  out  at  zero  airspeed  as  be¬ 
ing  a  torsion  mode. 

Correlation  Using  NACA  23012  l(R) 

Flutter  occurred  at  lerf  Mach  numbers.  This  case  was  run  because  results 
were  available  for  very  small  blade  pitch  angles.  Tha  calculated  and  ex¬ 
perimental  results  are  compared  in  Figure  23,  where  both  reduced  flutter 
speed  Vo.8R/bcn  and  the  ratio  of  the  flutter  frequency  to  blade  first 
natural  torsional  frequency  are  plotted  as  a  function  of  blade  pitch  angle. 
The  curve  of  reduced  flutter  speed  separates  the  stable  and  unstable  re¬ 
gions,  with  the  unstable  region  lying  above  the  flutter  curve.  As  the  blade 
pitch  angle  is  increased,  the  reduced  flutter  speed  drops  slightly  at  first 
and  then  rapidly  as  the  blade  apparently  begins  to  stall.  As  the  pitch 
angle  is  further  increased,  the  reduced  flutter  speed  decreases  until  some 
minimum  value  is  reached.  Further  increases  in  pitch  angle  result  in  a 
rather  sharp  rise  in  the  reduced  flutter  speed.  The  curve  of  frequency 
ratio  shows  that  a  reduction  in  the  value  of  the  reduced  flutter  speed  is 
accompanied  by  an  increase  in  flutter  frequency. 

The  portion  of  the  flutter  curve  corresponding  to  low  pitch  angles  defines 
the  region  of  classical  flutter,  ,'heieas  the  portion  of  the  curve  correspond¬ 
ing  to  high  pitch  angles  defines  the  region  of  stall  flutter.  As  shown  by 
the  frequency  ratio  curve,  classical  flutter  occurred  at  a  frequency  lying 
between  the  flapping  and  torsion  mode  frequencies  corresponding  to  the  fact 
that  the  classical  flutter  resulted  from  coupling  of  the  flapping  and  tor¬ 
sion  modes.  On  the  other  hand,  stall  flutter  occurred  at  a  frequency  very 
nearly  equal  to  the  first  torsional  natural  frequency  corresponding  to  the 
fact  that  stall  flutter  is  a  predominantly  torsional  oscillation. 

The  calculated  values  are  valid  for  small  pitch  angles  only,  and  the  results 
were  plotted  for  (#o.8r)s  =  0°  .  It  can  be  seen  that  reasonable  correla¬ 
tion  of  flutter  speed  was  obtained  but  that  the  flutter  frequency  was  con¬ 
siderably  lower  than  the  experimental  values.  This  discrepancy  will  be 
considered  again  later. 

Correlation  Using  NACA  23012  2(R) 

Figure  2k  presents  the  comparison  of  the  experimental  and  calculated  results. 
The  reduced  flutter  speed  now  includes  the  ratio  -Jp7po  which  was  introduced 
to  reduce  the  scatter  of  the  experimental  results  obtained  at  different  den¬ 
sities  of  the  testing  medium.  The  reduction  in  scatter  occurs  because  the 
experimental  flutter  speeds  occurred  at  approximately  constant  dynamic  pres¬ 
sure,  that  is  at  constant  equivalent  airspeed.  The  reduced  flutter  speed  is 
plotted  against  the  blade  tip  Mach  number,  variations  in  tip  Mach  number  be¬ 
ing  obtained  by  varying  the  speed  of  sound  in  the  testing  medium  from  approx¬ 
imately  500  fps  to  1000  fps.  The  calculated  results  show  the  reduced 
flutter  speed  decreasing  slowly  at  first,  then  more  rapidly  as  the  tip 
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Mach  number  increases  from  =  0.2  to  =  0,7 

Then  the  reduced  flutter  speed  decreases  more  slowly  until  a  minimum  occurs 
at  Mj  =  .9,  after  which  it  begins  to  increase  rapidly.  This  beneficial 
effect  of  increasing  Mach  number  is  due  to  the  rearward  shift  in  aerody¬ 
namic  center  as  the  Mach  number  increases  and  can  be  explained  as  follows. 
The  negative  aerodynamic  spring  acting  in  the  torsion  mode  is  composed  ol 
two  parts.  The  first  is  a  function  of  the  aerodynamic  center-to-elastic- 
axis  spacing  and  is  the  usual,  negative  aerodynamic  spring  encountered  in 
fixed-wing  flutter.  The  second  part  is  proportional  to  the  square  of  the 
rotor  RPM  and  is  a  function  of  the  aerodynamic-center-to-center-of -gravity 
spacing  (Reference  1,  pp  153  through  156).  With  the  2(R)  blade,  the 
elastic  axis  is  at  26.5%  chord;  for  low  Mach  numbers,  the  aerodynamic 
center  is  at  25%  chord  sc  that  the  negative  aerodynamic  spring  due  to  the 
aerodynamic-center-to-elastic-axis  spacing  is  negligible.  Hence,  the  only 
negative  aerodynamic  spring  in  the  torsion  mode  is  due  to  the  aerodynamic- 
center  -to-center-of-gravity  spacing.  This  is  large  since  the  center  of 
gravity  is  10.8#  chord  aft  cf  the  25#  chord.  Hence,  as  rotor  RPM  is  in¬ 
creased,  the  frequency  of  the  torsion  mode  decreases  as  can  be  seen  in 
Figure  25*  The  torsion  mode  becomes  coincident  in  frequency  with,  first, 
the  second  flatwise  mode;  then,  the  first  flatwise  mode;  and,  finally,  the 
flapping  mode.  As  the  torsion  mode  frequency  crosses  a  bending  mode  fre¬ 
quency,  it  has  the  opportunity  to  couple  with  that  mode  to  produce  flutter. 
Whether  it  does  or  not  is  a  function,  among  other  things,  of  the  dynamic 
coupling  between  the  bending  and  torsion  modes  and  the  amount  of  damping 
in  each  of  the  modes.  For  uhe  case  under  consideration,  flutter  occurred 
only  when  the  torsion  mode  frequency  crossed  the  flapping  mode  frequency, 
as  can  be  seen  in  Figure  26. 

At  transonic  Mach  numbers,  the  aerodynamic  center  moves  aft,  thereby  de¬ 
creasing  the  aerodynamic-center-to-center-of-gravity  spacing,  and  hence 
reducing  the  negative  aerodynamic  spring  in  the  torsion  mode.  Thus  the 
frequency  of  tne  torsion  mode  decreases  more  slowly  as  a  function  of  rotor  * 
RPM  as  the  aerodynamic  center  shifts  aft.  This  delays  to  higher  values  of 
rotor  RPM  the  point  at  which  frequency  coincidence  and,  hence,  flutter 
occur.  As  long  as  the  center  of  gravity  is  not  too  far  aft,  a  Mach  number 
will  eventually  be  reached  where  the  torsion  mode  frequency  no  longer  de¬ 
creases  with  increasing  rotor  RPM  and  may  even  begin  to  increase.  In 
either  ease,  if  the-  torsion  mode  frequency  has  not  come  close  to  the  flap¬ 
ping  mode  frequency,  then  flutter  is  eliminated.  From  Figure  2h ,  it  can 
be  seen  that  the  beneficial  effect  of  high  Mach  numbers  was  much  more  pro¬ 
nounced  with  the  experimental  results  than  with  the  calculated  values;  this 
indicates  that  perhaps  the  aerodynamic  center  moved  aft  faster  on  the  test 
blade  than  was  implied  by  the  theoretical  aerodynamic  coefficients. 

In  passing,  it  can  be  stated  that  the  above  mechanism  by  which  coupled 
bending-torsion  flutter  occurs  with  a  rotating  blade  is  quite  general.  The 
mechanism  can  be  described  in  the  following  way.  A  negative  aerodynamic 
spring  must  exist  in  the  torsion  mode  to  cause  its  frequency  to  decrease 
with  increasing  rotor  RPM.  A  negative  aerodynamic  spring  will  result  when 
the  elastic  axis,  or  the  center  of  gravity,  or  both,  are  aft  of  the  aero¬ 
dynamic  center.  As  the  torsion  mode  frequency  decreases,  it  will  cross  the 
bending  modes  from  the  highest  to  the  lowest  frequencies.  With  each  cross- 
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ing,  there  is  the  possibility  that  flutter  will  occur;  and  with  each 
successive  crossing,  this  possibility  increases  because  the  potential  to 
flutter  is  greatest  w'-th  the  low-frequency  bending  modes.  Also,  there  is 
quite  a  chance  that  after  the  torsion  mode  has  crossed  the  lowest  bending 
(or  flapping)  mode,  its  frequency  will  drop  to  zero,  at  which  point  static 
divergence  occurs ,  as  shown  in  Figure  25 •  Thus ,  if  flutter  occurs ,  there 
is  a  good  chance  that  static  divergence  will  occur  at  a  slightly  higher 
rotor  RPM.  In  fact,  if  flutter  does  not  occur  until  after  frequency  coin¬ 
cidence  has  occurred,  as  is  sometimes  the  case,  then  static  divergence  may 
occur  before  flutter. 

Returning  to  Figure  2k,  it  can  be  seen  that  reasonable  correlation  of  . 
flutter  speed  occurred  for  the  lower  blade  pitch  angle  (namely,  (#o.8r)s 
=  7.2°)  but  that  the  calculated  flutter  speed  is  unconservative  when  com¬ 
pared  with  the  experimental  results  for  the  larger  pitch  angle  (namely, 
(@0.8r)s  =  11.2°).  The  drop  in  the  measured  flutter  speed  at  the  higher 

pitch  angles  is  consistent  with  the  results  given  in  Figure  23  and  would 
indicate  that  the  flutter  mechanism  is  changing  from  classical  flutter  to 
stall  flutter.  The  frequency  curve  in  Figure  2k  shows  that  for  both  values 
of  blade  pitch  angle,  the  flutter  frequency  was  approximately  equal  to  the 
nonrotating  blade  torsion  frequency,  which  indicates  that  perhaps  the 
flutter  mechanism  for  the  instability  at  the  lower  pitch  angle  was  also 
more  in  the  nature  of  stall  flutter. 

Correlation  Using  i?ACA  23012  2(F) 

Blade  2(F)  is  very  similar  to  blade  2(R)  considered  in  the  previous  sec¬ 
tion.  The  main  difference  is  that  the  center  of  gravity  has  been  moved 
forward  so  that  it  is  2. 5%  chord  aft  of  the  elastic  axis  instead  of  10.8# 
chord.  The  comparison  of  experimental  and  calculated  flutter  speeds  is 
given  in  Figure  27*  Comparing  the  results  for  blade  2(F)  with  those  for 
blade  2(R)  given  in  Figure  2k  shows  that  decreasing  the  center-of-gravity 
offset  has  caused  the  reduced  flutter  velocity  for  the  experimental  res  llts 
to  increase  slightly, but  that  the  reduced  flutter  velocity  for  the  calcu¬ 
lated  results  has  almost  doubled.  Again,  the  measured  flutter  frequencies 
are  considerably  higher  than  the  calculated  values.  The  most  disturbing 
aspect  of  the  lack  of  correlation  is  that  the  predicted  flutter  speed  is 
considerably  higher  than  the  measured  values,  that  is,  the  predicted  flutter 
speed  is  unconservative.  The  discrepancies  between  the  measured  and  calcu¬ 
lated  values  as  outlined  above  could  be  attributed  to  the  model  blade's 
experiencing  stall  flutter,  since  the  transition  from  classical  to  stall 
flutter  is  marked  by  the  following  characteristics: 

1.  The  sensitivity  of  flutter  speed  to  center-of-gravity 
offset  decreases  and  practically  disappears. 

2.  The  flutter  frequency  increases  considerably. 

3.  The  flutter  speed  decreases  considerably. 

Items  2.  and  3.  can  be  seen  in  Figure  23. 


83 


Satisfactory  correlation  of  reduced  flutter  speed  was  obtained  for  blades 
l(R)  and  2(.R),  both  of  which  had  a  large  center-of-gravity  offset.  The 
lack  of  correlation  of  the  flutter  frequency  for  blades  2(R)  and  2(F)  and 
the  unconservative  reduced  flutter  speed  calculated  for  blade  2(F)  could 
be  explained  by  assuming  that  the  flutter  mechanism  was  stall  flutter 
rather  than  classical  flutter.  Such  an  explanation,  however,  cannot  be 
used  to  explain  the  lack  of  correlation  of  flutter  frequency  when  consid¬ 
ering  blade  l(R).  Hence,  it  was  decided  to  attempt  to  correlate  with  the 
experimental  flutter  results  stated  in  Reference  lU. 

CORRELATION  WITH  CORNELL  AERONAUTICAL  LABORATORY  EXPERIMENTAL  RESULTS 

The  test  equipment  consisted  of  a  single-blade  hydraulically  driven  rotor 
mounted  above  a  test  trailer  which  was  towed  along  an  airport  taxistrip. 
Blade  chordwise  center-of-gravity  position  was  adjusted  by  means  of  lead 
weights  housed  in  a  cylindrical  pod  at  the  blade  tip.  Two  features  of  the 
blade  were  somewhat  unusual.  One  was  the  location  of  the  blade  elastic 
axis  at  hl%  chord.  This  gives  a  large  negative  aerodynamic  spring  in  the 
torsion  mode  due  to  the  large  aerodynamic-center- to-elastic-axis  spacing. 

In  this  respect,  the  blade  is  more  akin  to  a  fixed  wing.  The  other  feature 
was  the  very  low  frequency  of  the  rigid  body  pitch  mode.  This  frequency 
was  varied  during  the  testing.  Otherwise,  the  Cornell  model  blade  was 
similar  to  the  NASA  Langley  model  blade.  The  blade  was  untwisted  and  ar¬ 
ticulated  in  flapping  with  a  zero  offset  flapping  hinge.  There  was  no  lag 
hinge  and  since  no  data  were  given  in  the  edgewise  direction,  the  blade 
was  assumed  to  be  rigid  for  edgewise  deformations.  With  respect  to  flat¬ 
wise  bending  and  torsion,  the  blade  was  also  assumed  to  be  rigid  out  to  a 
blade  radius  of  7  inches.  The  characteristic  blade  parameters  are  given 
in  Table  VII. 

From  the  blade  characteristics  it  appears  that  the  blade  had  small  center- 
of-gravity  offset,  as  the  center  of  gravity  was  only  1.5%  chord  aft  of 
elastic  axis.  However,  the  pitch  frequency  was  much  lower  than  the  canti¬ 
lever  torsion  frequency,  so  that  flutter  occurred  as  a  result  of  coupling 
between  the  flapping  and  rigid-body  pitch  modes.  Hence,  for  the  unbalanced 
blade,  with  the  feathering  axis  at  25%  chord,  the  center-of-gravity  offset 
was  16%  chord,  which,  in  combination  with  the  low  pitch  frequency  and  the 
large  negative  aerodynamic  spring  in  the  pitch  mode,  resulted  in  flutter 
at  low  forward  speeds  and  rotor  RPM.  Because  of  the  very  low  rotor  rota¬ 
tional  speeds  employed,  flutter  was  obtained  up  to  /x  =  0.7  at  relatively 
low  forward  velocities,  with  a  blade  tip  Mach  number  always  less  than  0.3. 

Correlation  was  attempted  with  the  results  given  in  Figures  10(a),  10(b), 
and  10(c)  of  Reference  1 6.  The  structural  parameters  varied  were  blade 
pitch  frequency  and  blade  chordwise  center-of-gravity  location,  the  latter 
being  varied  by  means  of  balance  weights  added  to  a  pod  at  the  blade  tip. 

A  diagram  of  the  pod  and  the  balance  weight  positions  is  given  in  Reference 
lU.  The  parameters  corresponding  to  the  correlation  cases  are  given  in 
Table  VIII. 


Qh 


where 


®*t 

“do 


is  the  nonrotating  first  flatwise  bending  frequency 

is  the  nonrotating  rigid-body  pitch  frequency 

is  the  distance  of  the  balance  weights  ahead  of 
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Xe  is  the  equivalent  center-of-gravity  location.  It  is 

c  defined  in  Reference  I**  as  the  chordwise  center  of 

gravity  location  of  a  rectangular  blade  having  a  uniform 
mass  distribution  such  that  the  ratio  of  the  pitching¬ 
flapping  product  of  inertia  to  the  flapping  inertia  is 
the  same  as  that  of  the  nonuniform  blade  being  studied. 
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The  experimental  d&ta  of  Figure  10(c)  in  Reference  16  is  for  the  case  of 
no  added  ’.eights,  and  represents  flutter  of  the  uniform  unbalanced  blade. 

The  results  given  in  Figure  28  pertain  to  data  for  the  flatwise  to  torsion 
frequency  ratio  aJ<£,/d ig0  =  .63,  with  uJ0o  =  132.5  radians/sec.  The 
rotor  speed  at  which  flutter  occurred  is  denoted  by  SI  ,  and  uJq0  is 
the  nonrotating  rigid-body  pitch  frequency.  Hence,  raising  the  flutter 
speed  Si  corresponds  to  decreasing  the  value  of  TDq0/SI  .  Nondimen- 
sionalizing  the  flutter  speed  SI  with  respect  to  the  pitch  frequency 
U)g0  allows  experimental  results  obtained  at  different  values  of  wq0 
to  be  plotted  together,  as  the  flutter  speed  is  roughly  proportional  to 
cn0o  .  It  can  be  seen  that  the  predicted  flutter  speeds  agree  very  well 
with  the  measured  values.  The  calculated  flutter  speed  corresponding  to  a 
particular  value  of  /z  was  obtained  by  simultaneously  increasing  forward 
speed  and  rotor  speed  so  as  to  keep  /u  constant.  For  example,  the  re¬ 
sults  corresponding  to  p.  -  0  are  given  in  Figures  )  and  30.  Figure  29 
is  a  plot  of  the  coupled  modal  frequencies  versus  Si  ,  and  Figure  30  is 
a  plot  of  the  corresponding  modal  decay  rates.  The  decay  rate  of  the  flap¬ 
ping  mode  changes  at  SI  =  3^5  RPM  from  being  convergent  to  divergent, 
giving  a  flutter  speed  of  3^5  RPM.  The  results  are  of  particular  interest, 
as,  from  the  frequency  and  damping  curves,  it  would  appear  that  the  flutter 
occurs  between  the  flapping  and  first  bending  modes.  This  was  the  only 
case  in  which  this  occurred,  and  it  was  also  the  only  case  for  which  the 
pitch  frequency  was  above  the  first  flatwise  bending  frequency.  A  detailed 
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investigation  was  made  of  this  case,  from  which  it  was  deduced  that: 

1.  The  only  modes  to  couple  together  to  give  flutter  were 
flapping  and  rigid-body  pitch.  The  resulting  flutter 
frequency  and  decay  rate  were  a  good  approximation  to 
the  five  mode  solution. 

2.  The  addition  of  first  flatwise  bending  to  the  flapping- 
pitch  binary  gave  a  solution  almost  identical  to  the 
five  mode  solution. 

3.  The  second  flatwise  bending  mode  and  cantilever  torsion 
modes  made  a  negligible  contribution  to  the  flutter 
mechanism  and  could  have  been  deleted  from  the  analysis. 

Thus,  flutter  was  due  to  coupling  of  the  flapping  and  rigid-body  pitch 
modes  as  was  to  be  expected,  and  this  was  found  to  be  the  flutter  mechanism 
for  all  of  the  cases  that  were  studied. 

The  results  given  in  Figure  31  pertain  to  data  for  the  lower  pitch  fre¬ 
quency;  they  can  be  compared  directly  with  Figure  10(c)  of  Reference  16. 

The  pitch  frequency  of  the  blade  has  been  reduced  from  132.5  rad/sec  to 
64.5  rad/sec  so  that  the  nonrotating  pitching  frequency  is  now  lower  than 
the  nonrotating  first  flatwise  bending  frequency,  with  the  uJ^/tugo  =  1.31. 
The  variations  of  the  coupled  modal  frequencies  and  decay  rates  with  in¬ 
creasing  airspeed  and  for  fj.  =0.5  are  given  in  Figures  32  and  33 

respectively.  These  variations  are  typical  of  the  variations  in  coupled 
modal  frequency  and  decay  rate  that  were  obtained  in  all  the  remaining 
cases.  It  can  be  seen  that,  due  to  the  large  negative  aerodynamic  spring 
in  the  pitch  mode,  the  pitch  mode  frequency  decreases  rapidly  with  in¬ 
creasing  airspeed  or  ii  .  Flutter  occurs  near  the  airspeed  at  which 
the  pitch  mode  frequency  crosses  the  flapping  mode  frequency.  It  can  be 
seen  that  the  flutter  is  extremely  violent,  the  decay  rate  of  the  flapping 
mode  changing  from  C/Cc  =  0.30  (stable)  to  C/Cc  =  0.0  (flutter) 

with  an  increase  in  rotor  speed  of  only  13  RPM.  This  is  to  be  expected 
with  pitch-flap  flutter,  especially  with  the  large  aerodynamic-center-to- 
center-of-gravity  and  pitch-axis-to-center-of-gravity  offsets.  After  fre¬ 
quency  coalescence  of  the  pitch  and  flapping  modes ,  the  frequency  of  the 
pitch  mode  continues  to  decrease  until  flapping  divergence  occurs.  In 
some  of  the  other  cases,  the  frequency  of  the  pitch  and  flapping  modes 
remained  almost  equal  after  flutter  occurred,  and  static  divergence  did 
not  take  place. 

In  Figure  31,  it  will  be  noted  that  the  calculated  results  are  now  somewhat 
conservative;  that  is,  the  calculated  flutter  speed  is  lower  than  the  ex¬ 
perimental  flutter  speed.  Also  plotted  is  the  flutter  speed  calculated  in 
Reference  16.  It  can  be  seen  that  the  calculated  flutter  speed  of  the 
present  study  lies  between  the  measured  value  and  the  value  calculated  in 
Reference  16. 

The  experimental  data  of  Figure  10(b)  in  Reference  16  is  for  the  partially 
balanced  case,  with  the  results  shown  in  Figure  34.  The  calculated  flutter 
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speed  is  now  lower  t..an  both  the  measured  value  and  the  value  calculated 
by  Cornell.  For  higl  advance  ratios,  both  calculated  values  are  extremely 
conservative  with  respect  to  tl.e  measured  flutter  speed.  For  example,  at 
f±  =0.5,  the  rotational  flutter  speed  calculated  in  Reference  14  is 
.50  times  the  measured  rotational  flutter  speed.  The  calculations  of  the 
present  study  give  a  rotational  flutter  speed  139  times  the  measured  rota¬ 
tional  flutter  speed 

The  Reference  16  results  were  obtained  by  using  an  analogue  solution,  and 
the  stability  of  the  rotating  blade  was  considered;  whereas  the  present 
study  uses  the  stopped-rotor  approach,  the  solutions  being  obtained  with 
the  rotor  stopped  at  an  azimuth  angle  of  \j/  =  90°,  which  is  considered 
to  be  the  most  unstable  azimuth  angle  with  respect  to  classical  flutter. 
Hence,  there  is  an  inherent  degree  of  conservatism  in  the  stopped  rotor 
approach;  ideally,  solutions  should  be  obtained  for  different  positions 
around  the  azimuth  angle  with  the  resulting  maximum  and  minimum  flutter 
speeds  bounding  the  actual  flutter  speed.  For  example,  the  region  of  in¬ 
stability  may  extend  over  an  azimuth  angle  of  60°.  While  the  blade  is 
passing  through  this  region,  its  oscillations  will  be  diverging.  For  the 
remaining  300°  of  rotation,  they  will  be  dying  away.  The  flutter  speed 
would  be  obtained  by  integrating  the  buildup  and  decay  of  the  oscillations 
for  one  revolution,  flutter  occurring  when  the  oscillation  amplitude  at 
the  end  of  the  360°  rotation  was  equal  to  the  amplitude  at  the  beginning. 
However,  this  approach  was  beyond  the  scope  of  the  study. 

Referring  again  to  Figure  34,  one  can  see  tlxvt  the  measured  flutter  speed 
occurred  at  a  constant  value  of  rotor  opeed  (for  uJ6o  constant)  inde¬ 
pendent  of  the  advance  ratio.  That  is,  the  measured  flutter  rotor  speed 
did  not  vary  with  forward  velocity,  which  is  a  very  surprising  result. 

The  coefficients  in  the  flutter  equations  are  obtained  by  integrating 
blade  structural  and  aerodynamic  properties,  including  the  mode  shapes, 
across  the  span.  With  a  fixed  wing,  this  integration  or  averaging  process 
for  a  simple  planform  results  in  a  good  definition  of  flutter  speed  through 
use  of  the  wing  properties  at  the  75%  span.  A  similar  effect  would  be  ex¬ 
pected  for  a  rotating  blade  for  which,  in  addition,  the  airspeed  varies 
spanwise.  The  flutter  speed  of  a  rotating  blade  would  then  be  a  function 
of  the  airspeed  at  some  spanwise  station,  say  Q0%  span,  and  would  be  in¬ 
dependent  of  the  combination  of  rotor  speed  and  forward  velocity  giving 
the  velocity  at  the  station.  Both  sets  of  calculated  results  behave  in 
this  manner.  As  forward  speed  increases,  the  rotor  speed  at  which  flutter  . 
occurs  decreases,  so  that  flutter  occurs  at  the  same  value  of  average  blade 
airspeed.  However,  the  measured  flutter  speeds  in  Figure  34  seem  to  in¬ 
dicate  that  the  flutter  mechanism  depends  predominantly  on  rotor  speed  and 
is  not  influenced  by  changes  in  the  average  airspeed  of  the  blade.  This 
trend  of  the  experimental  results  certainly  deserves  to  be  investigated 
thoroughly,  as  it  seems  to  indicate  that  some  aspect  of  the  flutter  mech¬ 
anism  is  missing  from  the  theoretical  flutter  calculations. 

The  results  given  in  Figure  35  pertain  to  the  almost  completely  balanced 
blade  of  Figure  10(a)  in  Reference  16.  The  flatwise  to  torsion  frequency 
ratio  is  oj^,/ajg0  =  1.03,  and  the  torsion  natural  frequency  is 

=  65.9  rad/sec.  Again,  the  test  results  show  that  flutter 
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occurred  at  constant  rotor  speed.  The  digital  solutions  are  so  conserva¬ 
tive  as  to  bear  little  relation  to  the  measured  results.  Flutter  still 
consists  of  pitch  flap  coupling,  and  the  question  arises  as  to  why  flutter 
occurred  between  flapping  and  pitch  when  the  inertia  coupling  between  these 
modes  is  almost  zero.  If  one  looks  at  the  change  in  the  measured  flutter 
speed  with  changing  mass  balance  as  given  in  Table  IX,  it  can  be  seen  that 
the  measured  flutter  speed  increased  by  a  factor  of  4  from  the  unbalanced 
to  the  fully  balanced  case,  whereas  the  calculated  flutter  speeds  increased 
only  slightly.  It  is  this  insensitivity  of  the  calculated  results  to 
changing  mass  balance  that  makes  the  correlation  so  poor  for  the  fully 
balanced  blade.  In  Figures  34  and  35  the  results  have  been  plotted  for 
the  case  where  the  mass  balance  coefficients  were  not  added  to  the  equa¬ 
tions  of  motion.  Again  the  insensitivity  of  the  calculated  flutter  speed 
to  changing  mass  balance  is  apparent. 

The  results  given  in  Figure  36  pertain  to  the  balanced  blade  of  Figure 
10(a)  in  Reference  1 6  with  the  flatwise  to  torsion  frequency  ratio 

/d>0o  =  1.44  and  t' natural  torsion  frequency  wgQ  =  46.9  rad/sec. 
Again  the  insensitivity  of  the  calculated  flutter  speed  to  changing  mass 
balance  has  caused  the  calculated  flutter  speed  to  be  extremely  conserva¬ 
tive. 

SUMMARY  OF  CORRELATION  WITH  NASA-I mGLEY  AND  CORNELL  TEST  RESULTS 


The  correlation  with  Cornell  test  results  was  begun  because  of  the  poor 
correlation  of  flutter  frequency  that  occurred  with  the  NASA-Langley  model 
results.  The  frequencies  at  which  the  Cornell  model  blade  fluttered  are 
not  stated,  and  so  a  comparison  of  calculated  and  measured  flutter  fre¬ 
quencies  cannot  be  made.  However,  it  was  hoped  to  achieve  satisfactory 
correlation  of  flutter  speed,  which  is  much  more  important  than  flutter 
frequency.  From  the  above  review  of  the  correlation  with  Cornell  data,  it 
can  be  seen  that  this  was  not  achieved.  The  overall  correlation  is  cer¬ 
tainly  encouraging,  but  in  some  areas  questions  remain  to  be  answered.  For 
cases  where  it  is  reasonably  certain  that  the  model  blade  was  undergoing 
classical  flutter,  the  calculated  flutter  speed  was  in  every  case  either 
equal  to  or  below  the  measured  flutter  speed.  This  conservative  nature  of 
the  calculated  flutter  speed  is  important.  It  may  cause  concern  unnec¬ 
essarily,  but  it  gives  assurance  that  classical  flutter  will  not  occur 
below  the  calculated  speed. 

The  emphasis  of  the  study  was  on  the  effects  of  compressibility  at  the 
high  advance  ratios  at  which  present-generation  uelicopters  are  operating. 
It  was  for  this  purpose  that  the  compressible  aerodynamic  coefficients 
were  assembled  and  used  in  the  computer  program  for  obtaining  flutter  solu¬ 
tions.  However,  because  of  the  lack  of  test  data,  the  correlation  had  to 
be  made  with  cases  involving  either  high  Mach  numbers  in  hover  or  very  low 
Mach  numbers  in  hover  and  in  forward  flight.  In  fact,  quite  a  few  of  the 
cases  were  for  the  hover  condition  involving  low  tip  Mach  numbers  for  which 
a  satisfactory  unsteady  aerodynamic  theory  exists  (Reference  15  ) •  It  is 
expected  that  use  of  the  Loewy  coefficients  in  these  cases  would  improve 
the  correlation  obtained. 


88 


It  is  unfortunate  that  test  results  were  not  available  for  the  high  Mach 
number,  high  advance  ratio  flight  regime;  not  only  is  this  the  area  where 
a  proven  method  is  needed  to  predict  classical  flutter  speeds,  but  it  is 
also  the  area  where  the  compressible  aerodynamic  coefficients  used  in  the 
computer  solution  are  considered  to  be  applicable.  On  the  other  hand,  the 
validity  of  these  coefficients  for  the  hover  case  with  low  tip  Mach  numbers 
is  questionable. 

It  is  felt  that  the  methods  developed  during  the  present  study  can  be  used 
as  a  design  tool  in  evaluating  blade  flutter  speeds  and  in  determining  the 
sensitivity  of  the  flutter  speed  to  changes  in  design  parameters.  It  is 
suggested  that  the  compressible  aerodynamic  coefficients  be  replaced  by 
Loewy  coefficients  for  cases  involving  low  tip  Mach  numbers  with  low  ad¬ 
vance  ratios. 


EFFECTS  OF  VARIATIONS  IN  PARAMETERS 


DESCRIPTION  OF  BASIC  CONFIGURATION 


It  was  decided  to  use  the  NACA  23012  2(F)  model  blade  described  in 
Reference  13  for  the  parametric  variations,  with  the  same  variations  being 
considered  for  the  two  cases  where  the  blade  root  is  considered  fixed  and 
articulated.  The  rotor  speed  and  chordwise  center-of-gravity  position 
were  adjusted  to  give  reasonably  strong  flutter  at  about  fx  =0.25. 
Otherwise,  if  flutter  is  weak,  small  changes  in  parameters  are  likely  to 
have  a  very  stabilizing  effect;  whereas  with  violent  flutter,  large  changes 
in  the  parameters  may  have  little  effect  in  reducing  the  flutter.  Initial 
solutions  attempted  to  adjust  the  speed  of  sound  so  that  at  fx  =  .50, 
the  Mach  number  was  unity.  This  is  equivalent  to  testing  the  model  blade 
in  a  gas  such  as  Freon,  for  which  the  speed  of  sound  is  much  less  than  in 
air.  However,  varying  the  speed  of  sound  had  a  confusing  effect  on  the 
flutter  speed,  so  this  attempt  was  not  pursued  and  the  solutions  were  run 
with  the  speed  of  sound  equal  to  1116  ft/sec.  Therefore,  the  blade  tip 
Mach  number  at  which  flutter  occurred  was  relatively  low.  Air  density  was 
set  at  its  standard  sea  level  value.  All  cases  studied  were  for  the  ad¬ 
vancing  blade. 

The  main  parameters  defining  the  nominal  case  for  the  articulated  and 
hingeless  blades  are  given  in  Table  X. 

The  blade  is  considered  to  be  rigid  for  the  first  8  inches  out  from  the 
rotation  axis.  For  the  remaining  38  inches,  the  blade  is  uniform  and  the 
structural  data  refer  to  the  uniform  part  of  the  blade.  The  parameter 
variations  that  were  considered  are  summarized  in  Table  XI.  In  the  case 
numbers,  A  refers  to  the  articulated  blade  and  R  to  the  nonart iculaced 
(rigid)  blade.  It  will  be  noticed  that  each  parameter  is  varied  individ¬ 
ually  from  the  nominal  case;  simultaneous  variations  are  not  considered. 

The  variations  for  the  articulated  blade  are  considered  first. 

The  results  are  presented  in  the  form  of  stability  boundaries  such  as  those 
shown  in  Figure  39*  The  advance  ratio  fx  at  which  instability  either 
begins  or  ends,  is  plotted  versus  the  parameter  under  consideration.  For 
a  particular  curve,  the  unstable  region  is  always  above  the  curve  except 
for  the  ct.  e  where  a  ioop  occurs.  In  this  case,  the  unzcable  region  lies 
inside  the  loop  so  that  as  advance  ratio  increases ,  instability  occurs  on 
crossing  the  lower  boundary  and  stability  is  regained  on  crossing  the  upper 
boundary.  The  stable  region  of  interest  is  the  area  lying  below  the  flutter 
curves  and  is  that  region  in  which  a  practical  helicopter  blade  must  oper¬ 
ate.  The  u^per  boundary  of  this  region  is  shown  shaded  in  the  diagrams. 
Some  of  the  curves  are  drawn  in  as  broken  lines  to  indicate  that  only  the 
form  of  the  variation  is  known,  not  the  exact  shape  and  position  of  the 
curves.  The  precise  location  of  the  curves  can  always  be  determined  by 
obtaining  additional  flutter  solutions.  This  was  not  attempted,  as  all 
information  required  from  the  parametric  variation  is  contained  in  the 
broken  curves. 
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RESULTS  FOR  THE  ARTICULATED  BLADE 

For  reference,  the  variation  of  the  solution  frequency  and  decay  rate  ver¬ 
sus  rotor  speed  is  gx/en  in  Figures  37  and  SS^espectivelyj  for  the  nominal 
blade  operating  in  hover.  The  blade  flutters  at  1190  RPM,  at  which  speed 
the  torsional  mode  and  flapping  mode  frequencies  are  almost  equal. 

Cases  1.A1  and  1.A2  of  Table  XI  consider  two  elastic  axis  positions  in 
addition  to  the  nominal  position;  the  results  are  shown  in  Figure  39* 

Figure  39  also  shows  a  number  of  additional  results ,  which  were  obtained 
to  permit  better  definition  of  the  instability  boundaries.  With  the 
elastic  axis  at  25#  chord,  the  major  contribution  to  the  negative  aero¬ 
dynamic  stiffness  is  the  aerodynamic-center-to-center-of-gravity  spacing, 
and  this  contribution  is  unchanged  as  the  elastic  axis  is  shifted  aft. 
However,  moving  the  elastic  axis  toward  the  blade  trailing  edge  increases 
the  aerodynamic-center-to-elastic-axis  spacing,  which  begins  to  make  a 
sizeable  contribution  to  the  negative  aerodynamic  stiffness,  thereby  caus¬ 
ing  a  reduction  in  the  torsional  divergence  speed.  The  reduction  in  the 
elastic-axis-to-center-of-gravity  spacing  from  13#  to  3#  chord  and  finally 
to  -7#  chord  (center  of  gravity  ahead  of  the  elastic  axis)  has  eliminated 
flutter  of  the  first  flatwise  mode,  raised  the  flutter  speed  of  the  flap¬ 
ping  mode,  and  eliminated  flutter  below  the  torsional  divergence  speed. 

As  the  elastic  axis  moves  from  37*1#  chord  to  37.2#  chord,  the  roles 
played  by  the  flapping  and  torsional  modes  are  interchanged,  but  there  is 
no  discontinuity  in  the  curves  as  a  result  of  this  switch.  On  the  root 
locus  plots  generated  by  the  program,  the  flapping  and  torsional  modes 
approach  each  other  when  the  advance  ratio  increases,  and  for  some  special 
value  of  elastic  axis  position  between  37-1%  and  37*2#  chord,  they  will 
meet  tangentially  from  opposite  directions.  At  this  point  the  flapping 
and  torsional  modes  have  identical  frequency,  damping,  and  motion.  Further 
increases  in  the  advance  ratio  cause  the  root  locus  paths  to  leave  the 
point  of  intersection  at  right  angles  to  the  incident  paths.  It  is  then 
arbitrary  as  to  which  of  the  retreating  paths  is  termed  the  flapping  or 
torsional  mode.  As  the  advarce  ratio  is  increased,  one  of  the  modes 
reaches  a  divergent  instability  and  then  the  other  reaches  a  flutter  in¬ 
stability.  The  onset  of  these  instabilities  provides  the  points  on  the 
boundaries  in  Figure  39  where  the  flutter  instability  and  the  divergence 
instability  modes  interchange.  With  the  elastic  axis  forward  or  aft  of 
this  special  position,  the  flapping  and  torsional  branches  each  break  away  ■ 
without  reaching  a  point  of  intersection.  In  this  case,  the  modes  can,  as 
usual,  retain  the  designations  assigned  for  fj.  =  0.  For  forward  positions 
of  the  elastic  axis  it  is  the  torsional  mode  which  meets  the  a)  -  0  axis 
and  therefore  suffers  static  divergence,  whereas  for  aft  positions  it  is 
the  flapping  mode  which  meets  the  =  0  axis. 

The  switch  in  the  stability  of  the  torsional  and  flapping  modes  is  a  result 
of  the  change  in  the  modal  coupling  that  occurs  as  the  elastic  axis  is 
shifted  aft.  This  effect  is  examined  in  Figures  40  and  4l,  in  which  the 
amplitude  ratio  qh  / c\a  and  the  phase  angle  are  plotted  versus  advance 
ratio,  for  both  the  flapping  and  the  torsional  modes.  Figure  !»0  considers 
the  elastic  axis  at  25#  chord,  and  it  can  be  seen  that  in  hover,  the  fla^- 
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ping  mode  begins  as  a  predominantly  flapping  motion  qh  /qa  =  16  while 
the  torsional  mode  beings  as  a  predominantly  torsional  motion  qh/qa  =0.6. 
As  advance  ratio  increases,  the  amount  of  torsional  motion  in  the  flapping 
mode  increases,  and,  simultaneously  the  amount  of  flapping  motion  in  the 
torsional  mode  increases  so  that  the  modal  ratios  of  the  two  modes  approach 
each  other.  At  jx  =  .35,  the  two  ratios  become  almost  identical.  A 
similar  effect  occurs  with  the  elastic  axis  at  k$%  chord.  Figure  4l;  how¬ 
ever,  in  this  case  the  amount  of  flapping  displacement  in  the  torsional 
mode  is  greater  in  hover  than  with  the  elastic  axis  at  25%  chord.  It  can 
be  seen  that  the  mode  with  the  smaller  phase  angle  at  hover  is  the  one 
which  suffers  static  divergence. 

The  effect  of  moving  the  elastic  axis  aft  may  be  summarized  by  saying  that 
the  effect  is  to  increase  the  flutter  speed  but  to  cause  a  decrease  in  the 
static  divergence  speed.  An  optimum  location  for  this  blade  occurred  with 
the  elastic  axis  in  the  range  32%  to  4($  chord. 

Cases  2.A1,  2.A2,  2. A3,  and  2.A4  of  Table  XI  consider  four  center-of- 
gravity  positions  in  addition  to  the  nominal  position;  the  results  are 
shown  in  Figure  1+2. 

As  the  center  of  gravity  moves  aft  the  static  divergence  speed  decreases 
uniformly,  corresponding  to  the  fact  that  increasing  the  aerodynamic-- 
center-to-center-of-gravity  spacing  causes  a  corresponding  increase  in  the 
negative  aerodynamic  stiffness  acting  in  the  torsion  mode. 

A  threshold  value  for  the  center-of-gravity  offset  occurs ,  below  which 
flutter  will  not  occur.  Once  the  threshold  value  is  reached,  the  flutter 
speed  decreases  steadily  as  the  center  of  gravity  moves  aft.  However,  the 
sensitivity  to  center-of-gravity  offset  decreases  at  the  same  time  so  that 
the  flutter  curve  flattens  out. 

Cases  3.A1  and  3.A2  of  Table  XI  consider  two  torsional  stiffnesses  in 
addition  to  the  nominal  value,  the  results  are  shown  in  Figure  43. 

Increasing  GJ  raises  the  torsional  divergence  speed  as  expected.  It  also 
has  the  effect  of  raising  the  flutter  speed  until  finally  the  flapping  mode 
appears  to  flutter  above  the  static  divergence  speed.  In  fact,  increasing 
GJ  may  have  eliminated  flapping  flutter  because  for  GJ/(GJ)nom  =  1.328, 
frequency  coalescence  of  the  torsional  and  flapping  modes  occurs  at  /x 
=  .366  yet  flutter  has  not  occurred  by  the  time  static  divergence  sets  in 
at  [X  =  .45.  Sometimes,  however,  flutter  does  not  begin  until  after 
frequency  coincidence  has  occurred.  For  example,  at  GJ/(GJ)nom  = 
frequency  coincidence  of  the  torsional  and  flapping  mode?  occurs  at  fx 
~  .235,  but  coupled  flapping-torsion  flutter  does  not  begin  until  [x  =  .27. 
Also,  the  difference  m  rotor  speed  at  which  frequency  coincidence  and 
static  divergence  occur  decreases  with  increasing  GJ,  as  can  be  seen  from 
Figure  UH.  Hence,  it  is  more  likely  that  the  flutter  sets  in  after  tor¬ 
sional  divergence  has  occurred  rather  than  that  flutter  has  been  eliminated. 

Cases  4.A1  and  4.A2  of  Table  XT  consider  two  blade  flatwise  bending  stiff¬ 
nesses  in  addition  to  the  nominal  value;  the  results  are  shown  in  Figure  45. 
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Increasing  El  has  eliminated  the  flutter  involving  the  first  flatwise 
bending  and  torsional  modes.  However,  the  other  flutter  mode  involving 
mainly  the  rigid-body  flapping  and  torsional  modes  shows  little  change 
with  increasing  El  ,  as  is  to  be  expected.  Hence,  increasing  the 
flatwise  bending  stiffness  will  raise  the  flutter  speed  and  perhaps  will 
eliminate  flutter  modes  involving  elastic  flatwise  bending;  however,  it 
cannot  help  toward  increasiiig  the  flutter  speed  of  flutter  modes  in  which 
flatwise  motion  is  due  principally  to  rigid-body  flapping. 

Cases  5-A1  and  5*A2  of  Table  XI  consider  two  blade  mass  variations  in 
addition  to  the  nominal;  the  results  are  shown  in  Figure  1*6. 

Changing  the  blade  mass  has  comparatively  little  effect  on  the  static 
divergence  speed.  The  main  effect  is  to  cause  a  change  in  the  flutter 
mode.  For  a  blade  mass  distribution  half  the  nominal  value,  the  flutter 
mode  involves  first  flatwise  bending,  while  for  a  blade  mass  distribution 
twice  the  nominal  value,  the  flutter  mode  involves  rigid  body  flapping. 
Based  on  the  results  obtained  in  cases  5-A1  and  5-A2,  it  is  expected  tnat 
the  flutter  speed  of  the  light  blade  (half  the  nominal  mass)  could  be 
raised  substantially  by  increasing  the  blade  bending  stiffness.  By  com¬ 
parison,  the  heavy  blade  should  be  relatively  insensitive  to  changes  in 
the  flatwise  bending  stiffness. 

Cases  6. A1  and  6.A2,  of  Table  XI  consider  two  variations  in  blade  pitching 
inertia,  in  addition  to  the  nominal  value;  the  results  are  shown  in  Figure 

1*7. 

Increasing  or  decreasing  la  from  its  nominal  value  reduces  the  stable 
operating  range  of  the  rotor.  In  fact  for  Ia/(Ia)(goM  <  0-4  and 

la / ( IoJnom  >  2  ,  the  rotor  suffers  flapping  flutter  for  /x  >  0  , 

so  that  there  is  no  stable  operating  region. 

RESULTS  FOR  THE  NONARTICULATED  BLADE 


For  reference,  the  variation  of  the  solution  frequency  and  decay  rate  ver¬ 
sa*  rotor  speed  is  given  in  Figures  1+8  and  1*9,  respectively,  for  the  nom¬ 
inal  blade  operating  in  hover.  The  blade  flutters  at  ll60  RPM,  followed 
by  static  divergence  at  1^00  RPM.  Also  shown  in  Figure  50  is  the  variation 
of  the  critical  advance  ratio  versus  rote:  speed.  As  expected,  the  higher 
the  rotor  speed,  the  lower  the  advance  ratio  at  which  instability  occurs, 
until  finally  the  Dlade  flutters  in  the  hover  condition  for  SI  =  1160 
RPM.  This  is  in  contrast  to  the  Cornell  experimental  flutter  data  in  which 
instability  occurred  at  approximately  constant  rotor  speed  irrespective  of 
the  advance  ratio. 

Cases  1.R1  and  1.R2  of  Table  XI  consider  tvo  elastic  axis  positions  in 
addition  to  the  nominal  position;  the  results  are  shown  in  Figure  51 •  The 
results  are  similar  to  those  of  the  articulated  case,  and  the  comments 
made  for  that  case  apply  also  to  the  hingeless  blade.  However,  because 
the  hingeless  blade  does  not  have  a  flapping  mode,  the  following  equiv¬ 
alence  of  modes  exists: 


93 


Articulated  Rotor 


Hingeless  Rotor 


Flapping 
First  Flatwise 
Second  Flatwise 


First  Flatwise 
Second  Flatwise 
Third  Flatwise 


Thus,  for  example,  the  behavior  of  the  nonarticulated  blade's  second  flat¬ 
wise  bending  mode  should  be  compared  to  the  first  flatwise  bending  mode  of 
the  articulated  blade,  since  the  mode  shapes  and  frequencies  will  tend  to 
be  approximately  the  same.  It  should  be  noted  that  the  results  shown  in 
Figures  39  and  51  are  applicable  to  the  blades  considered,  and  it  should 
not  be  concluded  that  the  results  are  typical  for  all  practical  ranges  of 
the  various  blade  parameters. 

Cases  2.R1,  2.R2,  2.R3,  and  2.R1*  of  Table  XI  consider  four  center-of- 
gravity  positions  in  addition  to  the  nominal  position;  the  results  are 
shown  in  Figure  52.  Again,  the  results  are  very  similar  to  those  of  the 
articulated  case,  and  the  comments  made  for  the  articulated  blade  apply 
also  to  the  hingeless  blade.  Only  one  static  divergence  speed  is  shown 
(namely,  at  13$  center-of-gravity  offset)  because,  for  all  other  center-of- 
gravity  offsets,  the  highest  value  of  p  for  which  solutions  were  ob¬ 
tained  was  below  the  divergence  speed.  If  solutions  had  been  obtained 
for  higher  values  of  p  ,  then  static  divergence  would  have  occurred  in 
every  case.  It  should  be  noted  that  no  flutter  instability  was  found,  for 
the  range  of  p  considered,  in  the  case  of  the  2.5$  center-of-gravity 
offset.  Therefore  only  the  solutions  for  the  remaining  U  cases  appear  in 
Figure  52. 

Cases  3.R1  and  3.R2  of  Table  XI  consider  two  torsional  stiffnesses  in 
addition  to  the  nominal  value;  the  results  are  shown  in  Figure  53.  The 
results  are  similar  to  the  corresponding  results  for  the  articulated  blade. 
The  difference  may  be  attributed  to  the  fact  that  the  torsional  frequency 
does  not  decrease  as  rapidly  with  increase  in  p  for  the  hingeless 
blade  as  it  does  for  the  articulated  blade.  This  result  applies  to  the 
particular  configuration  treated,  and  should  not  be  applied  to  blades  of 
conventional  stiffness . 

Cases  U.R1  and  4.R2  of  Table  XI  consider  two  blade  flatwise  bending  stiff¬ 
nesses  in  addition  to  the  nominal  value;  the  results  are  shown  in  Figure 
5U.  Increasing  the  El  raises  the  frequency  of  the  first  flatwise,  mode 
bringing  it  closer  in  frequency  to  the  torsional  mode.  As  a  result,  fre¬ 
quency  coincidence  of  these  two  modes  occurs  at  increasingly  lower  values 
of  p  as  El'  increases,  which  causes  a  reduction  in  the  flutter  speed. 
There  is  probably  no  simple  explanation  for  the  increase  in  the  torsional 
divergence  speed  with  increase  in  El,  since  torsional  divergence  occurs 
well  above  the  flutter  speed.  Hence,  by  the  time  it  occurs,  the  first 
flatwise  and  torsional  modes  are  so  strongly  coupled  that  one  cannot  in¬ 
terpret  changes  in  torsional  divergence  speed  in  terms  of  singxe-degree-of- 
freedom  results  for  a  torsional  mode. 
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Cases  5*R1  and  $.Rk  of  Table  XI  consider  two  blade  mass  variations  in 
addition  to  the  nominal;  the  results  are  shown  in  Figure  55  •  Changing 
blade  mass  again  causes  an  indicated  switch  in  the  stability  of  modes. 

The  results  are  considerably  different  in  detail  from  the  corresponding 
results  for  the  articulated  blade.  However,  the  practical  result  is  the 
same;  namely,  that  changing  blade  mass  does  not  significantly  alter  the 
stable  operating  r->nge  of  the  blade. 

Cases  6.R1  and  6.R2  of  Table  XI  consider  two  variations  in  blade  pitching 
inertia  in  addition  to  the  nominal  value;  the  results  art  shown  in  Figure 
56.  The  results  are  almost  identical  to  the  corresponding  results  for 
the  articulated  case  except  that,  as  noted  before,  the  second  flatwise 
mode  does  not  flutter  in  the  hingeless  case. 


SIGNIFICANCE  OF  CLASSICAL  FLUTTER  ANALYSIS  AND  RESULTS 


The  method  of  solution  developed  for  this  study  has  proved  to  be  versatile 
and  well  suited  to  design  use.  A  number  of  assumptions  and  approximations 
were  made  when  deriving  the  equations  of  motion.  It  was  considered  that 
the  resulting  simplified  equations  should  still  be  able  to  predict  flutter 
speeds  with  acceptable  accuracy.  It  was  initially  intended  to  check  these 
simplified  equations  by  comparing  the  solutions  with  experimental  flutter 
data  obtained  from  blades  operating  at  high  tip  Mach  numbers  and  advance 
ratios.  As  explained  previously  this  was  not  possible.  The  correlation 
that  was  attempted  proved  encouraging,  but  many  questions  remain  to  be 
answered.  What  is  still  needed  is  to  be  able  to  correlate  with  a  large 
body  of  accurate  experimental  test  data  at  high  advance  ratios  and  tip 
Mach  numbers.  The  estimated  accuracy  of  experimentally  determined  rotating 
blade  flutter  speeds  is  rarely  stated,  if  ever.  Thus,  in  a  correlation 
stuay,  the  experimental  results  have  to  be  assumed  to  be  accurate.  Yet, 
the  possibility  of  inaccuracies  existing  in  the  experimental  results  is 
considerable.  Examples  have  occurred  in  the  past  where,  for  the  simpler 
case  of  two-dimensional  airfoils  oscillating  in  incompressible  flow,  dis¬ 
crepancies  between  measured  and  calculated  flutter  speeds  were  removed  by 
improvements  in  the  experimental  technique.  In  addition,  there  is  the 
possibility  that  a  measured  flutter  speed  may  not  be  due  to  classical 
flutter.  It  would  be  more  satisfactory  if  the  solution  to  the  simplified 
equations  and  the  experimental  determination  cf  flutter  speeds  could  be 
made  part  of  one  test  program. 

The  limited  correlation  that  was  attempted  showed  that  the  calculated 
flutter  speed  was  always  conservative  for  those  cases  in  which  it  was 
reasonably  certain  that  the  flutter  mechanism  was  classical  flutter.  At 
times ,  the  results  were  so  conservative  as  to  be  of  little  use  in  predic¬ 
ting  flutter  speeds.  It  may  be  pointed  out  that  a  considerable  difference 
exists  between  attempting  to  predict  accurate  flutter  speeds  and  a  param¬ 
etric  study  in  which  the  effect  on  the  flutter  speed  due  to  changing 
parameters  is  investigated.  In  the  latter  case,  which  is  usually  the  major 
portion  of  a  flutter  investigation,  the  accuracy  of  the  solutions  is  of 
secondary  importance.  The  main  interest  is  the  detection  and  approximate 
location  of  any  regions  of  instability.  The  accurate  evaluation  of  their 
extent  is  of  secondary  importance.  The  major  requirement  for  parametric 
variations  is  that  the  equations  of  motion  always  give  conservative  flutter 
solutions.  Then  the  method  is  capable  of  locating  all  regions  of  insta¬ 
bility,  with  the  calculated  boundaries  enclosing  the  true  boundaries,  the 
difference  between  the  two  being  a  function  of  the  degree  of  conservatism 
of  the  flutter  solution.  The  possibility  of  missing  an  instability  is  then 
confined  to  incorrect  choice  of  the  parameters.  From  past  experience  and 
with  an  understanding  of  the  flutter  mechanism,  a  great  deal  of  information 
can  be  obtained  from  such  investigations ,  and  it  can  be  applied  toward 
obtaining  a  flutter  clearance.  This  use  of  indirect  methods  was  widespread 
before  the  introduction  of  electronic  computation  aids;  it  will  have  to  be 
relied  on  at  present  for  helicopter  blade  flutter  in  the  high-speed  flight 
condition  until  either  a  more  satisfactory  aerodynamic  theory  is  developed 
or  existing  approximate  theories  can  be  checked  experimentally  over  a  wide 
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range  of  flight  conditions.  Until  these  improvements  become  available,  it 
is  considered  that  the  simplified  equations  of  motion  developed  in  this 
study  represent  an  acceptable  method  for  investigating  the  flutter  charac¬ 
teristics  of  a  helicopter  blade  in  high-speed  flight;  it  is  believed  that 
this  method  is  suited  for  design  use  where  speed  of  solution  is  essential. 

The  parametric  study  should  be  used  with  caution,  as  it  has  limited  appli¬ 
cability.  In  this  respect,  it  may  be  noted  that  the  equations  of  motion 
have  many  types  of  coupling  terms,  and  any  one  may  contribute  critically. 
For  a  particular  class  of  blade  or  a  particular  flight  condition,  certain 
coupling  terms  will  predominate;  if  variations  are  made  within  this  class, 
the  variations  will  have  a  characteristic  form.  However,  for  a  different 
class  of  blade  or  flight  condition,  other  coupling  terms  will  predominate 
and  cause  the  parameter  variations  to  take  on  an  entirely  different  form. 
It  appears  that  the  center-of -gravity  offset  is  one  of  the  few  parameters 
that  are  important  for  most  classes  of  blades.  The  same  situation  applies 
to  fixed-wing  flutter.  Parameter  variations  have  been  published  consider¬ 
ing  straight,  high-aspect-ratio  wings  with  no  lumped  masses,  traveling  at 
low  subsonic  speeds.  However,  these  results  do  not  apply  if  the  wing  is 
traveling  at  high  subsonic  Mach  numbers  or  if  a  flexibly  mounted  store  is 
added  to  the  outboard  portion  of  the  wing.  These  limitations  do  not  de¬ 
crease  the  usefulness  of  parameter  variations  in  general,  but  they  do  mean 
that  a  new  set  of  curves  has  to  be  generated  when  investigating  a  new 
class  of  blade  or  flight  condition. 
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EXTENDED  NORMAL  MODE  TRANSIENT  ANALYSIS  CALCULATIONS 
CORRESPONDING  TO  CLASSICAL  FLUTTER  INVESTIGATIONS 


PURPOSES  AND  METHODS 


The  extended  Normal  Mode  Transient  Analysis  was  empl. yed  to  perform  blade 
response  calculations  corresponding  to  blade  configurations  and  flight 
conditions  for  which  classical  flutter  was  calculated. 

The  purpose  of  these  calculations  was  to  provide  comparative  results  and 
to  obtain  time  histories  of  blade  response  over  a  complete  rotor  revolu¬ 
tion,  part  of  which  was  unstable  from  the  standpoint  of  the  fixed  azimuth 
classical  flutter  calculation. 

It  should  be  noted  that  the  fixed  azimuth  classical  flutter  calculation 
and  the  extended  Normal  Mode  Transient  Analysis  are  quite  different  in 
their  basic  assumptions.  As  pointed  out  in  the  introduction  of  this  vol¬ 
ume,  these  fixed-azimuth  classical  flutter  calculations  consider  pertinent 
inertial  and  elastic  blade  structural  effects,  together  with  unsteady, 
compressible  aerodynamic  loadings.  These  aerodynamic  loadings  include 
phase  lag  effects  which  arise  because  of  blade  motions  in  an  assumed  steady 
velocity  field  appropriate  to  a  particular  azimuth  position.  This  assump¬ 
tion  permits  solution  of  the  flutter  stability  equations  in  a  manner  which 
is  mathematically  similar-  to  the  usual  fixed-wing  flutter  problem,  since 
the  coefficients  of  the  equations  are  constant  in  time.  The  extended 
Normal  Mode  Transient  Analysis,  in  contrast,  integrates  the  blade  differ¬ 
ential  equations  of  motion  in  a  step-by-step  fashion,  considering  pertinent 
inertial  and  structural  effects,  as  well  as  the  time  history  of  instantan¬ 
eous  aerodynamic  loadings,  including  stall  and  compressibility  effects. 
These  loadings,  however,  are  so-called  quasi-steady  loadings.  This  refers 
to  the  de  ermination  of  aerodynamic  loadings  by  steady-state  airfoil  data 
and  by  the  instantaneous  relative  velocity  vector  at  the  75#  chord  point  of 
each  blade  strip.  This  excludes  phase  lag  effects  between  motion  and  aero¬ 
dynamic  loading,  which  are  usually  included  in  flutter  calculations.  As 
pointed  out  earlier,  this  precludes  prediction  of  classical  flutter  in  some 
cases.  However,  the  type  of  classical  flutter  arising  because  of  aft 
center-of-gravity  location  on  a  fixed  wing  can  be  predicted  with  quasi¬ 
steady  aerodynamic  effects.  Hence,  it  could  be  expected  that  the  extended 
Normal  Mode  Transient  Analysis  results  would  reflect  classical  flutter-type 
instabilities  due  to  aft  blade  center-of-gravity  position. 

Therefore,  the  blade  configuration  chosen  was  that  used  for  the  classical 
flutter  parameter  study.  This  blade  had  a  flapping  hinge  only  for  the 
articulated  case,  and  was  assumed  to  be,  rigidly  built  in  for  the  nonartic- 
ulated  case.  Edgewise  motions  were  suppressed  in  both  cases.  In  the 
articulated  case,  rigid-body  flapping,  two  flatwise  bending  modes,  and  one 
torsional  mode  were  considered.  The  natural  frequencies  of  these  modes 
were  1.11,  2.69,  4.8l,  and  5-62  cycles  per  revolution,  respectively.  In 
the  nonarticulated  case,  three  flatwise  bending  modes  and  one  torsional 
mode  were  considered.  The  natural  frequencies  of  these  modes  were  1.43, 
4.12,  8.48,  and  6.14  cyr’es  per  revolution,  respectively.  In  order  to 
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make  the  results  comparable  to  those  for  the  classical  flutter  calcula¬ 
tions,  the  rotor  was  assumed  to  be  very  lightly  loaded,  with  a  collective 
pitch  setting  of  1°.  A  small  loading  of  this  magnitude  was  used  to  insure 
that  some  blade  excitation  was  present,  which  would  make  any  tendency 
toward  instability  more  apparent.  The  shaft  angle  of  attack  for  the  rotor 
was  effectively  zero.  In  common  with  the  classical  flutter  calculation, 
no  cyclic  pitch  inp1"'-  was  considered.  The  normal  full-scale  Reynolds 
number  aerodynamic  data  were  used,  appropriate  to  an  NACA  0012  airfoil. 
Since  blade  angles  of  attack  remained  low,  the  blade  response  Results  of 
interest  were  assumed  to  be  unaffected  by  Reynolds  number  and  airfoil 
camber. 

The  extended  Normal  Mode  Transient  Analysis  was  started  with  the  blade  at 
the  zero  azimuth  position.  Starting  values  for  articulated  blade  flapping 
were  estimated  to  be  steady  flapping  angles  in  each  case.  The  first  flat¬ 
wise  bending  amplitude  was  used  in  the  same  way  for  the  nonarticulated 
blade.  Starting  values  for  all  other  elastic  modal  velocities  and  dis¬ 
placements  were  zero. 

After  starting  in  the  above  manner,  the  numerical  integration  of  the  blade 
equations  of  motion  was  allowed  to  proceed  for  a  maximum  of  10  rotor 
revolutions.  In  many  cases,  the  blade  motions  became  cyclic  with  respect 
to  rotor  azimuth  prior  to  10  revolutions  after  starting.  In  these  cases, 
the  computer  program  automatically  terminated  the  solution,  since  each  of 
the  subsequent  revolutions  would  have  been  identical  if  they  had  been  cal¬ 
culated, 

CONFIGURATIONS  AND  FLIGHT  CONDITIONS  STUDIED 


The  method  described  above  was  employed  for  each  of  the  configuration 
flight-condition  combinations  shown  in  Table  XIII.  Pertinent  character¬ 
istics  of  the  blade  were  otherwise  identical  to  the  nominal  blade  config¬ 
uration  used  in  the  classical  flutter  parameter  studies. 

PRESENTATION  AND  DESCRIPTION  OF  NORMAL  MODE  TRANSIENT  ANALYSIS  RESULTS 


Time  histories  of  flapping  and  ^orsion  modal  displacements  are  shown  in 
Figures  57  and  58.  Results  are  shown  for  various  center-of-gravity 
positions  at  each  of  the  two  advance  ratios.  The  solutions  for  x/cb  =  0, 
-.065,  and ”.0975  became  cyclic  within  a  few  revolutions,  as  shown  in  Figure. 
57  for  x/cb  =“.0975.  All  these  solutions  were  similar,  with  somewhat 
smaller  flapping  and  torsional  amplitudes  for  the  s  nailer  center-of-gravity 
offsets.  The  solution  for  x/cb  =-.130,  shown  in  Figure  57,  did  not 
become  cyclic,  even  after  10  revolutions.  A  slowly  decaying  one-half -per- 
revolvtion  component  can  be  seen  in  the  torsion  modal  time  history,  with 
every  other  revolution  being  similar.  It  appears  that  the  motions  would 
not  grow  significantly  larger  if  the  solution  were  allowed  to  continue. 
Note,  however,  that  torsional  motions  are  very  much  larger  than  for  the 
x/cb  =”.0975  case.  Higher  frequency  components  of  approximately  four 
per  revolution  are  visible.  Figure  58  is  similar  to  Figure  57,  at  an  even 
lower  advance  ratio.  Note  that  the  x/cb  =-.130  case  converges  rapidly 
at  the  lower  advance  ratio.  Setting  the  center-of-gravity  position  further 
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aft,  to  x/cb  =~.l65  and-.20,  results  in  progressively  more  violent 
oscillations  of  a  complicated  nature.  These  apparently  grow  with  time  un¬ 
til  blade  destruction  occurs. 

The  nonarticulated  blade  cases  were  found  to  converge  rapidly,  even  with 
x/cb  =~.130.  The  first  flatwise  response  was  h6%  of  the  articulated 
blade  flapping  response,  although  the  amplitudes  of  the  nonarticulated 
blade  also  increased  rapidly  as  the  center  of  gravity  was  moved  aft. 

,  Figure  59  shows  blade  tip  deflections  during  the  last  revolution  of  the 
various  solutions  for  the  articulated  blade  at  advance  ratio  u.  =  .30. 
Comparison  with  Figure  57  shows  that  most  of  the  vertical  tip  deflection  is 
the  result  of  rigid-blade  flapping  motion.  Note  that  moving  the  center  of 
gravity  aft  results  in  higher  average  blade  twist,  which  acts  like  an  in¬ 
crease  in  collective  pitch. 

In  order  to  present  the  total  elastic  blade  deformations.  Figure  60  was 
drawn  for  the  articulated  blade  with  advance  ratio  fx  =  .30.  Note  that 
the  elastic  blade  deflections  are  essentially  cyclic,  with  a  large  component 
of  twice  rotor  frequency.  These  deformations  are  the  maximum  value  occur- 
ringon  the  blade  and  are  proportional  to  bending  and  twisting  stresses  and 
moments  at  that  section  of  the  blade. 

Figure  6l  is  a  similar  plot  for  the  articulated  blade  with  advance  ratio 
jj.  ~  .15  and  center-cf -gravity  positions  further  aft.  Note  that  for 
x/cb  =-.165  and-. 20,  the  elastic  deformations  are  noncyclic  with  azimuth 
angle  and  that  a  predominate  frequency  of  slightly  less  than  three  times 
rotational  frequency  is  present.  Note  that  an  enormous  increase  in  the 
amplitudes  of  the  elastic  deformations  takes  place  as  the  center  of  gravity 
is  moved  aft.  It  is  probable  that  the  deformations  would  continue  to  grow 
if  the  solution  were  continued  beyond  10  revolutions. 

In  Figure  62,  the  total  elastic  blade  deformations  are  shown  for  the  non¬ 
articulated  blade.  This  blade  was  stiff  enough  in  flatwise  bending  to  make 
its  behavior  distinctly  different  from  the  similar  articulated  blade.  The 
bending  deformation  d2w/dx2  consists  almost  entirely  of  one-per- 

revolution  bending.  The  blade  has  maximum  torsional  deflection  at  the  ad¬ 
vancing  blade  position.  Note  that  this  deflection  could  be  considered  an 
incipient  torsional  divergence  due  to  the  aft  center-of-gravity  position. 

In  Figures  63  and  64,  the  one-half  peak-to-peak  nondimensional  blade  defor¬ 
mations,  which  are  proportional  to  stress  or  elastic  moments,  are  plotted 
against  center-of-gravity  position  for  the  articulated  blade.  Figure  65  is 
a  similar  plot  for  the  nonarticulated  blade.  Note  that  a  rapid  rise  in 
peak-to-peak  deformations  occurs  for  the  articulated  configurations  as  the 
center  of  gravity  is  moved  aft  past  certain  limits.  The  corresponding  rise 
for  the  nonarticulated  blade  is  similar  but  much  less  pronounced. 
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DISCUSSION  OF  NORMAL  MODE  TRANSIENT  ANALYSIS  RESULTS 
CORRESPONDING  TO  CLASSICAL  FLUTTER  INVESTIGATIONS 


In  this  subsection,  appropriate  comparisons  will  be  made  of  the  extended 
Normal  Mode  Transient  Analysis  results  and  the  corresponding  fixed-azimuth 
classical  flutter  results. 

As  pointed  out  previously,  the  Normal  Mode  Transient  Analysis  and  the  fixed 
azimuth  flutter  calculation  are  quite  different  in  their  basic  intent  and 
analytical  approach.  In  spite  of  this,  some  comparisons  of  the  two  sets 
of  results  should  be  made.  The  aft  center-of-gravity  practical  operating 
boundaries  predicted  for  the  limited  number  of  cases  considered  here  do 
show  encouraging  agreement  between  the  two  methods. 

The  most  important  direct  comparison  is  thus  of  Figures  42,  63  and  64  for 
the  articulated  blade,  and  of  Figures  52  and  65  for  the  nonarticulated 
blade.  As  is  shown  on  Figure  42,  flatwise-bending-torsion  flutter  occurs 
at  an  advance  }jl  =  .30  for  a  center-of-gravity  offset  of  11$,  followed 
by  flapping-torsion  flutter  at  a  center-of-gravity  offset  of  about  12.2$. 
Figure  63  shows  that  the  flatwise  and  torsion  stresses  increase  rapidly  as 
the  center  of  gravity  is  located  aft  of  these  two  positions.  Reference  to 
Figures  57,  58,  and  59,  however,  shows  that  the  flapping  motion  of  the 
blades  increases  comparatively  little.  Inspection  of  Figure  60  shows  fre¬ 
quency  components  of  the  order  expected  for  the  flatwise-torsion  flutter. 
Thus,  some  manifestations  of  the  flatwise-torsion  flutter  predicted  for  the 
idealized  fixed-azimuth  system  appear  in  the  Normal  Mode  Transient  Analysis 
No  particular  evidence  of  a  flapping-torsion  flutter  can  be  discerned.  It 
would  be  expected  that  the  fixed-azimuth  assumption  would  be  more  valid  for 
the  higher  frequency  flutter  modes.  Similar  observations  can  be  made  for 
an  advance  ratio  of  jj.  =  .15,  for  which  flatwise-torsion  flutter  is  pre¬ 
dicted  at  a  center-of-gravity  offset  of  15$  and  flapping-torsion  flutter  at 
an  offset  of  17.5$.  The  flutter  boundary  for  the  nonarticulated  blade 
shown  in  Figure  52  corresponds  to  first  flatwise-torsion  flutter,  which 
corresponds  to  flapwise-torsion  flutter  for  the  articulated  blade.  The 
onset  of  this  comparatively  low  frequency  flutter  for  the  fixed-azimuth 
system  is  reflected  in  a  fairly  gentle  rise  in  stress  as  center-of-gravity 
offset  increases. 

On  the  basis  of  the  above  information,  it  can  be  expected  that  the  fixed- 
azimuth  flutter  calculations  will  be  most  accurate  for  classical  flutter 
frequencies  greater  than  two  times  rotor  frequency  when  the  advance  ratio 
is  greater  than  p.  =  .15.  When  flutter  frequencies  are  smaller,  the  pre¬ 
dictions  of  flutter  speed  are  expected  to  be  increasingly  conservative. 


Figure  3.  A  Typical  V-g  Plot. 


Figure  1j.  A  Typical  Variation  of  Decay  Rate  C/C 


Notation  for  Blade  Displacements. 
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Figure  7-  Radial  Displacements  Resulting  From  Edgewise  Bending. 


Figure  8.  Radial  Displacements  Resulting  from  Flatwise  Bending. 
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9.  Air  Velocity  Relative  to  the  Blade. 
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Figure  il.  Arbitrary  Chordvise  Displacement  of  the  Mass 

Balance  Center  of  Gravity  and  the  Elastic  Axis. 


Figure  12.  Arbitrary  Chordwise  Displacement  of 
vhe  Aerodynamic  Center. 
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Figure  15  •  Notation  Used  for  a  Two-Dimensional  Chordvise 
Rigid  Airfoil  in  Direct  Flow  and  Reverse  Flow. 
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Figure  16.  Example  of  the  Use  of  the  Frequency  Versus 
Decay-Rate  Digital  Plot. 
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Figure  17-  Example  of  the  Digital  Plot  of  Decay  Rate  and 
Frequency  Versus  Forward  Velocity. 


FREQUENCY  ,  «  ,  CPM 


Figure  19-  Comparison  of  Sikorsky  Flutter  Analysis  Results  for 
Flutter  Frequency  with  the  Solution  of  Reference  12 
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Figure  20.  Comparison  of  Sikorsky  Flutter  Analysis  Results  for 
Decay  Rate  With  the  Solution  of  Reference  12. 
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Figure  21.  Comparison  of  the  Sikorsky  Flutter  Analysis  Results 
for  Flutter  Frequency  with  the  Solution  of  Reference 
12  Using  Simplified  Aerodynamics. 
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Figure  22.  Comparison  of  the  Sikorsky  Flutter  Analysis  Results 
for  Decay  Rate  with  the  Solution  of  Reference  12 
Using  Simplified  Aerodynamics. 
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Figure  23.  Comparison  of  Calculated  and  Measured  Flutter 
Speeds  for  the  Model  Helicopter  Blade 
NACA  23012  l(R). 
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Figure  2k.  Comparison  of  Calculated  and.;  Measured  Flutter 
Speeds  for  the  Model  Blade  MCA  23012  2(R)  , 


FREQUENCY,  w  ,  CPM 


Figure  25.  Variation  of  Solution  Frequency  With  Change  in 
Rotor  Speed  for  the  Model  Blade  NACA  23012  2(R) ; 
Tip  Mach  Number  <  0,2. 
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Figure  2 6.  Variation  of  Solution  Decay  Rate  With  Change  in 
Rotor  Speed  for  the  Model  Blade  NACA  23012  2(R) ; 
Tip  Mach  Number  <  0.2. 
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Figure  27.  Comparison  of  Calculated  and  Measured  Flutter 
Speeds  for  the  Model  Blade  NACA  23Q12  2(.b  \ . 
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FLUTTER  FREQUENCY  RATIO 
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INVERSE  FLUTTER  SPEED  RATIO 
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Figure  28.  Comparison  of  Calculated  and  Measured  Flutter  Speeds 
for  the  CAL  Blade  Model  1* ;  7g/C  =.139  ,  S5^1/u5g0  -  .63  . 
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Figure  29.  Variation  of  Solution  Frequency  With  Change  in  Rotor 
f.peed  for  the  CAL  Blade  Model  If;  Xg/C  =  .139  , 

=  .63  ,  fj,  -  0  . 
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Figure  30 r  Variation  of  Solution  Decay  Rate  With  Change  in  Rotor 
Speedjfor  the  CAL  Blade  Model  4;  JE/C  =  .139  , 
Uft/UBo  =  -63  ,  fJ.  =  0  . 
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of  Calculated  and  Measured  Flutter  Speeds 
,  Blade  Model  k\  xE/C  =  .139  ,  5ty,/5Jeo  =  1.3!  . 
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Figure  32.  Variation  of  Solution  Frequency  With  Change  in  Rotor 
Speed  for  the  CAL  Blade  Model  k;  Xg/C  =  .139  , 

<*><pi/&60  =  -63  »  H-  -  -5  • 
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Figure  33.  Variation  of  Solution  Decay  Rate  With  Change  in  Rotor 
Speedy  for  the  CAL  Blade  Model  4;  Xc/C  -  .139 
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Figure  34.  Comparison  of  Calculated  and  Measured  Flutter 
Speeds  for  the  CAL  Blade  Model  4;  Xe/C  =  .044, 
®<#>i/®60  =  1«°5 
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Figure  36.  Comparison  of  Calculated  and  Measured  Flutter 
Speeds  for  the  CAL  Blade  Model  4;  Xg/C  =  .0036 , 
=  L44  • 
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FREQUENCY,  w  ,  CPM 
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Figure  37.  Variation  of  Solution  Frequency  With  Change 
in  Rotor  Speed  for  the  Nominal  Articulated 
Blade;  fi  =  0. 
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Figure  38.  Variation  of  Solution  Decay  Rate  With  Change 
in  Rotor  Speed  for  the  Nominal  Articulated 
Blade;  fJ.  =0. 
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ELASTIC  AXIS  POSITION  , 

(  PERCENT  OF  CHORD  AFT  OF  LEADING  EDGE  ) 


Figure  39*  Varying  the  Chordwise  Position  of  the  Elastic 
Axis  for  an  Articulated  Blade. 


135 


00 


504 


o  20 
%  10 


E 

l«- 


2 

5 

a: 

-i 

< 

Q 

O 

2 

a 

UJ 

N 


< 

2 

a: 

o 

z 


1.0 

.8 

.6 

.4 

.2 

0 


FLAPPING  MODE 
(  MODE  I  ) 


TORSION"  MODE 
(MODE  4) 


N 

\ 

T( 

)RSIO 

/ERGE 

) 

NAL  . 

:nce 

DP 

0  .2  .4 

ADVANCE  RATIO,  p 


Figure  1*0.  Variation  of  Modal  Amplitude  Ratio  and  Phase 
Angle  With  Change  in  Advance  Ratio  for  the 
Nominal  Articulated  Blade. 
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Figure  42.  Varying  the  Chordwise  Position  of  the  Center 
of  Gravity  for  an  Articulated  Blade. 


ADVANCE  RATIO 


Figure  43.  Varying  the  Blade  Torsional  Stiffness  for  an 
Articulated  Blade;  GJN0M  =  12,436  Lb-In?. 
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Figure  44.  Typical  Variation  of  Torsion  Frequency  With 
Change  in  Advance  Ratio  for  Large  Center-of- 
Gravity  Offset. 


Figure  45.  Varying  the  Blade  Flatwise  Bending  Stiffness  for 
an  Articulated  Blade;  EIznqm  =  29,500  Lh-In.  . 
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ED  OF  FLAPPING  MODE 


Figure  U$>.  Variation  of  Solution  Decay  Rate  With  Change 
in  Rotor  Speed  for  the  Nominal  Nonarticulated 
Blade;  fj.  =0. 
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Figure  50.  Variation  of  Critical  Advance  Ratio  With  Change 
in  Rotor  Speed  for  the  Nominal  Nonarticulated 
Blade . 
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ELASTIC  AXIS  POSITION  , 

(PERCENT  OF  CHORD  AFT  OF  LEADING  EDGE) 


Figure  51*  Varying  the  C'nordwise  Position  of  the  Elastic 
Axis  for  a  Honarticulated  Blade. 
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Figure  52.  Varying  the  Chordvise  Position  of  the  Center 
of  Gravity  for  a  Nonarticulated  Blade. 
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Figure  53. 


Varying  the  Blade  Torsional  Stiffness  for  a 
Nonarticulated  Blade;  GJjjqm  =  12,  ^36  Lb-In?. 
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Figure  5k.  Varying  the  Blade  Flatwise  Bending  Stiffness  for  a 
Nonarticulated  Blade:  EIZnom  =  29,500  Lb-In.2. 
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Figure  55*  Varying  the  Blade  Mass  for  a  Nonarticulated 
Blade;  m  jjqm  =  *0162  Slug  /Ft. 
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Figure  56.  Varying  the  Blade  Pitching  Inertia  for  a  Non- 

articulated  Blade;  Ia  N0M  =  ,7l»2  x  10"1*  Slug-Ft2/Ft. 
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Figure  58.  Continued. 
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Figure  60.  Articulated  Blade  Elastic  Deformations 

Azimuth  Angle  With  Aft  Center-of -Gravity 
Position;  tt=.30,  A|S=0°  ,  B1S  =  0°  , 


Figure  60.  Concluded. 
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Figure  6l.  Articulated  Blade  Elastic  Deformations  Versus 
Azimuth  Angle  With  Aft  Center-of-Gravity 
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Nonarticulated  Blade  Deformation  Versus  Azimuth 
Angle  With  Aft  Center-of-Gravity  Position; 

/i=  -30  ,  AIS=  0°  ,  BIS  =  0°  ,  0.75R  =  1.0°  , 
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Figure  63.  One-Half  Peak-to-Peak  Articulated  Blade 

Deformations  Versus  Aft  Center-of-Gravity 
Position;  u  =  .30,  A|S=0°,  Bis=0°  , 
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Figure  64.  One-Half  Peak-to-Peak  Articulated  Blade 

Deformations  Versus  Aft  Center-of-Gravity 
Position;  u.  =  .15  ,  Ais=0°  ,  Bis  =  0°  , 
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Figure  65. 


)ne-Half  Peak-to-Peak  Nonart iculated  Blade 
Deformations  Versus  Aft  Center-of' -Gravity- 
Position;  fj.  =  -30  ,  Ais"  0°  ,  Bis  =  0°  , 
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TABLE  I.  EXAMPLE  OF  THE  APPLICATION  OF  REVERSE-FLOW 
RELATIONS 


1-1N  k  Reference 


Tabulated 


Calculated  Using 
Reverse-Flow  Relations 


0 

0.2 

2 

.72758 

.017730 

.72756 

.017723 

7 

.72750 

.017750 

.7361,7 

.017720 

9 

.72760 

.017720 

.72760 

.017720 

o.5 

0.2 

2 

.761,26 

.030701 

.761,13 

.030723 

It 

.761,25 

.030700 

.761,29 

.030682 

0 

7 

.76375 

.029850 

.77862 

.029860 

0.8 

0.1 

2 

1.00268 

.038675 

1.002Q9 

.038706 

7 

1.01550 

.039200 

1.01,030 

.039720 

0.8 

0.7 

2 

.69096 

.068596 

.6901,1 

.068628 

7 

.63607 

.059300 

.66570 

.059325 

1.0 

0.2 

2 

.68508 

.091657 

.68510 

.0?i61|8 

3 

.68513 

.09161,8 

.68511 

.091650 

l.n 

0.1,75 

i, 

.52501 

.090162 

.52501 

.090161 

8 

.52501 

.090163 

.52501 

.090162 

TABLE  II.  EXAMPLE  OF  ERRORS  IN  THE  TABULATED 
COEFFICIENTS 

Mn  -0.8 


(  REF.  2  ) 


.58728 


.9U856 

.99870 

.96130 


(  REF.  7  ) 


1.6250 


1.U928 


1.16787 

1.07725 
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TABLE  III.  AIRFOIL  DATA  USED  FOR  THE  TEST  CASE 


Item 

Symbol 

Value 

Units 

Blade  Chord 

2b 

1.0 

ft 

Blade  Mass  per  Unit  Span 

m 

. 16786 

slug/ft 

Radius  of  Gyration  About  Pitch  Axis 

ra 

.28 

ft 

Aerodynamic  Center  Chordwise 
Position  Aft  of  Leading  Edge 

AC 

25. 

%  chord 

Elastic  Axis  Chordwise  Position 

Aft  of  Leading  Edge 

EA 

30. 

%  chord 

Center-oi-Gravity  Chordwise 
Position  Aft  of  Leading  Edge 

CG 

42.5 

%  chord 

Vertical  Translation  Natural 
Frequency 

wh, 

218. 

cpm 

Pitching  Natural  Frequency 

a>d| 

1090 

cpm 

Lift  Curve  Slope 

CL  a 

2  t r 

rad"1 

Air  Density 

P 

.  002378  slug/ft3 

Rotational  Speed 

SI 

0.0 

rad/sec 

Azimuth  Angle 

* 

90. 

deg 
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TABLE  VI.  NACA  MODEL  BLADE  PROPERTIES 


First  Flatwise  Bending  Frequency 
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First  Cantilever  Torsion  Mode 
Structural  Viscous  Damping  Coefficient 


TABLE  VII.  CORNELL  AERONAUTICAL  LABORATORY  MODEL  BLADE  PROPERTIES 


Item 

Symbol 

Value 

Unit 

Blade  Number  (  Reference  ll)  ) 

— 

h 

— 

Chord 

2b 

3.5 

in. 

Flapping  Hinge  Radius 

e 

0 

in. 

Outboard  End  of  Rigid  Section 

7 

in. 

Rotor  Radius 

R 

1)8 

in. 

Shaft  Axis  Chordvd.se  Position 

25 

%  chord 

Feathering  Axis  Chordvri.se  Position 

25 

%  chord 

Center -of -Gravity  Chordvri.se 

Position  Aft  of  Leading  Edge 

CG 

1)2.5 

%  chord 

Elastic  Axis  Position  Chordwise 
Position  Aft  of  Leading  Edge 

EA 

1)1.0 

%  chord 

First  Flatwise  Bending  Frequency 
(Calculated) 

wh, 

76.9 

rad/sec 

Second  Flatwise  Bending  Frequency 
(Calculated) 

wh2 

92k 

rad/sec 

First  Cantilever  Torsion  Frequency 
(Calculated) 

wat 

6bo 

rad/sec 
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TABLE  VII.  -  Concluded 


Item 


Symbol  Value  Unit 


Rigid  Body  Pitch  Frequency 

we0 

h 6. 9-132.5  rad/sec 

Blade  Mass  Density  Ratio 

H-m 

h6.6 

Flanping  Mechanical  Viscous 

DsmDing  Coefficient  gh  0 

"o 

Rigid  Body  Pitching  Mode 

Structural  Viscous  Damping  Coefficient  gaQ  0 

First  Flatwise  Bending  Mode 

Structural  Viscous  Damping  Coefficient  g^  .Oh 

Second  Flatwise  Bending  Mode 

Structural  Viscous  Damping  Coefficient  p  .Oh 

'h2 

First  Cantilever  Torsion  Mode 

Structural  Viscous  Damping  Coefficient  g„  .Oh 
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TABLE  VIII.  CASES  USED  FOR  CORRELATION  WITH  CORNELL  DATA 


Reference  16 
Figure  No. 

U80 

RAO/SEC 

Cb 

ape0 

V aw  oe0 e0 

10(a) 

1.03 

65.9 

.0036 

-.00317 

10(a) 

1.1*1* 

1*6.9 

.0036 

-.00317 

10(b) 

1.05 

69.0 

.OlO* 

-.101 

10(c) 

.63 

132.5 

.139 

-.338 

10(c) 

1.31 

61*.  5 

.139 

-.338 

TABLE  IX. 

EFFECT  ON  FLUTTER  SPEED  OF  VARYING  MASS 

BALANCE 

a0Bo 

16  r~z — ■ — 

s80 

H- 

a/ud0 

Figure  No. 

v  Q0/3  %B0 

MEASURED  j  E 

4$0  SOLUTION 

10(c) 

-.338 

.63 

0.0 

.25 

.27 

10(c) 

-.338 

1.31 

0.0 

.33 

.32 

10(b) 

-.101 

1.05 

0.0 

.56 

10(b) 

-.101 

1.05 

0.2 

mm 

.29 

10(a) 

-.00317 

1.03 

0.0 

1.0 

.35 

10(a) 

-.00317 

1.10* 

0.0 

1.0 

- 

10(a) 

-.00317 

1.10* 

0.2 

.33 

TABLE  X.  BASIC  PFOPERTIES  OF  THE  BLADE  USED  FOR 
THE  PARAMETER  VARIATION  STUDY 


Item 

- - - -  -  -  -  .  0L 

Symbol 

Value 

Unit 

Chord 

2b 

k 

in. 

Articulated  Blade  Flapping  Hinge  Radius 

e 

2.5 

in. 

Outboard  End  of  Rigid  Section 

8 

in. 

Rotor  Radius 

R 

h6 

in. 

Articulated  Blade  Rotational  Speed 

a 

95 0 

rpm 

Nonart iculated  Blade  Rotational  Speed 

a 

870 

rpm 

TABLE  XT.  CASES  CONSIDERED  IN  THE  VARIATION  OP  PARAMETERS 
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C\  CN 

cvj  cvj 

8888 

XAlAlA^A 
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88888888 

1AlAV\lA"lA\AtA1A 
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TABLE  XIII. 

BLADE  CONFIGURATIONS  AND  FLIGHT  CONDITIONS 
FOR  EXTENDED  NORMAL  MODE  ANALYSIS 
CALCULATIONS  CORRESPONDING  TO  CLASSICAL 
FLUTTER  INVESTIGATIONS 

Case 

No. 

Articulated 

Or 

Nonarticulated 

X 

% 

1 

Articulated 

.30 

0.0 

2 

Articulated 

.30 

-.06^ 

3 

Articulated 

.30 

-.0975 

U 

Articulated 

.30 

-.130 

5 

Art-!  culated 

.1? 

0.0 

6 

Articulated 

.15 

-.130 

7 

Articulated 

.15 

-.165 

8 

Articulated 

.15 

-.200 

9 

Nonart i culated 

.30 

0.0 

10 

Nonarticulated 

.30 

-.0975 

-  11 

Nonarticulated 

.30 

-.130 

CONCLUSIONS 


The  simplified  equations  of  motion  predicted  flutter  speeds  which 
were  conservative  when  compared  with  experimental  flutter  data. 

Torsional  divergence  and  classical  bending-torsion  flutter  of  a 
rotating  blade  are  closely  related,  with  static  divergence  occxirring 
usually  at  a  slightly  higher  speed  than  flutter. 

As  the  blade  tip  Mach  number  approaches  unity,  the  rearward  shift 
of  the  aerodynamic  center  has  a  stabilizing  effect.  The  increase 
in  the  flutter  or  divergence  speeds  was  considerably  greater  with 
the  experimental  data  than  that  which  occurred  with  the  calculated 
speeds.  This  is  attributed  to  the  rearward  shift  of  the  aerodynamic 
centers  being  greater  in  practice  than  is  implied  by  the  aerodynamic 
coefficients. 
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RECOMMENDATIONS 


ft  V 

The  correlation  study  should  he  extended  to  help  define  in  more 
detail  any  limitations  that  exist  in  the  simplified  equations 
of  motion. 

A  program  should  he  undertaken  which  consists,  of  experimental 
measurement  of  rotor  blade-  classical  flutter  speeds  at  high  ad¬ 
vance  ratios  and  tip  Mach  numbers ,  together  with  the  calculation 
of  the  corresponding  flutter  speeds  through  the  use  of  the  sim¬ 
plified  equations  of  motion. 

Aerodynamic  coefficients  taken  from  Reference  15  should  he  used 
in  place  of  the  table  of  compressible  aerodynamic  coefficients 
when  the  blade  tip  Mach  number  and  advance  ratio  are  low. 

The  effect  on  the  aerodynamic  stiffness  and  damping  of  bending 
modes  of  the  spanwise  air  flow  along  the  blade  should  be  invest¬ 
igated. 

The  effects  of  blade  twist  and  coning  angle  should  be  incorporated 
into  the  blade  structural  coefficients. 

The  conservatism  inherent  in  the  fixed  azimuth  approach  should 
be  investigated. 

As  the  next  step  toward  developing  a  more  accurate  unsteady  aero¬ 
dynamic  theory  applicable  to  high  Mach  numbers  and  high  advance 
ratios,  an  attempt  should  be  made  to  modify  the  existing  two- 
dimensional  fixed-wing  unsteady  compressible  flow  theory  to  in-  . 
-elude  in  an  approximate  fashion  the  effect  of  a  nonplanar  wake, 
the  effect  on  the  circulation  of  the  presence  of  more  than  one 
wake,  and  the  effect  of  a  time-varying  dynamic  pressure. 

The  experimental  occurrence  of  flutter  phenomena,  which  appear  to 
be  a  combination  of  stall  flutter  and  classical  flutter  has  been 
noted.  This  type  of  combined  flutter  should  be  studied  by  analysis 
and  test. 
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